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ABSTRACT 


Heat  transfer  tests  were  carrisd  out  using  a  small  tube  bundle  of 
Turbo-B  tubes  in  a  pool  of  different  R-114/oil  mixtures.  By  accurately 
instrumenting  five  tubes  within  the  bundle,  both  the  convective  and 
nucleate  boiling  regions  were  studied  in  detail,  with  emphasis  on  the 
'bundle  effect'  (ie.  the  effect  of  the  lower  tubes  in  operation  on  the 
upper  tubes  within  the  bundle) .  In  addition,  the  influence  of  increased 
amounts  of  oil  on  the  tube  bundle  was  studied  to  see  how  this  affected  the 
overall  heat  transfer  and  in  particular,  the  shape  of  the  hysteresis  loop. 
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I.  INTRODUCTION 


A.  BACKGROUND 

One  of  today's  major  environmental  concerns  is  the  depletion  of  the 
earth's  protective  ozone  layer.  In  1987,  an  international  conference  was 
held  in  Montreal,  Canada  to  address  the  problems  caused  by 
Chlorofluorocarbons  (CFCs)  to  the  earth's  ozone  layer.  CFCs  are  manmade 
chemicals  of  chlorine,  fluorine,  and  carbon  and  are  unique  in  that  they 
have  a  combination  of  desirable  properties:  low  in  toxicity,  non¬ 
flammable,  non-corrosive,  non-explosive,  extremely  stable  and  compatible 
with  many  other  materials.  This  extreme  stability  is  what  causes  problems 
to  the  ozone  layer  due  to  the  fact  that  CFCs  only  break  down  in  the  upper 
atmosphere  when  subjected  to  intense  ultraviolet  radiation.  This  break 
down  produces  chlorine  which  has  been  linked  to  the  depletion  of  the 
earth's  ozone  layer. 

In  September  1987,  24  nations  representing  the  United  Nations 
Environment  Program  (UNEP)  met  and  signed  the  Montreal  Protocol.  They 
discussed  the  substances  that  deplete  the  ozone  layer  [Ref.l]  and  called 
for  a  near-term  freeze  on  the  production  and  consumption  of  these 
substances.  The  agreement  required  production  of  these  chemicals  to  be 
cut  back  to  1986  levels  followed  by  a  two-phased  reduction  culminating  in 
cutbacks  of  50%  by  mid-1998;  this  came  into  effect  on  July  1,  1989.  In 
1990,  a  progress  meeting  was  held  in  London  where  UNEP  delegates  agreed  to 
completely  phase  out  all  CFCs  by  the  year  2000  [Ref.  2].  In  the  spring  of 


1 


1992,  President  Bush  pushed  up  the  complete  phase  out  od  CFC's  by  the  year 
1995. 

The  U.S.  Navy  uses  a  number  of  different  CFCs  (designated  by  "R"  for 
refigerants)  for  various  refrigeration  and  air  conditioning  (AC)  needs. 
Presently,  the  U.S.  Navy  has  approximately  1850  shipboard  AC  plants  using 
both  R-12  (in  reciprocating  compressor)  and  R-114  (in  centrifugal 
compressor)  plants.  To  comply  with  the  Montreal  Protocol  and  U.S. 
legislation,  the  Mechanical  Systems  Branch  (Code  2772)  at  the  Naval 
Surface  Warfare  Center  (NSWC)  is  pursuing  research  in  the  elimination  of 
shipboard  use  of  CFCs.  As  mentioned  by  Chilman  [Ref.  3]  this  research  is 
to  be  completed  in  three  phases: 

1.  To  identify  in  the  short  term  suitable  alternative  ozone-safe 
chemicals  to  replace  R-114  and  R-12.  To  accomplish  this  task,  the  heat 
transfer  characteristics  must  be  similar  to  the  existing  refrigerants  in 
place  and  hence  the  need  for  a  database  exists  for  current  refrigerants 
(R-114  and  R-12)  so  that  they  can  be  compared  to  the  new  proposed 
refrigerants  (HFC-124  and  HFC-134A  respectively) . 

2.  In  the  longer  term,  to  research,  develop,  and  test  substitute 
chemicals  and  alternative  technologies  to  replace  existing  CFC  uses . 

3.  To  implement  new  cooling  system  technologies  into  the  fleet  which 
do  not  depend  on  CFCs  or  their  replacement. 

This  thesis  is  a  continuation  of  the  previous  work  at  NPS  and 
supplements  NSWC's  research  on  alternatives  to  CFCs  by  establishing 
baseline  nucleate  pool  boiling  data  of  pure  R-114  and  R-114/oil  mixtures 
from  a  small  bundle  of  enhanced  tubes  representing  a  section  of  a  flooded 
evaporator.  Emphasis  is  placed  on  the  natural  convection  and  boiling 
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regions,  hystersis  phenomena,  and  analysis  of  various  oil  concentrations 
on  the  overall  heat  transfer  performance. 

B.  OBJECTIVES 

The  objectives  of  this  thesis  are  as  follows: 

1.  Understand  in  greater  detail  both  the  convection  and  nucleate  pool 
boiling  phenomena  and  hysteresis  effects  within  a  small  Turbo-B  tube 
bundle . 

2.  Obtain  data  using  a  Turbo-B  tube  bundle  for  increasing  and 
decreasing  heat  flux  for  R-114/oil  mixture  with  oil  concentrations  of  0, 
1 ,  2 ,  3 ,  6 ,  and  10  percent . 

3.  Compare  data  with  earlier  studies  at  the  Naval  Postgraduate  School 
using  R-114/oil  mixtures  for  other  enhanced  tube  surfaces. 
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II.  LITERATURE  SURVEY 


A.  SMOOTH  TUBE  BUNDLES 

In  recent  years  significant  progress  has  been  made  in  understanding 
nucleate  boiling  heat  transfer  phenomena  on  the  shell  side  of  flooded 
evaporators.  Extensive  work  on  smooth  tube  bundles  has  been  reported  by 
Cornwell  (Leong  and  Cornwell  [Ref.  4],  Cornwell  et  al.  [Ref.  5],  Cornwell 
and  Scoones  [Ref.  6],  Cornwell  [Ref.  7]).  Cornwell  and  Schuller  [Ref.  8] 
conducted  a  photographic  study  of  boiling  R-113  in  a  smooth  tube  bundle  at 
one  atmosphere.  One  of  their  conclusions  was  that  bubbles  leaving  the 
lower  tubes  in  the  bundle  impacted  and  caused  a  sliding  motion  around  the 
upper  tubes.  Cornwell  and  Schuller  observed  the  two-phase  flow  patterns 
and  deduced  that  sliding  bubbles  from  lower  tubes  on  upper  tubes  could 
account  for  significant  heat  transfer  in  the  top  part  of  the  bundle. 
Cornwell  [Ref.  7]  later  found  that  in  the  nucleate  boiling  region,  sliding 
bubbles  and  liquid  forced  convection  could  account  for  all  the  heat 
transfer  in  the  top  of  the  bundle. 

The  influence  of  tube  position  within  a  bundle  of  smooth  tubes  using 
R— 11 ,  R-12,  R-22  and  R-113  has  been  studied  extensively  by  Wallner  [Ref. 
8],  Fujita  et  al.  [Ref.  9],  Chan  and  Shourki  [Ref.  10],  Rebrov  et  al. 
[Ref.  11],  and  Marto  and  Anderson  [Ref.  12].  Using  both  in-line  and 
staggered  tube  arrangements  with  various  tube  pitch-to-diameter  ratios 
between  1.2  and  2.0,  their  work  verified  that  the  influence  of  the  lower 
tubes  in  a  bundle  can  significantly  increase  the  heat  transfer  performance 
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of  upper  tubes  at  low  heat  fluxes  due  to  two-phase  convective  effects.  At 
high  heat  fluxes  (typically  >  50  kW/m2)  in  the  fully  developed  nucleate 
boiling  region,  the  data  for  all  the  tubes  merged  onto  a  single  curve. 
This  is  representative  of  a  single  smooth  tube  and  shows  that  there  is  no 
'bundle  effect'  (ie.  no  improvement  over  a  single  tube  under  similar 
conditions)  in  the  high  heat  flux  (nucleate  boiling)  region. 

Chan  and  Shoukri  [Ref.  10]  studied  the  boiling  characteristics  of  a 
smooth  in-line  tube  bundle  in  R-113.  They  concluded  that  at  lower  heat 
fluxes,  the  heat  transfer  process  is  strongly  influenced  by  two-phase 
convection  effects,  resulting  in  higher  heat  transfer  coefficients  on  the 
upper  tubes.  At  high  heat  fluxes,  however,  they  found  that  the  dominant 
mode  of  heat  transfer  was  nucleate  boiling  from  the  upper  tubes  and  that 
convective  effects  from  below  were  insignificant.  At  these  high  fluxes, 
the  bundle  performance  was  similar  to  the  trends  of  a  single  tube  in  a 
single  tube  apparatus.  Fujita  et  al.  [Ref.  9]  also  found  that  the  heat 
transfer  at  low  heat  fluxes  using  a  smooth  tube  bundle  in  R-113  was 
enhanced  by  convection  induced  by  rising  bubbles  (ie.  a  steady  increase  in 
performance  of  the  upper  tubes  as  additional  lower  tubes  were  activated) . 
They  attributed  this  enhancement  to  the  "positive  bundle  effect".  A  high 
heat  fluxes,  this  enhancement  disappeared. 

Anderson  [Ref.  13]  found  similar  effects  as  above  for  a  smooth  tube 
bundle  in  pure  R-114.  Furthermore,  he  reported  that  the  presence  of  up  to 
3%  oil  (by  mass)  actually  improved  the  heat  transfer  performance.  This  is 
similar  to  data  reported  for  a  single  smooth  tube  by  Vanniarachchi  et  al. 
[Ref.  14].  Furthermore,  at  an  oil  concentration  of  10%,  only  a  slight 
degradation  in  the  heat  transfer  (compared  to  pure  R-114)  was  found.  He 
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obtained  a  maximum  heat  transfer  performance  for  the  bundle  at  an  oil 
concentration  of  2%. 

B.  ENHANCED  TUBE  BUNDLES 

Much  less  work  has  been  done  on  enhanced  tubes  (Enhanced  means  any 
surface  that  is  not  smooth).  However,  the  effects  at  high  and  low  heat 
fluxes  mentioned  above  are  similar  to  those  obtained  for  finned  tube 
bundles  by  Yilmaz  and  Palen  [Ref.  15],  Muller  [Ref.  16],  and  Hahne  and 
Muller  [Ref.  17].  Stephen  and  Mitrovic  (Ref.  18]  looked  at  R-12  and  R-114 
boiling  from  a  GEWA-T  tube  bundle .  Apart  from  the  magnitude  of  the  heat 
transfer  coefficient  varying  with  fluid,  the  trends  are  very  similar  to 
those  mentioned  above  for  smooth  and  finned  tube  bundles. 

For  porous  coated  surfaces,  Czikk  et  al.  [Ref.  19]  found  no  'bundle 
effect'  over  a  wide  range  of  heat  flux  (1-100  kW/m2)  and  the  bundle  data 
agreed  closely  with  single  tube  data.  Arai  et  al.  [Ref.  20]  found  that 
the  'bundle  effect'  for  a  Thermoexcel  tube  bundle  was  smaller  than  that 
found  for  a  smooth  or  finned  tube  bundle.  However  as  before,  any  'bundle 
effect'  was  eliminated  at  high  heat  fluxes  where  the  data  for  all  the 
tubes  agreed  closely  with  single  Thermoexcel  tube  results.  These  effects 
are  similar  to  those  found  by  Czikk  et  al.  [Ref.  19]  for  the  porous 
coated. 

Chilman  [Ref.  3]  reported  experiments  with  a  Turbo-B  tube  bundle  in 
pure  R-113,  conducting  both  increasing  and  decreasing  heat  flux  tests.  He 
concluded  that  in  the  natural  convection  region,  heated  lower  tubes  do  not 
have  any  influence  on  the  heat  transfer  from  upper  tubes.  Also,  Chilman 
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reported  the  presence  of  heated  lower  tubes  within  a  bundle  reduced  the 
incipient  boiling  point. 

Stephan  and  Mitrovic  [Ref.  18]  reported  the  influence  of  oil  on  the 
boiling  heat-transfer  coefficient  of  R-12  using  a  T-shaped  finned  tube 
(Gewa-T)  bundle.  They  reported  that  the  ratio  of  oil  to  no  oil  heat- 
transfer  coefficients  decreased  with  the  mass  fraction  of  oil  for  all 
except  the  highest  heat  flux  (22  kW/m2)  where  an  increase  in  heat  transfer 
was  noted  for  oil  concentrations  between  1  and  6%  .  They  concluded  that 
the  influence  of  oil  on  heat  transfer  was  mainly  due  to  the  thermal 
properties  of  the  specific  oil  used  in  the  experiments  and  its  interaction 
with  the  refrigerant. 

Schlager  et  al.  [Ref.  21]  summarized  the  influence  of  oil  on 
refrigerant  in  pool  boiling.  They  stated  that  under  certain  conditions 
(typically  low  pressure  and  high  heat  flux),  the  heat  transfer  coefficient 
increased  at  low  oil  concentration.  Stephan  [Ref.  22]  first  pointed  this 
out  and  attributed  the  phenomenon  to  foaming.  Burkhardt  and  Hahne  [Ref. 
23]  for  a  finned  tube  bundle  found  that  the  maximum  heat  transfer 
coefficient,  which  was  10%  to  15%  above  the  oil-free  value,  occurred  at  a 
concentration  of  about  4% . 

Heimbach  [Ref.  24]  conducted  experiments  with  R-12/oil  mixtures  on  a 
finned  tube  bundle.  He  reported  that  the  presence  of  up  to  2%  oil,  did 
not  affect  the  heat  transfer  performance  significantly.  However  at  higher 
concentrations  (3%  to  7%),  an  increase  in  the  heat  transfer  was  observed. 
He  also  attributed  this  to  foaming  and  postulated  that  changes  in  the 
properties  of  the  mixture  might  facilitate  the  formation  of  bubbles. 
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Anderson  [Ref.  13]  and  Akcasayar  [Ref.  25]  conducted  experiments  with 
finned  (19  fpi)  and  High  Flux  (porous  coated)  tube  bundles  in  pure  R-114 
and  R-114/oil  mixtures  in  the  same  apparatus.  Anderson  reported  a  maximum 
heat  transfer  performance  at  an  oil  concentration  of  3%  for  the  finned 
tube  bundle.  Akcasayar  also  reported  that  the  finned  tube  bundle 
performance  increased  1.65  times  with  3%  oil  concentration  (compared  with 
pure  R-114)  at  the  maximum  heat  flux  level.  For  6%  and  10%  oil 
concentrations,  the  performance  of  the  bundle,  when  compared  with  lower 
oil  concentrations,  decreased.  This  was  especially  significant  at  a  10% 
oil  concentration.  Compared  with  the  finned  tube  bundle,  the  High  Flux 
tube  bundle  had  a  1.5  times  better  heat  transfer  performance  at  a  heat 
flux  of  30  kW/m2  in  pure  R-114.  However,  these  performance  ratios 
decreased  with  increased  oil  concentrations  such  that  the  finned  bundle 
outperformed  the  High  Flux  bundle  at  6%  oil  concentration.  This  was 
especially  true  at  the  highest  heat  fluxes  where  the  High  Flux  bundle 
performance  was  not  much  better  than  a  smooth  tube  bundle  of  similar  size. 
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III.  EXPERIMENTAL  APPARATUS 


A.  TEST  APPARATUS  OVERVIEW 

The  experimental  apparatus  including  the  auxiliary  equipment  and  the 
evaporator/condenser  is  shown  in  Figure  1.  The  following  is  only  a 
general  description  of  the  whole  experimental  apparatus.  A  more  detailed 
look  at  the  condenser  and  evaporator  is  provided  in  section  C.  Further 
information  about  the  apparatus  is  provided  by  Murphy  [Ref.  26]  and 
Anderson  [Ref.  13]. 

The  experimental  apparatus  is  essentially  made  up  of  three  closed 
loops.  The  first  loop  consists  of  an  8  ton  refrigeration  unit  located 
outside  the  laboratory  which  is  used  to  cool  an  ethylene  glycol/water 
mixture.  The  second  loop  is  the  ethylene  glycol/water  mixture  flowing 
through  the  condenser.  This  mixture  is  contained  within  a  large  sump 
withir  the  laboratory.  The  flow  rate  through  the  condenser  is  delivered 
by  two  pumps  which  can  be  operated  independently  or  together;  this  coolant 
mixture  condenses  the  refrigerant  vapor  in  the  condenser  and  maintains 
system  pressure  and  temperature .  Pump  number  one  provides  coolant  flow 
through  the  four  test  condenser  tubes  as  well  as  to  one  of  the  auxiliary 
condenser  coils  (bottom  coil) .  Pump  number  two  provides  coolant  through 
a  manifold  which  distributes  the  coolant  to  the  remaining  four  auxiliary 
condenser  coils  within  the  condenser.  The  third  loop  is  the  evaporator 
and  condenser  itself  designed  for  reflux  operation.  The  vapor  generated 
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in  the  evaporator  flows  upward  and  condenses  in  the  condenser;  the 
condensate  then  returns  to  the  evaporator  via  gravity. 

B .  AUXILIARY  EQUIPMENT 

1.  28  kW  Refrigeration  Unit 

This  unit  was  used  to  cool  a  1.8  m3  reservoir  sump  of  an 
ethylene  glycol/water  mixture  (coolant)  to  the  desired  temperature  needed 
to  condense  the  refrigerant  vapor.  For  R-114,  the  temperature  control  was 
set  to  maintain  the  sump  at  -15  °C.  The  refrigeration  unit  had  a  cooling 
capacity  of  28  kW  (8  tons). 

2  .  Ethylene  Glycol/Water  Mixture 

The  coolant  used  was  a  54%  (by  weight)  ethylene  glycol/water 
mixture.  This  refrigerated  mixture  was  used  to  control  the  system 
pressure  and  temperature  by  circulating  coolant  through  the  auxiliary 
condenser  coils  and/or  condenser  test  tubes  at  different  flow  rates. 

3 .  Pumps 

Two  pumps  were  available  to  circulate  the  coolant  from  the  sump 
through  the  condenser.  Pump  number  one  fed  four  test  condenser  tubes  and 
one  of  the  auxiliary  condenser  coils.  Pump  number  two  fed  the  other  four 
auxiliary  condenser  coils.  Pump  number  one  was  the  primary  pump  used  at 
low  heat  fluxes;  pump  number  two  was  started  (as  necessary)  at  high  heat 
fluxes  to  maintain  the  desired  saturation  pressure  in  the  evaporator. 

4 .  Flowmeters 

Five  calibrated  float-type  flowmeters,  connected  to  pump  number 
one,  were  used  to  measure  the  flow  rates  passing  through  the  four  test 
condenser  tubes  and  one  auxiliary  condenser  coil.  One  additional 
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flowmeter  (connected  to  pump  number  two)  was  used  to  measure  the  total 
flow  to  the  remaining  four  auxiliary  condenser  coils.  Each  of  the  five 
auxiliary  condenser  coils  had  a  globe  valve  to  regulate  (or  shut  off) 
flow  as  desired.  For  the  four  test  condenser  tubes,  the  coolant  flow  was 
regulated  by  a  flowmeter  valve. 

C .  EVAPORATOR/CONDENSER 

An  overall  view  of  the  evaporator  and  condenser  is  shown  in  Figure  2. 
The  evaporator  was  designed  to  simulate  a  small  portion  of  a  refrigerant 
flooded  evaporator.  Front  and  side  views  of  the  evaporator  are  shown  in 
Figures  3  and  4.  It  was  fabricated  from  stainless  steel  plate  and  formed 
into  a  short  cylinder,  610  mm  in  diameter  and  241  mm  long.  Electrically- 
heated  tubes  were  cantilever-mounted  from  the  back  wall  of  the  evaporator 
to  permit  viewing  along  the  axis  of  the  tubes  through  the  lower  of  two 
viewing  windows.  A  plexiglas  plate  was  attached  to  the  front  of  the  tube 
bundle  to  ensure  tube  alignment  during  experiments.  Each  viewing  window 
had  a  layer  of  glass  and  plexiglas,  the  glass  being  used  on  the 
refrigerant  side  in  order  to  prevent  surface  cracking  of  the  plexiglas. 
The  plexiglas  gave  the  glass  added  strength  and  served  as  a  safety  barrier 
in  case  the  glass  broke  during  pressurization. 

The  electric  power  can  be  applied  separately  to  each  set  of  heaters 
using  a  STACO  240  V,  23.5  kVA  rheostat  controller  shown  in  Figure  5. 
Also,  the  desired  number  of  instrumented  tubes,  active  tubes,  simulation 
or  auxiliary  heaters  can  be  individually  activated  by  using  circuit 
breakers.  The  five  simulation  heaters,  each  with  a  maximum  rating  of  4 
kW,  were  mounted  below  the  test  bundle  in  order  to  simulate  15  additional 
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tube  rows  in  a  larger  bundle  and  to  provide  an  inlet  vapor  quality  into 
the  bottom  of  the  test  bundle  as  suggested  by  Webb  [Ref.  27] .  The  liquid 
pool  was  maintained  at  2.2  °C  (corresponding  to  a  saturation  pressure  of 
1  atmosphere)  by  passing  coolant  through  the  condenser. 

Figure  6  is  a  schematic  sectional  view  of  the  evaporator  that  shows 
the  four  kinds  of  heated  tubes  installed  in  the  evaporator.  These  were 
instrumented,  active,  auxiliary,  and  simulation.  For  this  study  only  the 
instrumented,  active,  and  simulation  heaters  were  used;  the  auxiliary 
heaters  are  needed  for  experiments  either  at  higher  pressures  or  for  other 
refrigerants  which  have  a  higher  normal  boiling  point  (such  as  R-113) . 
Table  1  gives  the  power  rating  for  these  heaters  and  the  number  used  in 
the  evaporator. 

The  tube  bundle  itself  consists  of  instrumented,  active,  and  dummy 
tubes.  The  location  of  each  tube  is  represented  by  the  respective  letter 
I,  A,  and  D  as  shown  in  Figure  6.  The  test  bundle  consists  of  two  types  of 
heated  tubes:  12  active  tubes  (marked  "A")  which  contained  1  kW  cartridge 
heaters,  and  5  instrumented  tubes  (marked  "I")  which,  in  addition  to  the 
1  kW  cartridge  heaters,  contained  six  wall  thermocouples  each. 

In  measuring  boiling  heat  transfer  coefficients,  great  care  must  be 
exercised  with  the  cartridge  heater  and  temperature  measuring 
instrumentation  to  ensure  good  accuracy.  The  instrumented  test  tubes  were 
fabricated  using  the  same  method  as  that  used  by  Hahne  and  Muller  [Ref. 
17]  and  Wanniarachchi  et  al.  [Ref.  14).  The  exact  procedure  can  be  found 
in  Eraydin  [Ref.  28].  Figure  7  is  a  cross-sectional  sketch  of  an 
instrumented  tube,  showing  the  construction  details  and  the  location  of 
the  wall  thermocouples.  The  thermocouples  were  embedded  in  the  wall  at 
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different  circumferential  and  longitudinal  positions  along  the  heated 
section  of  the  tube  shown  in  Figure  7 . 

The  five  instrumented  tubes  were  located  along  the  centerline  of  the 
tube  bundle,  forming  a  vertical  in-line  column.  All  the  instrumented  and 
active  tubes  were  Turbo-B  tubes  made  by  Wolverine  Tube  Co.  (see  Section 
E) .  These  tubes  were  nominally  15.9  mm  in  outside  diameter  and  were 
arranged  in  an  equilateral  triangular  pitch  (ie.  centerline-to— centerline 
spacing)  of  19.1  mm,  giving  a  pitch-to-diameter  ratio  of  1.35.  The 
thickness  of  the  Turbo-B  enhancement  was  0.85  mm  giving  a  diameter  to  the 
base  of  the  enhancement  of  14.2  mm. 

The  bundle  also  contained  a  number  of  unheated  dummy  smooth  tubes 
(marked  "D")  that  were  used  to  guide  the  two  phase  mixture  through  the 
bundle.  The  dummy  tubes  were  made  from  commercially  available  15.9  mm  OD 
smooth  copper  tubing.  Two  vertical  baffle  plates  made  of  aluminum  were 
used  on  either  side  of  the  bundle  to  restrict  circulation  into  and  out  of 
the  bundle  at  the  sides.  A  dummy  rack  (Figure  8)  consisting  of  12  solid 
rods  made  of  aluminum  (15.9  mm  OD  and  spaced  19.1  mm  from  centerline-to- 
centerline)  was  placed  below  the  tube  bundle.  This  rack  had  a  triangular 
pitch  arrangement  with  vapor  retaining  plates  on  the  sides  and  was 
designed  for  two  purposes:  to  collect  all  rising  two  phase  flow  generated 
by  the  simulation  heaters  and  direct  it  into  the  test  bundle  and  to 
simulate  vapor  passing  through  a  large  bundle  before  reaching  the 
instrumented  tubes.  A  small  open  space  (approximately  5  mm  in  height)  was 
left  between  the  bundle  and  dummy  tube  rack.  This  space  allowed  some 
refrigerant  from  below  to  enter  the  bundle  and  replace  the  vapor  being 
generated  in  the  bundle.  However,  there  was  also  a  space  below  the  dummy 
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rack  that  allowed  the  majority  of  the  circulation  to  occur.  Thus, 
liquid/vapor  circulation  was  vertically  upward  over  the  five  instrumented 
test  tubes  with  no  net  horizontal  component.  Most  of  the  liquid-vapor 
mixture  after  passing  through  the  bundle  was  separated  when  it  reached  the 
pool  surface.  However,  due  to  the  strong  circulation  patterns  set  up 
within  the  liquid  pool,  some  vapor  bubbles  remained  trapped  in  the  liquid 
and  circulated  around  the  pool. 

The  condenser  included  four  test  tubes  (each  1.219  m  in  length  and 
15.9  mm  OD)  in  a  vertical  in-line  column  and  five  auxiliary  copper  coils. 
For  the  boiling  experiments,  these  tubes  were  used  to  regulate  the 
pressure  and  temperature  in  the  evaporator.  The  condenser  was  designed  to 
permit  independent  condensation  studies  of  small  in-line  tube  bundles, 
using  the  evaporator  as  a  source  of  vapor.  Details  of  the  condenser  can 
be  found  in  Mazzone  [Ref.  29]. 

D.  DATA  ACQUISITION  SYSTEM/INSTRUMENTATION 

As  described  by  Akcasayar  [Ref.  25]  ,  a  Hewlett  Packard  HP-3497A  Data 
Acquisition  System,  HP-9125  computer  and  HP-701  printer  were  used  for  data 
acquisition,  data  reduction  and  data  printing  respectively.  Although  an 
HP-9826  computer  and  HP-7470A  plotter  can  be  used  for  final  graph 
printing,  a  Macintosh  computer  (using  Criketgraph)  was  utilized.  As 
described  by  Anderson  [Ref.  13],  type-T  copper-constantan  thermocouples 
were  used  to  make  temperature  measurements  on  the  HP-3497A  using  a  relay 
multiplexer  assembly  equipped  with  thermocouple  compensation.  A  20- 
channel  relay  multiplexer  card  was  used  to  measure  the  voltage  output  from 
both  the  voltage  and  current  sensors.  The  voltage  measurements  were  taken 
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from  separate  sensors  that  measured  the  voltage  going  to  the  tube  bundle, 
the  simulation  heaters  and  the  auxiliary  heaters.  The  total  current 
going  through  the  auxiliary  and  simulation  heaters  was  measured  using  an 
American  Aerospace  Control  (ACC)  current  sensor.  The  current  to  each 
instrumented  tube  heater  was  measured  using  five  identical  current 
sensors.  The  voltage  supplied  to  the  other  active  tubes  was  also 
measured,  but  the  current  of  each  active  tube  was  not.  Instead,  the  total 
current  for  a  pair  of  active  tubes  was  measured,  and  this  was  sufficient 
since  these  tubes  each  had  the  same  power  output  (1000  W)  as  the 
instrumented  test  tube  heaters  and  there  was  no  apparent  reason  to  monitor 
each  active  tube  individually.  Computer  channel  assignments  for  data 
acquisition  and  array  assignments  are  given  in  Table  2. 

E.  GEOMETRY  OF  TURBO -B  TUBE 

The  Turbo-B  tube ,  manufactured  by  Wolverine  Tube  Inc . ,  contained  an 
enhanced  surface  geometry.  The  exterior  boiling  enhancement  is  made  by 
raising  low  integral  fins,  cutting  diagonally  across  these  fins,  and  then 
rolling  the  fins  to  compress  them  to  form  mushroom-like  pedestals  [Ref. 
30].  This  process  forms  numerous  re-entrant  passageways.  Figure  9  shows 
the  surface  of  the  tube  at  25  times  its  actual  size.  The  tube  is 
currently  available  in  copper,  cupro-nickel ,  and  low  carbon  steel. 

The  relative  dimensions  of  the  tube  used  in  this  study  are  as  follows: 
Tube  material  -  Copper 
Nominal  Outside  diameter  -  15.9  mm 
Enhanced  surface  length  -  203 . 2  mm 
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Thickness  of  Enhancement  -  0.85  mm 
Diameter  to  Base  of  Enhancement  -  14.2  mm 
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THERMOCOUPLE  **  1  ^  THERMOCOUPLE  **  2 
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Figure  4.  Side  View  of  Evaporator 
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DUMMY  TUBE  RACK 


Figure  7.  Thermocouple  Locations  on  an  Instrumented  Bolling  Tube  and 
Tube  Section  View 
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Figure  9.  Close-up  View  of  Turbo-B  Tube  Surface  (25  X) 
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Table  1.  EVAPORATOR  HEATERS 


Heater  Type 

Number 

Power  Rating  per  Heater 

Instrumented  Tube  Heaters 

5 

1000W 

Active  Tube  Heaters 

12 

1000W 

Auxiliary  Heaters 

4 

4000W 

Simulation  Heaters 

5 

4000W 

Table  2.  COMPUTER/DATA  ACQUISITION 
ASSIGNMENT 


Amperage  Sensor  Description 

Channel 

Array 

Tube  1 

30 

Amp(O) 

Tube  2 

31 

Amp(1) 

Tube  3 

32 

Amp(2) 

Tube  4 

33 

Amp(3) 

Tube  5 

34 

Amp(4) 

Active  Heater  Group  1 

35 

Amp(5) 

Active  Heater  Group  2 

36 

Amp(6) 

Active  Heater  Group  3 

37 

Amp(7) 

Active  Heater  Group  4 

38 

Amp(8) 

Active  Heater  Group  5 

39 

Amp(9) 

Auxiliary  Heaters 

25 

Amp(10) 

Simulation  Heaters 

26 

Amp(1 1) 

Voltage  Sensor  Description 

Channel 

Array 

Instrumented/Active 

27 

Volt(O) 

Simulation  Heaters 

28 

Volt(1) 

Auxiliary  Heaters 

29 

Volt(2) 

Table  2.  COMPUTER/DATA  ACQUISITION  ASSIGNMENT  (CONT. ) (cont . ) 


Thermocouple  Description 

Channel 

Vapor  1-Top  of  Condenser 

00 

Vapor  2-Top  of  Condenser 

01 

Vapor  3-Top  of  Evaporator 

02 

Liquid  1-Top  of  bundle 

03 

Liquid  2-Top  of  bundle 

04 

Liquid  3-Bottom  of  bundle 

05 

Tube  1,No.  1 

40 

Tube  1,No.  2 

41 

Tube  1,No.  3 

42 

Tube  1,No.  4 

43 

Tube  1,No.  5 

44 

Tube  1,No.  6 

45 

Tube  2, No.  1 

46 

Tube  2,No.  2 

47 

Tube  2,No.  3 

48 

Tube  2, No.  4 

49 

Tube  2,No.  5 

50 

Tube  2,No.  6 

51 

Tube  3,No.  1 

52 

Tube  3,No.  2 

53 

Tube  3,No.  3 

54 

Tube  3,No.  4 

55 

Tube  3,No.  5 

56 

Tube  3,No.  6 

57 

Tube  4,No.  1 

58 

Tube  4,No.  2 

59 

Tube  4,No.  3 

60 

Tube  4, No.  4 

61 

Tube  4, No.  5 

62 

Tube  4,No.  6 

63 

Tube  5, No.  1 

64 

Tube  5, No.  2 

65 

Tube  S.No.  3 

66 

Tube  5,No.  4 

67 

Tube  5, No.  5 

68 

Tube  5, No.  6 

69 
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IV.  EXPERIMENTAL  PROCEDURES 


A.  REMOVAL  OF  THE  TUBE  BUNDLE  AND  BUNDLE  DISASSEMBLY 

Before  starting  the  removal  of  the  tube  bundle  from  the  evaporator, 
the  front  glass  viewing  windows  were  carefully  removed.  Next,  all 
thermocouple  wires  and  tube  heater  electrical  connections  were 
disconnected.  After  this  was  completed,  the  nuts  securing  the  backing 
plate  and  support  block  were  removed,  and  the  tube  bundle  was  taken  out 
from  the  back  of  the  evaporator. 

When  the  bundle  needed  to  be  disassembled,  it  was  ensured  that  there 
was  a  clean  working  surface.  The  first  task  was  to  remove  the  plexiglas 
plate  attached  to  one  end  of  the  aluminum  baffle  plates  (ie.  the  front  of 
the  bundle  assembly)  by  four  screws.  The  ten  screws  on  the  side  of  each 
aluminum  baffle  plate  were  then  removed  (these  were  attached  to  the  dummy 
tubes  down  each  side  of  the  bundle) .  The  aluminum  plates  were  then  pulled 
off  the  bundle.  The  four  corner  dummy  tubes  (two  top  and  two  bottom) 
remained  attached  to  the  tube  bundle  support  block  as  they  were 
countersunk  into  the  block.  The  six  outer  smooth  tubes  (three  per  side) 
could  be  easily  pulled  from  the  aluminum  plates  as  they  were  attached  only 
by  the  screws  already  removed.  The  other  ten  smooth  tubes  were  then 
unscrewed  from  the  tube  bundle  support  block  as  seen  in  Figure  10.  The 
smooth  tubes  were  engraved  to  ensure  proper  identification  during 
reassembly.  With  these  tubes  and  aluminum  baffle  plates  removed,  only  the 
instrumented  and  active  enhanced  heater  tubes  remained.  These  tubes  were 
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removed  from  the  support  block  by  loosening  the  outer  O-ring  compression 
plate,  disconnecting  the  active  heater  tube  wired  pairs,  and  pulling  the 
tubes  from  the  block.  Reassembly  of  the  tube  bundle  is  done  by  reversal 
of  this  procedure. 

B.  SYSTEM  CLEAN-UP 

If  the  system  had  been  previously  operated  with  refrigeration  oil  (or 
contaminated  from  some  other  source)  it  had  to  be  thoroughly  cleaned.  To 
accomplish  this,  the  entire  apparatus  had  to  be  taken  apart  and  cleaned  in 
the  following  manner. 

After  removal  of  the  refrigerant  (R-113  by  directly  draining  into  5 
gallon  drums  via  drain  valve  R-5  (see  Figure  2  in  Chapter  3)  at  the  bottom 
of  the  evaporator  or  R-114  by  boiling  off  into  the  storage  tank  by  opening 
R— 1  and  R-7)  and  with  the  system  at  atmospheric  pressure,  all  electrical 
connections  to  the  bundle  were  disconnected  and  the  front  viewing  glass 
windows  were  removed.  An  electric  fan  was  used  for  safety  to  ensure 
proper  ventilation.  The  tube  bundle  was  then  removed  and  disassembled  as 
described  in  section  A;  the  dummy  tube  rack  was  also  removed. 

Having  removed  the  tubes  from  the  tube  bundle,  they  were  individually 
washed  with  warm  water,  rinsed  and  then  wiped  down  with  acetone.  The 
smooth  tubes  were  cleaned  with  Copper  Brite  (a  commercial  copper  cleaning 
product)  to  remove  any  oxidation.  They  were  also  wiped  down  with  warm 
water  and  then  with  acetone.  The  same  procedure  was  followed  for  the 
Turbo-B  tubes  except  they  were  not  cleaned  with  Copper  Brite  for  fear  of 
clogging  the  channels .  During  the  cleaning  process ,  a  soft  bristled 
toothbrush  was  used  to  ensure  the  enhanced  surface  was  cleaned  properly, 
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exercising  care  not  to  interfere  with  the  tube  surface.  The  evaporator 
shell  was  cleaned  in  a  similar  manner,  using  warm  water  and  acetone. 

C.  INSTALLATION  OF  THE  TUBE  BUNDLE 

Once  the  tube  bundle  had  been  cleaned  and  reassembled  (see  section  A) , 
and  before  tightening  the  backing  plate  nuts ,  the  whole  assembly  was 
carefully  guided  back  into  the  evaporator  section,  ensuring  the  plexiglas 
viewing  cover  of  the  tube  bundle  was  not  damaged.  After  the  bundle  was  in 
position,  it  was  ensured  that  the  dummy  tube  rack  was  properly  positioned 
below  the  bundle  and  that  the  vapor  thermocouple  positions  were  still  1.75 
cm  above  the  bundle.  Then,  all  the  nuts  were  tightened  equally  on  opposite 
sides  to  give  equal  compression  on  the  gasket.  To  replace  the  front 
window,  very  small,  equal  torques  (using  a  torque  wrench)  were  applied 
circumferentially  to  each  nut  on  the  outer  ring  support  in  turn.  After 
the  window  was  in  place,  each  tube  (which  extended  through  the  outer  0- 
ring  compression  plate)  was  lightly  tapped  forward  so  as  to  touch  the 
front-viewing  window.  The  backing  plate  was  then  tightened  and  the 
individual  tube  0-rings  compressed,  providing  a  good  seal  for  the  system. 
The  compression  plate  had  grooves  for  the  tube  0-rings  to  sit  in  to  help 
with  proper  alignment  and  ensure  a  good  seal . 

D.  SYSTEM  LEAKAGE  TEST 

After  the  system  was  isolated  from  the  atmosphere  and  system  integrity 
was  restored,  a  Seargent  Welch  10  SCFM  vacuum  pump  was  connected  to  the 
apparatus  (via  valves  R-l  and  R-8)  and  the  pressure  taken  down  to  25  inHg 
vacuum.  Valves  R-l  and  R-8  were  then  secured  and  the  system  was  left 
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untouched  for  at  least  24  hours  to  see  if  there  was  any  air  leakage  in. 
If  there  was  significant  leakage  (>1  inHg  over  24  hours),  then  the  vacuum 
was  broken  by  cracking  open  valve  R-2  slowly  (this  ensured  that  no 
moisture  entered  the  system) .  The  system  was  pressurized  (with  air)  to  15 
psig  through  valve  R-2.  Large  leaks  could  then  be  detected  by  simply 
listening  to  the  air  issuing  from  the  system;  small  leaks  were  detected  by 
spraying  a  soapy  water  solution  to  all  surfaces  where  leaks  were  most 
likely  to  occur  (front  viewing  glass  gaskets,  backing  plate  gasket,  all 
fittings/valves  coming  off  the  condenser/evaporator,  O-ring  seals  of  the 
bundle  tubes  etc) .  Extreme  care  must  be  taken  to  ensure  no  moisture 
enters  the  inside  of  the  heated  tubes  where  the  heater  wires  protrude. 
After  all  leaks  were  detected  and  corrected,  the  system  was  again 
subjected  to  a  vacuum  for  another  24  hour  period.  If  the  vacuum  held, 
then  the  system  was  ready  to  receive  refrigerant.  If  not,  the  above  leak 
correction  test  was  repeated. 

E.  REFRIGERANT 
1.  Fill 

a.  From  System  Storage  Tank 

A  refrigerant  storage  tank  was  used  to  store  R-114  during 
modification/repairs  to  the  system.  The  storage  tank  prevented  discharge 
of  the  R-114  into  the  atmosphere  and  made  the  experimentation  less  costly. 
To  fill  the  evaporator  with  R-114  from  the  storage  tank,  the  ethylene 
glycol/water  coolant  temperature  was  first  reduced  to  -15  °C.  The  system 
pressure  was  then  maintained  below  the  storage  tank  pressure  (vapor 
pressure  of  R-114  at  20  °C  is  approximately  15  psig)  by  circulation  of  the 
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coolant  through  the  condenser  test  tubes  and  auxiliary  coils.  Valves  R-6 
and  R-4  were  then  opened  to  uraw  the  R-114  from  the  storage  tank  to  the 
evaporator.  The  amount  of  refrigerant  that  was  transferred  was  controlled 
by  throttling  valve  R-6  to  obtain  the  desired  level  in  the  evaporator.  If 
required,  additional  R-114  could  be  transferred  from  a  68  kg  storage 
cylinder  to  the  system  using  valve  R-2  (see  section  l.b). 
b.  From  Refrigerant  Storage  Cylinder 

To  fill  the  apparatus  from  the  68  kg  storage  cylinder,  the 
system  pressure  was  reduced  in  the  same  way  as  above.  A  hose  assembly 
containing  a  Drierite  gas  purifier  was  connected  between  the  storage 
cylinder  and  valve  R-2.  A  gas  purifier  was  used  not  only  to  remove  all 
impurities,  but  also  to  remove  any  water  from  the  refrigerant.  Once  in 
place,  both  the  storage  cylinder  valve  and  R-2  were  opened  until  the 
desired  refrigerant  level  was  reached  in  the  evaporator. 

2.  Removal  to  the  Storage  Tank 

For  tube  replacement,  system  maintenance  or  system  clean  up,  the 
R-114  was  transferred  to  the  storage  tank.  The  ethylene  glycol/water 
coolant  temperature  flowing  through  the  storage  tank  was  cooled  to  -15  °C; 
valves  R-7  and  R-8  were  opened  and  the  vacuum  pump  was  turned  on  to  put 
the  storage  tank  under  vacuum.  Once  the  storage  tank  was  under  a  20  inHg 
vacuum,  valve  R-8  was  shut  and  valve  R-l  was  opened.  Next  R-114  was  slowly 
boiled  off  to  the  storage  tank  by  using  the  tube  bundle,  simulation  and 
auxiliary  heaters  at  a  heat  flux  of  600  kW/m2  (slow  boiling  is  important 
to  ensure  minimum  transfer  of  oil).  As  the  refrigerant  level  decreased, 
individual  heaters  were  turned  off  to  ensure  none  were  uncovered.  Once 
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below  the  level  of  the  heaters,  the  final  few  cm  of  R-114  was  boiled  off 
using  heat  from  the  atmosphere.  Once  all  of  the  R-114  was  transferred, 
valves  R-l  and  R-7  were  shut. 

F.  OPERATION 

1.  System  Startup,  Securing  and  Emergency  Procedures 

See  Appendix  D 

2 .  Normal  Operation 

The  evaporator  was  filled  with  R-114  to  a  level  of  approximately 
10  cm  above  the  top  tubes  in  the  bundle.  Prior  to  operating  the  system, 
the  8  ton  refrigeration  unit  was  run  for  approximately  an  hour  to  reduce 
the  ethylene  glycol/water  coolant  in  the  sump  to  a  temperature  of  -15  °C. 
The  pressure  in  the  evaporator/condenser  was  usually  12  to  15  psig  if  the 
system  had  been  secured  overnight.  As  the  sump  was  brought  to 
temperature,  the  data  acquisition  system  and  computer  were  turned  on. 
This  allowed  the  temperature  in  the  system  to  be  monitored  during  cool¬ 
down  to  saturation  conditions.  With  this  and  pump  number  one  running, 
one  auxiliary  condenser  coil  and  the  four  condenser  test  tubes  were  used 
to  bring  the  pool  down  to  a  subcooled  condition  (for  R-114,  approximately 
1  °C  on  all  three  pool  thermocouples).  Subcooling  of  the  refrigerant  was 
done  to  ensure  the  pool  had  an  evenly  distributed  temperature  prior  to 
starting  a  run.  After  reaching  this  subcooled  condition,  all  coolant 
supply  to  the  condenser  was  secured.  The  pool  was  then  allowed  to  'heat 
up'  by  conduction  from  the  surroundings.  Once  a  saturation  temperature  of 
2.2  °C  was  reached,  the  instrumented  tube(s)  (and  simulation  heaters  for 
test  7)  was/were  switched  on  and  set  to  the  desired  heat  flux  value. 
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This  lengthy  procedure  was  done  to  prevent  the  tubes  from  prematurely 
nucleating.  The  heat  flux  of  the  instrumented  tubes  was  then  slowly 
increased  at  desired  intervals  by  adjusting  the  rheostat.  For  increasing 
heat  flux,  the  data  was  taken  with  very  small  heat  flux  increments  (every 
1000  kW/m2)  ,  waiting  at  least  5  minutes  to  attain  steady  state  conditions 
at  each  heat  flux.  At  all  regions  of  the  boiling  curve  (and  especially 
near  the  onset  of  nucleate  boiling) ,  two  readings  were  taken  at  each  heat 
flux  to  ensure  accuracy.  The  bundle  was  continuously  visually  monitored 
through  the  observation  windows.  Figure  11  shows  the  tube  bundle 
arrangements  used  during  the  experimentation.  Test  one  was  with  only  one 
instrumented  tube  turned  on  at  any  position  within  the  bundle.  Test  two 
was  with  instrumented  tubes  one  and  two  active  in  the  bundle.  Test  three 
was  with  instrumented  tubes  one,  two  and  three  active  in  the  bundle.  Test 
four  was  with  instrumented  tubes  one,  two,  three,  and  four  active  in  the 
bundle.  Test  five  was  with  all  five  instrumented  tubes  active  in  the 
bundle.  Test  six  was  all  five  instrumented  tubes  plus  all  five  pairs  of 
active  enhanced  heater  tubes  active  in  the  bundle.  Test  seven  was  the 
same  as  test  six  with  the  addition  of  all  five  simulation  heaters  active. 
For  each  data  set,  the  five  simulation  heaters  had  the  same  heat  flux  as 
the  tubes  within  the  bundle. 

G.  OIL  ADDITION 

During  the  bundle  experiments,  successive  amounts  of  York-C  oil  were 
added  into  the  evaporator.  Since  the  weight  of  the  refrigerant  in  the 
evaporator  was  60.3  kg,  the  amount  of  oil  corresponding  to  1%  by  weight 
was  measured  as  670  ml,  2%  1340  ml  etc.  The  oil  was  syphoned  into  the 
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evaporator  via  a  funnel/hose  connection  through  valve  R-3  by  reducing  the 
pressure  in  the  evaporator  to  less  than  15  inHg  vacuum.  Ensuring  that  no 
air  entered  the  system,  valve  R-3  was  promptly  shut  when  the  desired 
amount  of  oil  had  been  added. 

U.  DATA  REDUCTION  PROCEDURES 

The  data  reduction  program  "DRP4RH"  was  used  during  the  experiments 
for  processing  the  data  collected  (  see  Appendix  E  for  listing).  The 
program  was  written  in  HP  Basic  3.01  and  run  on  an  HP-9000  series 
computer.  The  characteristics  and  capabilities  of  this  software  are 
similar  to  those  provided  by  Anderson  [Ref.  13] .  The  following 
modifications  were  made: 

1.  Correction  for  pool  height  by  Chilman  [Ref.  3] 

2.  Installation  of  new  thermocouple  at  the  bottom  of  the  liquid  pool 
(bundle  inlet  temperature)  by  Chilman  [Ref.  3] 

3.  The  ability  to  obtain  data  from  one  instrumented  tube  at  any  position 
in  the  bundle . 

4.  The  natural  convection  correlation  of  Churchill  and  Chu  [Ref.  31]  for 
a  single  horizontal  cylinder  in  an  'infinite'  liquid  pool  was  added  for 
comparison  with  experimental  data. 
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Figure  11.  Tube  Bundle  Arrangements  used  During  Experimentation 
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V.  RESULTS  AND  DISCUSSION 


A.  INTRODUCTION 

The  results  are  presented  in  four  sections  with  subs-sections  as 
appropriate.  The  first  section  discusses  the  preliminary  experiments 
which  led  to  modifying  the  experimental  start-up  procedure  to  include 
subcooling.  The  second  section  discusses  the  natural  convection  effects, 
nucleate  pool  boiling  phenomena  and  hysteresis  effects  within  the  tube 
bundle  in  pure  R-114.  The  third  section  discusses  similar  phenomena  for 
R-114/oil  mixtures  and  their  effects  on  the  above.  The  fourth  section 
shows  comparisons  of  data  taken  during  this  thesis  with  previously 
obtained  data  at  the  Naval  Postgraduate  School. 

A  list  of  data  files  taken  during  this  investigation  may  be  found  in 
Appendix  A.  All  data  files  used  in  this  thesis  use  the  following  filename 
sequence.  Each  file  is  composed  of  five  sets  of  alpha-numeric  characters 
used  to  describe  the  experiment. 

First  set  (2  char.)  Tube  Type  TB  (Turbo-B) 

Second  set  (1  char.)  Heat  Flux  I  (Increasing) 

D  (Decreasing) 

Third  set  (2  char.)  Oil  Percent  00  (0%) 

01  (1%) 

02  (2%) 

03  (3%) 

06  (6%) 

10  (10%) 

Fourth  set  (2  char.)  Test  Type  01  (test  1) 

02  (test  2) 

03  (test  3) 

04  (test  4) 
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05  (test  5) 
06  (test  6) 
07  (test  7) 


Fifth  Set  (1  char.)  Additional  tests  A-Z  (If  conducted) 

To  give  an  example,  the  filename  TBI0107  means  "Turbo-B  tube  bundle, 
increasing  heat  flux,  1%  R-114/oil  mixture  and  test  number  7".  If  more 
detail  is  desired  about  a  specific  data  set,  see  Appendix  A.  All  plot 
filenames  are  similar  to  data  filenames  except  they  start  with  the  letter 
"P" .  The  test  numbers  are  shown  in  Figure  11  in  Chapter  IV. 

All  graphs  are  plotted  showing  heat  flux  (W/m2)  along  the  ordinate  (y 
axis)  as  a  function  of  wall  superheat  (K)  along  the  abscissa  (x  axis) . 
The  wall  superheat  is  defined  as  the  difference  between  the  'corrected' 
average  tube  wall  temperature  (ie.  having  accounted  for  depth  of 
thermocouple  burial)  and  the  local  liquid  saturation  temperature 
(corrected  for  hydrostatic  head  within  the  bundle) .  The  heat  flux  was 
corrected  to  account  for  the  heat  lost  through  the  unheated  tube  ends. 
The  heat  flux  was  varied  from  600  to  100,000  W/m2  for  increasing  heat  flux. 
To  ensure  greater  detail  at  the  point  of  incipience,  the  heat  flux  was 
increased  in  small  steps;  these  settings  varied  from  test  to  test.  The 
heat  flux  values  for  decreasing  experiments  were  taken  at  prescribed 
settings  for  easy  comparison  with  past  experiments  and  future  reference. 
These  heat  flux  settings  were  1x10s,  7.5x10*,  5x10*,  3x10*,  2x10*,  1.5x10*, 
1x10*,  7x103,4x103,  2x103,  1x103  W/m2.  Approximately  30-40  data  points  were 
taken  for  each  increasing  heat  flux  run  and  20-25  points  were  taken  for 
each  decreasing  heat  flux  run. 
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B.  PRELIMINARY  EXPERIMENTS 


After  the  pure  R-113  was  removed  following  Chilman's  [Ref.  3] 
experiments  and  the  apparatus  and  system  cleaned,  pure  R-114  was  added 
from  the  storage  cylinders.  Five  tests  were  conducted  using  this  R-114 
and  these  are  shown  in  Figure  12.  The  first  test  (TBI0001A)  was  test  one 
with  the  top  tube  activated.  The  procedure  described  by  Eraydin  [Ref.  28] 
was  followed,  but  the  plot  of  test  TBI0001A  and  Eraydin' s  data  for  test 
one  (also  shown  in  Figure  12)  produced  significantly  different  results  in 
the  natural  convection  (NC)  region.  The  data  of  Eraydin  show  a  greater 
heat  transfer  coefficient  (lower  wall  superheat)  than  test  TBI0001A  which 
show  results  closer  to  the  Churchill/Chu  [Ref.  31]  correlation  (C/C)  for 
natural  convection.  The  only  difference  in  the  apparatus  between  test 
TBI0001A  and  Eyradin's  experiments  was  the  addition  of  a  third 
thermocouple  at  the  bottom  of  the  pool.  For  test  TBI0001A,  all  three  pool 
temperature  readings  were  within  +/-  0.1  °C  prior  to  recording  data.  For 
the  data  of  Eraydin,  only  the  temperature  at  the  top  of  the  pool  could  be 
checked. 

Test  TBI0001B  was  thought  to  be  a  repeat  of  test  TBI0001A,  but  for  the 
bottom  tube  in  the  bundle  (tube  5  only).  However,  upon  observation  and 
investigation,  it  was  found  that  two  tubes  (tube  one  and  tube  five)  were 
activated  due  to  the  way  the  program  DRP4  was  set  up.  Hence  the  data 
presented  is  for  tube  one  (the  top  tube)  with  tube  five  (the  bottom  tube) 
activated  as  well.  The  program  was  then  modified  to  obtain  data  for  a 
single  tube  (test  one)  at  any  position  within  the  bundle. 

Test  TBI0001C  was  conducted  using  the  bottom  tube  (tube  five)  as  a 
single  tube.  Again  following  Eraydin's  [Ref.  28]  procedure,  partial 
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nucleation  was  observed  immediately ,  explaining  why  the  data  lie  well  to 
the  left  of  the  Churchill/Chu  correlation  (ie.  a  higher  heat  transfer 
coefficient)  .  It  can  be  seen  that  test  TBI0001C  is  in  good  agreement  with 
the  single  tube  of  Eraydin  (top  tube) .  The  reason  for  this  is  probably 
due  to  the  start-up  conditions,  which  were  not  carefully  monitored  for  any 
of  these  tests. 

Test  TBI0001D  was  conducted  using  only  tube  three  and  also  followed 
the  procedure  of  Eraydin.  This  displayed  similar  behavior  as  test 
TBI0001A  (tube  one)  except  that  nucleation  was  delayed,  ie.  occurred  at  a 
higher  wall  superheat. 

All  experiments  thus  far  were  conducted  following  Eraydin' s  procedure 
for  pure  R-114.  It  was  next  decided  to  vary  his  procedure  slightly.  For 
test  TBI0001E,  the  pool  was  first  subcooled  slightly  to  1°C.  This  ensured 
that  any  nucleation  sites  were  deactivated  and  that  the  whole  pool  was  at 
an  even  temperature.  The  pool  was  then  brought  up  to  the  required 
saturation  temperature.  However,  partial  nucleation  was  still  observed 
during  the  run.  The  plot  of  TBI0001E  is  similar  to  test  TBI0001C  and 
Eraydin 's  test  one. 

With  all  of  these  confusing  results,  it  was  decided  to  empty,  clean 
and  recharge  the  evaporator  with  fresh  R-114.  Upon  boiling  off  the  'old' 
R-114,  a  small  quantity  of  oil  contamination  was  found  in  the  bottom  of 
the  evaporator.  This  could  have  been  either  vacuum  pump  oil  (which  may 
have  leaked  into  the  apparatus)  or  refrigerant  oil  that  entered  the  system 
with  the  R-114  charge  (which  should  be  minimal).  A  third  possibility 
could  be  oil  from  previous  refrigerant/oil  mixture  experiments.  However, 
this  is  not  likely  since  the  system  had  been  completely  stripped  and 
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cleaned  twice  since  the  last  mixture  tests  had  been  conducted.  Samples  of 
each  type  of  oil  (vacuum  pump  oil,  miscible  refrigerant  oil  (York-C)  and 
a  sample  of  the  contamination)  were  sent  away  for  analysis  and  the  results 
were  still  forthcoming  at  the  time  of  writing  this  thesis.  From  the  color 
of  the  contamination,  it  would  appear  to  have  originated  from  the  vacuum 
pump.  This  would  seem  to  be  more  logical  since  Chilman  had  experienced 
problems  with  the  vacuum  pump. 

The  system,  including  the  bundle,  were  cleaned  thoroughly  and  fresh 
R-114  was  added.  However,  a  gas  purifier  was  utilized  to  ensure  that  only 
pure,  clean  R-114  (with  no  moisture)  was  added  (see  Chapter  IV  for  further 
details).  Test  TBI0001F  then  was  conducted  following  Eraydin's  [Ref.  28] 
procedure  with  no  subcooling  and  'ignoring'  the  bottom  pool  thermocouple 
value .  The  data  showed  good  agreement  with  Eraydin  in  the  natural 
convection  region,  but  a  lower  heat-transfer  coefficient  in  the  boiling 
region.  A  possible  explanation  for  this  lower  heat  transfer  is  the  fact 
that  each  individual  tube  was  fully  cleaned  prior  to  adding  the  new  R-114 
and  the  surface  characteristics  may  have  been  modified  in  some  way. 

Test  TBI0001G  was  conducted  using  the  same  procedure  except  the  pool 
was  initially  subcooled  down  to  1°C  for  30  minutes.  This  ensured  the  pool 
had  an  even  temperature  distribution  throughout.  The  data  were  now  much 
closer  to  the  Churchill/Chu  [Ref.  31]  correlation.  However,  some 
premature  nucleation  was  still  observed.  Test  TBI0001J  repeated  the  above 
test  with  the  pool  temperature  subcooled  to  1°C  for  at  least  one  hour  to 
further  deactivate  any  remaining  nucleation  sites  within  the  bundle  prior 
to  starting  experiments.  The  data  then  agreed  with  the  Churchill/Chu 
correlation  as  seen  in  Figure  12.  It  became  apparent  that  premature 
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nucleation  could  affect  the  natural  convection  data  significantly. 
Therefore,  for  all  subsequent  tests,  this  same  procedure  was  adopted  with 
strict  observation  of  the  bundle  to  ensure  no  premature  nucleation 
occurred. 

C.  PURE  R-114  TURBO -B  TUBE  BUNDLE  EXPERIMENTS 
1.  Test  One  for  Different  Tube  Positions 

The  first  set  of  experiments  conducted  were  performance  tests  in 
pure  R-114.  Figure  13  shows  increasing  heat  flux  for  a  single  tube  within 
the  bundle  at  different  positions  (positions  1,  3,  and  5)  while  Figure  14 
shows  corresponding  data  (including  typical  uncertainties)  for  decreasing 
heat  flux.  All  three  tube  positions  agree  closely  with  the  Churchill/Chu 
[Ref.  31]  correlation  in  the  natural  convection  region.  The  difference 
between  position  1  and  position  3  and  5  may  be  that  position  1  is  affected 
by  the  fact  that  the  flow  is  'free'  to  expand  after  leaving  the  bundle. 
This  difference  may  also  be  due  to  wall  temperature  uncertainty  due  to 
differences  in  the  fabrication  process  (see  uncertainty  analysis  Appendix 
C) ;  however  in  the  natural  convection  region,  this  uncertainty  is  low  due 
to  relatively  high  values  of  wall  superheat.  Figure  13  also  shows  that 
tube  position  within  the  bundle  may  influence  the  point  of  incipience. 
Bergles  and  Rohsenow  [Ref.  32]  have  studied  the  incipient  point  in  more 
detail.  They  concluded  that  nucleation  was  controlled  by  a  nucleation 
parameter,  N,  given  by 

(<0(T8)(vfB) 

N  -  - 

(PvXhfg) 
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which  creates  an  incipient  boiling  superheat  given  by 

(Th  -  Ts)  -  2(N)/r . 

In  the  above  expression,  r  is  the  local  bubble  radius.  Calculation  of  the 
nucleation  parameter  for  pure  R-114  showed  that  as  saturation  pressure 
increases,  the  nucleation  parameter  decreases.  Assuming  that  the  radius 
of  curvature  of  the  forming  bubbles  is  constant  (which  is  reasonable  for 
a  Turbo-B  surface  which  has  large,  regularly  spaced  cavities),  then  (Tw  - 
Ts)  also  decreases  and  nucleation  may  be  expected  to  occur  earlier.  This 
certainly  seems  to  be  verified  in  Figure  13  where  the  incipient  point 
occurs  earlier  (lower  wall  superheat)  for  a  lower  tube  ie.  where  there  is 
an  increase  in  the  local  saturation  pressure.  You  et  al.  [Ref.  33]  also 
showed  a  decrease  (approx.  30%)  in  the  average  incipient  superheat  as 
pressure  was  increased  from  1  to  1.5  bar  for  pool  boiling  of  FC-72  on  a 
single  tube,  offering  some  other  experimental  verification  for  this 
conclusion. 

Chilman  [Ref.  3]  conducted  test  one  for  R-113  using  the  top  tube 
only  and  varied  the  local  saturation  pressure  by  varying  the  pool  height 
in  the  evaporator.  He  found  that  the  point  of  incipience  was  delayed  when 
increasing  the  liquid  pool  height  (ie.  the  hydrostatic  pressure  head). 
Boundary  layer  effects  due  to  different  liquid  circulation  patterns  may 
have  caused  this  delay  in  nucleation  and  more  research  is  certainly  needed 
to  fully  understand  the  influence  of  pressure  on  the  point  of  incipience. 
One  experiment  that  could  be  conducted  would  be  to  vary  the  pool  height, 
but  keep  the  local  pressure  at  each  tube  constant  by  simultaneously 
varying  the  vapor  pressure  above  the  pool. 
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Once  nucleacion  occurs,  Figure  13  shows  that  the  single  tube 
experiments  merge  onto  a  single  boiling  curve.  In  this  region,  there 
appears  to  be  no  effect  of  hydrostatic  pressure  head.  Figure  14  shows  cbe 
corresponding  decreasing  heat  flux  data  for  a  single  tube  at  the  same 
three  positions  within  the  bundle.  It  shows  no  significant  influence  of 
tube  position  in  the  bundle  for  decreasing  heat  flux.  Note  that  at  low 
heat  flux,  the  experimental  uncertainty  in  heat  flux  and  wall  superheat  is 
the  largest  (see  Appendix  C) . 

2.  Test  Two  to  Test  Seven 

Figures  15  to  20  show  test  two  to  test  seven  for  increasing  heat 
flux  with  pure  R-114.  Also  shown  for  comparison  in  each  figure  is  the 
Churchill/Chu  [Ref.  31]  correlation  although  this  is  only  truly  valid  for 
a  single  tube  in  an  infinite  pool.  Figure  15  shows  good  agreement  with 
the  Churchill/Chu  correlation  in  the  natural  convection  (NC)  region  and 
shows  no  effect  of  the  lower  tube  on  the  upper  tube.  The  incipient  point 
occurs  approximately  at  the  same  wall  superheat  for  both  tubes ;  once 
boiling  the  lower  tube  has  the  higher  heat-transfer  coefficient.  This  is 
contrary  to  the  results  obtained  by  Chilman  [Ref.  3]  and  Anderson  [Ref. 
13]  with  pure  R-113,  where  the  higher  tube  had  the  better  heat  transfer. 
The  reason  for  this  difference  is  not  known,  but  may  be  due  to  the 
explosive  nature  of  incipience  for  R-114  compared  to  the  more  gradual 
partial  incipience  for  R-113.  For  R-114  experiments,  the  pool  was 
subcooled  by  1  °C  for  only  about  hour  while  for  the  R-113  experiments,  the 
pool  was  left  in  a  subcooled  state  since  the  previous  experiment  (for  R- 
113  experiments,  the  pool  is  'heated  up'  to  saturation  conditions).  This 
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difference  in  subcooling  may  significantly  affect  the  nature  of  the 
observed  incipience.  Further  research  should  be  conducted  at  the 
incipient  point  to  address  some  of  the  questions.  Figure  16  for  three 
tubes  activated  shows  similar  behavior  as  Figure  15  (ie.  during  nucleate 
boiling,  the  lowest  tube  has  the  highest  heat-transfer  coefficient). 
Also,  Figure  16  shows  that  the  lowest  tube  seems  to  nucleate  last. 
Figures  17  and  18  for  four  and  five  tubes  activated  show  similar  trends  in 
both  the  NC  and  boiling  regions.  It  appears  that  the  tubes  nucleate  in 
order  down  the  bundle  (ie.  the  top  tube  nucleates  at  the  lowest  wall 
superheat  and  the  bottom  tube  nucleates  at  the  highest  wall  superheat) . 
For  tests  six  and  seven  (Figures  19  and  20  respectively) ,  the  maximum 
controllable  heat  flux  was  less  than  tests  one  through  five  due  to  the  use 
of  smooth  tubes  in  the  condenser  limiting  the  condensate  rate  (and  hence 
pressure)  in  of  the  vapor  space  (if  enhanced  tubes  had  been  used,  this 
could  have  been  increased).  Figure  19  (test  six)  is  consistent  with  the 
above  trends .  Furthermore ,  the  effect  of  activating  the  whole  bundle 
seems  to  cause  the  lowest  tube  to  nucleate  at  a  lower  heat  flux  and  wall 
superheat.  It  also  appears  that  there  may  be  some  influence  of  tube 
position  in  the  boiling  region;  however,  this  is  probably  due  to 
inaccuracy  in  the  wall  temperature  measurements  (see  Appendix  C) .  When 
the  simulation  heaters  are  also  activated  (test  seven,  Figure  20)  the 
trends  in  the  NC  region  are  similar  (ie.  no  effect  of  lower  tubes  on  upper 
tubes).  Full  nucleation  of  the  bundle  however,  seems  to  occur  earlier. 

Figure  21  shows  the  data  from  tube  one  for  all  seven  tests  with 
increasing  heat  flux.  This  figure  is  of  more  fundamental  interest  as  it 
shows  the  same  tube  under  different  bundle  conditions.  It  therefore  gives 
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a  better  direct  comparison  of  the  effect  of  heated  lower  tubes  as  any 
uncertainty  in  the  tube  wall  measurements  (the  largest  error  in  the 
experimental  data)  is  the  same  for  each  test  (ie.  any  effects  seen  in  the 
data  are  bundle  effects).  In  the  NC  region,  tube  one  alone  (test  one)  is 
somewhat  different.  This  may  be  due  to  'expansion'  of  the  flow  as  it 
leaves  the  top  of  the  bundle  (ie.  where  the  velocity  of  the  flow  has 
slowed  down)  or  due  to  convective  effects  by  the  addition  of  another  tube. 
For  test  two,  Figure  21  shows  an  effect  of  the  lower  tube  on  tube  one 
performance  in  the  NC  region.  For  test  three  and  all  subsequent  tests  no 
further  improvement  is  seen.  It  appears,  therefore  that  in  the  natural 
convection  region,  an  upper  tube  is  affected  by  a  lower  tube  directly 
below;  however,  when  additional  lower  tubes  are  heated,  there  is  no 
further  increase  in  performance  of  the  top  tube.  There  is  also  no  effect 
on  the  incipient  point  (apart  from  test  one  mentioned  above) .  In  the  high 
heat  flux  boiling  region,  there  is  also  little  enhancement  due  to  the 
lower  tubes.  This  is  to  be  expected  in  an  enhanced  tube  bundle,  where  the 
total  heat  transfer  at  high  heat  fluxes  is  primarily  due  to  nucleation 
from  the  tube  surface  itself,  rather  than  from  convection  around  the 
surface  from  the  tubes  below. 

Figures  22  to  27  show  tests  two  to  seven  for  decreasing  heat  flux 
in  pure  R-114.  When  comparing  tube  one  to  tube  two,  Figure  22  shows  a 
increase  in  heat  transfer  performance  of  tube  one  by  tube  two  in  the 
boiling  region  while  Figure  IS  (increasing  heat  flux)  showed  the  opposite 
effect.  The  most  probable  reason  for  this  crossover  is  that  these  two 
experiments  (TBI0002)  and  (TBD0002)  were  conducted  on  different  days  and 
startup  procedures  were  slightly  different.  For  TBI0002  (increasing  heat 
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flux),  the  test  was  conducted  as  outlined  in  Chapter  IV  section  F  (ie. 
subcooled  to  1  °C,  gradually  heated  up  with  data  taken  over  a  period  of 
approximately  4  hrs).  For  TBD0002  (decreasing  heat  flux),  the  pool  was 
subcooled  to  1  °C,  and  then  tube  one  and  two  were  turned  on  to  the  highest 
heat  flux  (105  W/m2)  and  allowed  to  heat  up  for  30  minutes.  The  total  time 
boiling  for  TBI0002  at  the  highest  heat  flux  was  therefore  less  than  for 
TBD0002.  For  all  other  experiments,  increasing  followed  by  decreasing 
runs  were  conducted  on  the  same  day  approximately  15  minutes  apart.  More 
research  should  be  conducted  to  investigate  nucleation  site 
activation/deactivation. 

At  high  heat  fluxes  for  test  three,  Figure  23  shows  no  heat 
transfer  performance  improvement  between  the  tubes.  However,  at  low  heat 
fluxes ,  there  does  appear  to  be  an  improvement  on  tube  one  and  two  from 
tube  three.  Figures  24  to  27  show  similar  trends  at  high  and  low  heat 
fluxes  (ie.  the  top  tubes  are  further  enhanced  by  lower  tubes  at  low  heat 
fluxes).  It  should  also  be  noted  that  the  lowest  tube  in  any  specific 
decreasing  heat  flux  test  had  the  worst  performance. 

Figure  28  compares  tube  one  only  for  tests  one  to  seven  for 
decreasing  heat  flux.  As  stated  above,  there  appears  to  be  a  definite 
tube  enhancement  at  low  heat  fluxes  with  little  or  no  enhancement  at  high 
heat  fluxes.  This  is  probably  due  to  convective  effects  which  tend  to 
increase  the  heat  transfer  performance  of  the  upper  t..jes  due  to  the 
presence  of  lower  tubes  (ie.  bubbles  coming  from  the  lower  tubes  impinge 
and  slide  over  the  upper  tubes  and  increase  the  heat  transfer) .  At  high 
heat  fluxes,  on  the  other  hand,  all  the  tubes  are  nucleating  so  vigorously 
that  these  'sliding'  bubbles  have  little  or  no  noticeable  affect  on  the 
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overall  performance.  This  supports  the  hypothesis  of  Cornwell  [Ref.  7] 
that  total  heat  transfer  in  a  bundle  is  due  to  a  summation  of  convective 
and  nucleation  heat  transfer  phenomena.  Similar  trends  were  found  by 
Anderson  [Ref.  13]  and  Akcasayar  [Ref.  25]  for  smooth  and  finned  tube 
bundles  respectively  (using  the  same  apparatus)  and  also  by  Aral  et  al. 
[Ref.  20]  for  a  Thermoexcel-E  tube  bundle.  However,  Akcasayar  [Ref.]  did 
not  find  such  an  enhancement  effect  for  a  High  Flux  tube  bundle  indicating 
that  at  low  heat  fluxes,  a  porous  coated  surface  already  has  a  significant 
number  of  active  nucleation  sites  such  that  impinging  bubbles  from  below 
have  little  or  no  added  effect.  Turbo-B  is  more  similar  to  a  Thermoexcel- 
E  surface  and  at  low  heat  fluxes ,  these  two  types  of  surface  obviously 
exhibit  different  nucleation  characteristics  to  those  of  a  porous  coated 
tube . 


D.  R- 114/OIL  MIXTURES  TURBO-B  TUBE  BUNDLE  EXPERIMENTS 
1.  Tests  with  1%  and  2%  oil 

Only  four  experiments  were  conducted  with  a  1%  and  2%  R-114/oil 
mixture.  These  were  tests  one  and  seven  for  both  increasing  and 
decreasing  heat  flux;  no  experiments  were  conducted  for  tests  two  through 
six.  Figure  29  shows  test  one  at  1%  oil  concentration  for  increasing  and 
decreasing  heat  flux,  clearly  showing  a  hysteresis  pattern.  Compared  with 
pure  R-114  (Figures  13  and  14) ,  Figure  29  shows  no  apparent  effect  of  oil 
on  the  heat  transfer  in  either  the  NC  or  boiling  regions. 

Figure  30  shows  test  seven  for  increasing  heat  flux  (1%  oil) . 
Compared  with  Figure  20  (pure  R-114),  there  are  similar  trends  (ie.  no 


49 


effect  in  the  NC  or  boiling  regions).  The  tubes  again  appear  to  be 
nucleating  'in  order'  down  the  bundle,  as  found  with  test  seven  in  pure  R- 
114.  For  decreasing  heat  flux,  Figure  31  shows  similar  trends  to  test 
seven  in  pure  R-114  (Figure  27).  Thus  a  1%  oil  concentration  appears  to 
have  little  or  no  effect  on  bundle  performance  for  both  It  creasing  and 
decreasing  heat  flux. 

Figure  32  shows  test  one  for  a  2%  oil  concentration  for 
increasing  and  decreasing  heat  flux.  Figures  33  and  34  show  test  seven 
for  increasing  and  decreasing  heat  flux  respectfully  for  2%  oil 
concentration.  All  three  graphs  (Figures  32  to  34)  are  similar  to  those 
for  pure  and  1%  oil  concentrations,  showing  that  2%  oil  also  has  little 
effect  on  overall  bundle  heat  transfer  performance. 

2.  Tests  with  3%  oil 

Nine  experiments  were  conducted  with  a  3%  R-114/oil  mixture.  In 
addition  to  test  one  and  seven  (conducted  for  both  increasing  and 
decreasing  heat  flux  as  before) ,  tests  two  through  six  were  conducted  for 
decreasing  heat  flux  only.  Figure  35  shows  test  one  with  3%  oil  for  both 
increasing  and  decreasing  heat  flux.  Again,  the  figure  clearly  shows 
hysteresis  effects  between  the  increasing  and  decreasing  experiments.  As 
with  previous  oil  percentages  (Figures  13,  14,  29,  and  32),  it  shows  there 
is  no  apparent  effect  of  oil  in  the  NC  region,  similar  to  previous  test 
one  data  for  other  oil  percentages.  However,  at  the  highest  heat  flux 
(100  kW/m2)  there  is  an  increase  in  the  heat  transfer  of  about  10%.  This 
is  similar  to  the  increases  found  by  Burkhardt  and  Hahne  [Ref. 23]  in  a 
finned  tube  bundle  and  Aral  et  al.  [Ref.  20].  Figure  36  shows  test  seven 
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for  increasing  heat  flux.  The  Churchill/Chu  [Ref.  31]  correlation  for 
pure  refrigerant  is  plotted  for  comparison  only.  Agreement  is  good, 
demonstrating  that  in  the  NC  region,  in  addition  to  there  being  no  effect 
of  tube  position,  there  is  also  no  apparent  effect  of  oil  concentration  on 
the  heat-transfer  coefficient.  As  before,  the  tubes  appear  to  be 
nucleating  'in  order'  (ie.  top  tube  nucleates  first  with  the  bottom  tube 
nucleating  last) . 

Figures  37  to  42  show  data  from  tests  two  to  seven  for  decreasing 
heat  flux  only.  All  show  no  effect  of  lower  tubes  on  upper  tubes  (within 
the  bundle)  in  the  boiling  region  at  high  heat  fluxes  (the  small  amount  of 
scatter  is  probably  due  to  inaccuracies  in  the  wall  temperature 
measurements).  Each  successive  figure  shows  that  the  lowest  tube  has  the 
lowest  heat-transfer  coefficient;  this  tube  is  then  enhanced  by  the 
activation  of  tubes  below  it.  Again  it  should  be  noted  that  the 
experimental  uncertainty  is  larger  at  low  heat  fluxes.  At  all  oil 
concentrations,  tube  five  is  seen  to  have  the  lowest  heat  transfer 
performance.  According  to  Chilman  [Ref.  3],  tube  five  had  the  highest 
uncertainty  in  the  wall  temperature  measurements  and  this  might  be  the 
cause  of  this  discrepancy. 

If  one  compares  Figure  42  for  3%  oil  with  Figure  28  for  pure  R- 
114,  it  can  again  be  seen  that  there  is  a  small  increase  in  the  bundle 
heat-transfer  coefficient  for  the  R-114/oil  mixture  at  the  highest  heat 
fluxes.  For  all  tests  with  oil  added,  significant  foaming  was  observed  at 
the  pool  surface  and  this  may  be  the  cause  of  this  increase  in  heat 
transfer.  Schlager  et  al.  [ref.  21]  in  their  review  article  point  out 
that  for  certain  conditions  (typically  low  pressure  and  high  heat  flux), 
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the  heat-transfer  coefficient  increases  at  low  oil  concentrations;  they 
attributed  this  to  foaming.  Figure  43  compares  test  one  to  seven  for  tube 
one  for  decreasing  heat  flux.  As  before  with  pure  R-114  (Figure  28), 
there  appears  to  be  a  definite  increase  in  performance  of  the  upper  tubes 
by  lower  tubes  at  the  low  heat  fluxes  due  to  convection  effects,  with 
little  or  no  such  increase  at  high  heat  fluxes. 

3.  Tests  with  6%  oil 

The  same  nine  tests  as  with  3%  oil  were  conducted  with  a  6%  R- 
114/oil  mixture.  Figure  44  shows  test  one  for  both  increasing  and 
decreasing  heat  flux.  It  clearly  shows  a  hysteresis  'loop'  between 
increasing  and  decreasing  experiments.  In  comparison  with  other  oil 
concentrations,  the  point  of  incipience  occurs  at  a  slightly  lower  heat 
flux.  There  also  appears  to  be  a  small  degradation  in  performance  (10-15% 
compared  with  3%  oil  concentration)  at  the  highest  heat  flux  (100  kW/m2) 
due  to  the  oil,  but  there  is  no  apparent  effect  in  the  NC  region.  Figure 
45  shows  test  seven  for  increasing  heat  flux.  As  before,  there  is  no 
apparent  effect  of  oil  on  the  heat  transfer  in  the  NC  region  and  the  tubes 
appear  to  be  nucleating  'in  order'.  The  point  of  incipience  also  seems 
unaffected  by  the  presence  of  the  oil. 

Figures  46  to  51  show  tests  two  to  seven  for  decreasing  heat 
flux.  At  the  highest  heat  fluxes,  there  is  a  similar  small  degradation 
in  the  heat  transfer  as  found  with  test  one  (10-15%)  when  compared  with  a 
3%  oil  concentration  (Figures  37  to  42) .  When  compared  with  pure  R-114 
(Figures  22  to  27),  there  is  neither  enhancement  nor  degradation, 
indicating  that  any  enhancement  provided  by  3%  oil  is  offset  by  6%  oil. 
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At  low  heat  fluxes,  the  data  are  not  only  very  similar  to  that  for  pure 
refrigerant,  but  also  to  the  other  R-114/oll  mixtures  (ie.  at  low  heat 
fluxes,  there  is  no  effect  on  heat  transfer  at  any  oil  concentration). 
Figure  52  compares  tests  one  to  seven  for  tube  one  for  decreasing  heat 
flux.  As  before,  there  appears  to  be  the  same  convective  enhancement  at 
low  heat  fluxes  with  no  enhancement  (due  to  the  successive  activation  of 
lower  tubes  within  the  bundle)  at  high  heat  fluxes. 

4.  Tests  with  10%  oil 

The  same  nine  tests  were  repeated  for  an  R-114/oil  mixture  with 
10%  oil.  Figure  53  shows  test  one  for  both  increasing  and  decreasing 
heat  flux.  Incipience  occurred  at  a  slightly  heat  higher  flux  than  both  3% 
and  6%  oil  concentrations  indicating  that  there  appears  to  be  no 
systematic  increase  or  decrease  in  this  point  with  increase  in  oil 
concentration.  More  importantly,  there  is  a  significant  decrease  in  the 
heat  transfer  at  the  highest  heat  fluxes  (20%)  when  compared  with  pure  R~ 
114.  This  is  probably  due  to  the  re-entrant  channels  becoming  'clogged' 
with  oil  as  the  R-114/oil  mixture  is  'transported'  to  the  surface  at  a 
high  rate.  Figure  54  shows  test  seven  for  increasing  heat  flux.  As 
before,  this  shows  that  the  NC  region  is  unaffected  by  either  oil 
concentration  or  lower  tubes  in  the  bundle.  At  the  highest  heat  fluxes 
available  (40  kW/m2)  there  appears  to  be  little  decrease  in  the  bundle 
performance  (when  compared  to  pure  R-114)  due  to  the  oil.  This  indicates 
that  at  typical  evaporator  operating  heat  fluxes,  the  presence  of  oil  does 
not  significant  effect  the  heat  transfer  enhancement  process.  At  higher 
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heat  fluxes,  however,  the  effect  of  oil  appears  to  be  very  significant  as 
seen  in  Figure  53. 

Figures  55  to  60  show  test  two  to  test  six  for  decreasing  heat 
flux  with  10%  oil.  At  low  heat  fluxes,  there  seems  to  be  no  effect  of  the 
oil  on  the  local  heat  transfer  performance.  However,  at  high  heat  fluxes, 
there  is  a  significant  decrease  in  performance.  Interestingly,  if  one 
compares  Figures  53,  and  Figures  55  to  58  at  high  heat  fluxes,  the  lowest 
activated  tube  in  the  bundle  is  significantly  degraded.  The  effect  of 
activating  a  lower  tube  significantly  enhances  the  heat  transfer  from  the 
tube  directly  above  and  (to  a  lesser  degree)  the  tubes  even  higher  in  the 
bundle .  This  may  be  due  to  the  vigorous  boiling  action  of  lower  tubes 
partly  'scouring'  the  oil  rich  layer  which  'blankets'  the  upper  tubes. 
This  effect  was  also  noticeable  with  the  High  Flux  bundle  (Akcasayar  [Ref. 
25]). 

Figure  61  compares  tube  one  for  tests  one  to  seven  for  decreasing 
heat  flux.  If  one  compares  Figures  26  (0%),  43  (3%),  52  (6%)  and  61  (10%) 
for  tube  one  for  all  seven  tests,  it  is  clear  that  at  low  heat  fluxes,  the 
heat  transfer  coefficient  is  similar,  regardless  of  oil  concentration. 
Furthermore,  convective  effects  are  consistent  and  provide  similar 
enhancements  in  heat  transfer  performance  for  all  concentrations.  At  high 
heat  fluxes  at  0,  3,  and  6%  oil  concentrations,  there  is  little 
enhancement  due  to  activation  of  lower  tubes.  However,  at  10%,  there  does 
appear  to  be  a  small  heat  transfer  enhancement  due  to  activation  of  lower 
tubes.  This  was  attributed  above  to  increased  'scouring'  of  the  oil  from 
the  vicinity  of  the  Turbo-B  surface  by  the  increase  in  bubble  activity  as 
more  tubes  are  activated  within  the  bundle.  However,  for  a  practical 
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operating  heat  flux  range  between  15  and  30  kW/m2,  there  is  no  significant 
degradation  in  heat  transfer  performance  for  an  oil  concentration  of  up  to 
10%. 


E.  COMPARISON  OF  R-114/0IL  MIXTURE  EXPERIMENTS 

Figures  62  and  63  compare  tube  one  from  test  one  for  increasing  and 
decreasing  heat  flux  for  all  oil  concentrations.  Figure  62  shows  no 
effect  of  oil  in  the  NC  region,  but  some  degradation  in  the  boiling 
region.  The  correlation  of  Churchill/Chu  is  included  for  comparison.  The 
incipient  point  appears  relatively  random,  indicating  no  early  or  delayed 
nucleation  caused  by  the  presence  of  oil.  Figure  63  shows  similar 
degradation  with  a  significant  effect  of  the  oil  (20%  decrease  in  the  heat 
transfer  from  0%  to  10%  oil)  at  the  highest  heat  fluxes. 

Figures  64  and  65  compare  tube  one  from  test  seven  for  increasing 
and  decreasing  heat  flux  for  all  oil  concentrations.  As  compared  to  tube 
one  test  one  (Figures  62  and  63) ,  the  presence  of  oil  has  some  degradation 
effect  on  the  heat  transfer  performance  (15%)  in  the  NC  region;  this  may 
be  due  to  a  change  in  the  mixture  properties  which  would  tend  to  increase 
the  wall  superheat  slightly  as  shown.  As  with  Figure  62,  Figure  64  shows 
that  the  incipient  point  appears  relatively  random.  Figure  65  shows 
similar  trends  to  Figure  63  (ie.  no  effect  of  oil  at  high  heat  fluxes), 
but  shows  convective  effects  at  low  heat  fluxes.  This  was  expected  and 
previously  reported  (Figures  28,  43,  52,  and  61). 

Figures  66  and  67  show  the  average  bundle  heat-transfer  coefficient 
(ie.  an  average  of  all  five  instrumented  tubes)  as  a  function  of  heat  flux 
for  test  seven  for  increasing  and  decreasing  heat  flux  respectively  at  all 
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oil  concentration.  The  data  are  from  the  same  data  set  as  that  shown  in 


Figures  64  and  65.  Comparing  Figures  64  and  66,  degradation  is  seen  in 
the  NC  region  (15%)  due  to  the  change  in  mixture  properties  as  mentioned 
above.  However,  at  a  practical  operating  heat  flux  range  between  15  and 
30  kW/m2,  the  presence  of  up  to  10%  oil  causes  no  degradation  in  bundle 
performance  as  seen.  These  trends  over  this  heat  flux  range  were  similar 
to  that  found  for  a  High  Flux  tube  bundle  (Akcasayar  [Ref.  25]). 
Comparing  Figures  65  and  67,  similar  trends  (ie.  no  significant  effect  of 
oil)  are  found.  However,  due  to  the  limit  in  maximum  controllable  heat 
flux  for  test  seven,  data  at  'higher'  heat  fluxes  (  up  to  105  W/m2)  could 
not  be  obtained  (as  mentioned  earlier).  As  shown  in  Figures  62  and  63, 
there  may  be  a  significant  degradation  in  the  heat  transfer  performance  at 
these  higher  heat  fluxes. 

F.  COMPARISON  WITH  PREVIOUS  NPS  DATA 

Figure  68  shows  a  comparison  between  the  present  data  for  a  Turbo-B 
bundle,  the  data  of  Anderson  [Ref.  13]  for  a  smooth  tube  bundle  and  the 
data  of  Akcasayar  [Ref.  25]  for  both  a  19  fpi  and  High  Flux  tube  bundle  in 
R-114.  For  clarity,  only  test  one  (tube  one)  for  a  decreasing  heat  flux 
in  pure  R-114  has  been  shown.  Figure  68  shows  that  the  Turbo-B  tube  has 
a  significantly  lower  heat-transfer  coefficient  than  the  High  Flux  tube  at 
all  heat  flux.  This  is  surprising  since  Sugiyama  [Ref.  34]  showed  that  in 
the  single  tube  apparatus,  the  Turbo-B  tube  was  the  best  performer.  The 
reason  for  this  difference  in  behavior  is  not  known.  Also,  all  of  the  data 
appear  to  be  parallel  to  each  other;  one  may  expect  the  enhanced  tubes  to 
have  a  different  slope  to  a  smooth  tube  due  to  the  greater  amount  of 
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nucleation.  At  the  highest  heat  fluxes,  Turbo-B  and  the  19  fpi  tube 
appear  to  have  a  similar  heat  transfer  performance.  However,  the  finned 
tube  heat  flux  is  based  on  the  root  diameter;  if  the  actual  finned  area 
had  been  used,  then  the  heat  flux  would  be  significantly  lower.  Figure  68 
shows  that  the  heat  transfer  enhancement  given  by  the  Turbo-B  tube  when 
compared  to  a  smooth  tube  is  about  three  at  high  heat  fluxes  and  increases 
to  about  five  at  low  heat  fluxes. 

Figures  69  and  70  compare  the  average  bundle  heat-transfer  coefficient 
(ie.  test  seven)  for  a  given  oil  percentage  to  that  with  no  oil  for  all 
four  tube  bundles  at  heat  fluxes  of  15  and  30  kW/m2  respectively.  15  and 
30  kW/m2  were  chosen  as  being  representative  of  the  lower  and  upper  limits 
of  heat  fluxes  used  in  practical  Naval  evaporators.  At  15  kW/m2,  Figure 
69  shows  large  enhancements  for  the  smooth  and  finned  tube  bundles  for  all 
oil  concentrations,  especially  at  the  lower  oil  concentrations.  However, 
the  Turbo-B  and  High  Flux  tube  bundles  show  a  degradation  in  the  heat- 
transfer  coefficient  at  all  oil  concentrations  (approximately  5-10%  at  low 
oil  concentrations  dropping  to  nearly  25%  at  10%  oil  for  the  High  Flux 
bundle).  At  the  higher  heat  flux  (30  kW/m2) ,  Figure  70  shows  similar 
trends  as  Figure  69.  However,  the  High  Flux  bundle  now  exhibits  a  40% 
decrease  in  the  average  bundle  heat-transfer  coefficient  at  10%  oil. 

Figures  71  and  72  compare  the  average  bundle  heat-transfer  coefficient 
(ie.  test  seven)  for  each  enhanced  tube  to  that  for  the  smooth  tube  bundle 
(tested  by  Anderson  [Ref.  13])  for  all  oil  concentrations  at  heat  fluxes 
of  15  and  30  kW/m2  respectively.  At  15  kW/m2  with  pure  R-114,  Figure  71 
shows  an  enhancement  factor  of  3.7  for  the  Turbo-B  tube  bundle.  This 
enhancement  decreases  slowly  with  increasing  oil  percentage  to  a  factor  of 
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about  2.5  at  10%  oil.  This  agree  very  closely  with  the  19  fpi  bundle. 
The  High  Flux  bundle  exhibits  much  larger  enhancements,  from  over  6  at  0% 
oil  to  just  over  3  at  10%  oil.  At  30  kW/m2  with  pure  R-114,  Figure  72 
shows  an  enhancement  of  3.8  for  the  Turbo-B  tube  bundle,  decreasing  to 
about  3  at  10%  oil.  This  again  agrees  closely  with  the  19  fpi  bundle.  It 
should  be  noted  that  the  High  Flux  bundle  enhancement  has  decreased  to  a 
similar  value  and  gets  worse  than  the  other  bundles  as  the  heat  flux  is 
further  increased. 
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Figure  12.  Performance  of  Test  One  For  Preliminary  Experiments 
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Figure  13.  Performance  of  Test  One  at  Various  Tube  Positions  for 
Increasing  Heat  Flux 
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Figure  15.  Performance  of  Tubes  1  and  2  for  Increasing  Heat  Flux  in  Pure 
R-114 
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Figure  16.  Performance  of  Tubes  1,  2,  and  3  for  Increasing  Heat  Flux  In 
Pure  R-114 
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Figure  18.  Performance  of  All  Five  Tubes  for  Increasing  Heat  Flux  in 
Pure  R-114 
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Figure  19.  Performance  of  All  Five  Tubes  with  Active  Pairs  for 
Increasing  Heat  Flux  in  Pure  R-114 
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Figure  20.  Performance  of  the  Bundle  with  Simulation  Heaters  for 
Increasing  Heat  Flux  in  Pure  R-114 


67 


100000 


o 

o 


1- 

< 

CO 


< 

£ 

h 


Figure  21.  Comparison  of  Tests  One  to  Seven  for  Tube  1  for  Increasing 
Heat  Flux  in  Pure  R-114 
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Figure  23.  Performance  of  Tubes  1,  2,  and  3  for  Decreasing  Heat  Flux  in 
Pure  R-114 
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Figure  24.  Performance  of  Tubes  1,  2,  3,  and  4  for  Decreasing  Heat  Flux 
in  Pure  R-114 
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Figure  25.  Performance  of  All  Five  Tubes  for  Decreasing  Heat  Flux  in 
Pure  R-114 
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Figure  26.  Performance  of  All  Five  Tubes  with  Active  Pairs  for 
Decreasing  Heat  Flux  in  Pure  R-114 
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Figure  27.  Performance  of  the  Bundle  with  Simulation  Heaters  for 
Decreasing  Heat  Flux  in  Pure  R-114 
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Figure  28.  Comparison  of  Tests  One  to  Seven  for  Tube  1  for  Decreasing 
Heat  Flux  in  Pure  R-114  6 
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Figure  30.  Performance  of  the  Bundle  with  Simulation  Heaters  for 
Increasing  Heat  Flux  In  R-114  with  1%  Oil 
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Figure  33.  Performance  of  the  Bundle  with  Simulation  Heaters  for 
Increasing  Heat  Flux  in  R-114  with  2%  Oil 
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Figure  34.  Performance  of  the  Bundle  with  Simulation  Heaters  for 
Decreasing  Heat  Flux  in  R— 114  with  2%  Oil 
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Figure  35.  Performance  of  Tube  1  for  Increasing/Decreasing  Heat  Flux  in 
R-114  with  3%  Oil 
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Figure  36.  Performance  of  the  Bundle  with  Simulation  Heaters  for 
Increasing  Heat  Flux  in  R-114  with  3%  Oil  Heaters  for 
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Figure  40.  Performance  of  All  Five  Tubes  for  Decreasing  Heat  Flux  in  R- 
114  with  3%  Oil 
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Figure  41.  Performance  of  All  Five  Tubes  with  Active  Pairs  for 
Decreasing  Heat  Flux  in  R-114  with  3%  Oil 
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Figure  42.  Performance  of  the  Bundle  with  Simulation  Heaters 
Decreasing  Heat  Flux  in  R-114  with  3%  Oil 
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Figure  43.  Comparison  of  Tests  One  to  Seven  for  Tube  1  for  Decreasing 
Heat  Flux  in  R-114  with  3%  Oil 
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Figure  44.  Performance  of  Tube  1  for  Increasing/Decreasing  Heat  Flux  in 
R— 114  with  6%  Oil 
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Figure  45.  Performance  of  the  Bundle  with  Simulation  Heaters 
Increasing  Heat  Flux  in  R-114  with  6%  Oil 
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Figure  46.  Performance  of  Tubes  1  and  2  for  Decreasing  Heat  Flux  in  R- 
114  with  6%  Oil 
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Figure  48 .  Performance  of  Tubes  1  9  ^  9n>i  /  _ 

in  R-114  with  6%  Oil  -  •  3.  and  4  for  Decreasing  Heat  Flux 
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Figure  49.  Performance  of  All  Five  Tubes  for  Decreasing  Heat  Flux  in  R- 
114  with  6%  Oil 
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Figure  50.  Performance  of  All  Five  Tubes  with  Active 
Decreasing  Heat  Flux  in  R-114  with  6%  Oil 
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Figure  51.  Performance  of  the  Bundle  with  Simulation  Heaters  for 
Decreasing  Heat  Flux  in  R— 114  with  6%  Oil 
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Figure  53.  Performance  of  Tube  1  for  Increasing/Decreasing  Heat  Flux  in 
R-114  with  10%  Oil 
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Figure  54.  Performance  of  the  Bundle  with  Simulation  Heaters  for 
Increasing  Heat  Flux  in  R-114  with  10%  Oil 
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Figure  56.  Performance  of  Tubes  1,  2,  and  3  for  Decreasing  Heat  Flux  in 
R-114  with  10%  Oil 
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Figure  57.  Performance  of  Tubes  1,  2,  3,  and  4  for  Decreasing  Heat  Flux 
In  R-114  with  10%  Oil 
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Figure  58.  Performance  of  All  Five  Tubes  for  Decreasing  Heat  Flux  in  R- 
114  with  10%  Oil 
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Figure  59.  Performance  of  All  Five  Tubes  with  Active  Pairs  for 
Decreasing  Heat  Flux  in  R-114  with  10%  Oil 
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Figure  60.  Performance  of  the  Bundle  with  Simulation  Heaters  for 
Decreasing  Heat  Flux  in  R-114  with  10%  Oil 
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Figure  61.  Comparison  of  Tests  One  to  Seven  for  Tube  1  for  Decreasing 
Heat  Flux  in  R-114  with  10%  Oil 
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Figure  62.  Comparison  of  Test  One  for  Increasing  Heat  Flux  in  R-114  /Oil 
Mixtures 
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Figure  64.  Comparison  of  Tests  One  to  Seven  Tube  One  for  Increasing  Heat 
Flux  in  R-114/Oil  Mixtures 
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Figure  66.  Mean  Bundle  Heat-Transfer  Coefficient  for  Increasing  Heat 
Flux  in  R-114/Oil  Mixtures 
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Figure  67.  Mean  Bundle  Heat-Transfer  Coefficient  for 
Flux  in  R— 114/Oil  Mixtures 
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Figure  68.  Test  One  Comparison  of  Turbo-B,  Smooth,  Finned,  and  High  Flux 
Tube  Bundles  for  Decreasing  Heat  Flux  in  Pure  R-114 
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Figure  69.  Ratio  of  Mean  Bundle  Heat-Transfer  Coefficient  With  Oil  to 
Heat-Transfer  Coefficient  Without  Oil  for  Different  Oil  Percentages  at 
a  Heat  Flux  of  15  kW/m2 
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Figure  70.  Ratio  of  Mean  Bundle  Heat-Transfer  Coefficient  With  Oil  to 
Heat-Transfer  Coefficient  Without  Oil  for  Different  Oil  Percentages  at 
a  Heat  Flux  of  30  kW/m2 
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Figure  71.  Ratio  of  Mean  Bundle  Heat-Transfer  Coefficient  of  Enhanced 
Tube  to  Heat-Transfer  Coefficient  of  Smooth  Tube  for  Different  Oil 
Percentages  at  a  Heat  Flux  of  15  kW/m2 
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Figure  72.  Ratio  of  Mean  Bundle  Heat-Transfer  Coefficient  of  Enhanced 
Tube  to  Heat-Transfer  Coefficient  of  Smooth  Tube  for  Different  Oil 
Percentages  at  a  Heat  Flux  of  30  kW/o2 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

Nucleate  boiling  data  of  R-114  at  atmospheric  pressure  were  obtained 
using  a  small  bundle  of  Turbo-B  copper  tubes.  The  data  were  obtained  for 
both  increasing  and  decreasing  heat  flux  and  at  different  oil 
concentrations.  Based  upon  the  results  pertaining  to  this  particular 
bundle  and  apparatus,  the  following  conclusions  may  be  made: 

1.  Natural  Convection  Region 

a.  For  a  single  upper  tube,  a  second  lower  tube  directly  below 
when  turned  on  does  increase  the  heat  transfer  performance  of  the  upper 
tube,  however  when  additional  lower  tubes  are  heated  no  net  increase  in 
performance  occurs. 

b.  The  presence  of  heated  lower  tubes  in  the  bundle  reduces  the 
incipient  boiling  point  of  the  upper  tubes  and  the  tubes  tend  to  nucleate 
'in  order'  (ie.  top  tube  first,  bottom  tube  last). 

c.  The  effect  of  adding  oil  to  the  refrigerant  (up  to  10%) 
reduces  the  heat-transfer  coefficient  slightly  (approximately  10-15%)  due 
to  changes  in  the  fluid  properties. 

2.  Boiling  Region 

a.  For  pure  R-114,  the  presence  of  heated  lower  tubes  on  the  top 
tube  causes  no  enhancement  at  high  heat  fluxes  (>  20  kW/m2)  ,  but  at  low 
heat  fluxes  (<  20kW/mz) ,  there  is  a  significant  enhancement  due  to 
convective  effects. 


120 


b.  At  very  low  heat  fluxes  (<2  kW/m2)  ,  the  presence  of  oil  has 
little  effect  on  the  heat  transfer  performance  of  the  top  tube  in  the 
bundle.  At  higher  heat  fluxes  (>  2kW/m2)  ,  the  performance  is  enhanced  by 
10-15%  at  low  concentrations,  but  is  degraded  up  to  20%  at  10%  oil 
concentration  at  the  highest  heat  fluxes. 

c.  At  typical  operating  heat  fluxes  (15-30  kW/m2) ,  the  bundle 
performance  is  reduced  between  5-15%  with  oil. 

B.  RECOMMENDATIONS  FOR  FUTURE  WORK 

1.  Conduct  experiments  with  varying  pool  height,  but  keep  the  local 
pressure  at  each  tube  constant  by  simultaneously  varying  the  vapor 
pressure  above  the  pool. 

2.  Additional  experiments  with  R-113  and  R-114  should  be  conducted  to 
investigate  explosive  (R-114)  and  partial  (R-113)  incipience  at  the  onset 
of  nucleation  varying  the  time  at  the  incipience. 

3.  Some  instrumentation  should  be  added  such  that  the  flowrates 
through  the  bundle  can  be  determined.  From  these  measurements,  vapor 
quality  can  be  determined. 

4.  Metal  guide  plates  should  be  manufactured  and  placed  on  each  side 
between  the  simulation  tube  bundle  and  the  tube  bundle  itself.  This 
further  channels  the  flow  of  refrigerant  thru  the  bundle  at  high  heat 
fluxes . 

5.  Attention  needs  to  be  given  to  the  question  of  refrigerant 
disposal.  There  are  reclamation  projects  undertaken  by  most 
manufacturers;  however,  a  method  still  needs  to  be  found  to  remove  the 
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refrigerant  from  the  apparatus  into  a  container  suitable  for  such 
reclamation. 

6.  A  high  speed  camera  should  be  used  to  study  the  nucleation  process 
and  circulation  patterns  in  more  detail  in  the  bundle.  Neutrally  buoyant 
particles  might  be  placed  in  the  pool  to  facilitate  study  of  circulation 
patterns  within  the  bundle. 
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APPENDIX  A:  LIST  OF  DATA  FILE 


Table  3.  DATA  FILE  NAMES  FOR  TURBO-B 
TUBE  BUNDLE  EXPERIMENTS 


FILE  NAME  NUMBER  OF  NUMBER  OF 


PERCENT  NUMBER  OF  NUMBER  OF 


TBI0001H 


DATA  POINTS  INSTRUMENTED  TUBE  OF  OIL  ACTIVE  SIMULATION 

PAIRS  HEATERS 


FILE  NAME 

TBI0007 

TBD0007 

TBI0101 

TBD0101 

TBI0107 

TBD0107 

TBI0201 

TBD0201 

TBI0207 

TB00207 

TBI0301 

TB00301 

TB00302 

TBD0303 

TB00304 

TB00305 

TB00306 

TBI0307 

TB00307 


TBI0601 


32 


1 
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FILE  NAME 

NUMBER  OF 

DATA  POINTS 

NUMBER  OF 
INSTRUMENTED  TUBE 

PERCENT 

OF  OIL 

NUMBER  OF 
ACTIVE 

PAIRS 

NUMBER  OF  1 
SIMULATION  | 
HEATERS  1 

TBI0001E 

33 

1 

0 

0 

0 

TBI0001F 

38 

1 

0 

0 

0 

TB 1000 1C 

38 

1 

0 

0 

0 
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APPENDIX  B:  SAMPLE  CALCULATIONS 


Data  set  number  1  Tube  1  of  experiment  TBD1005  (Turbo-B  tube, 
decreasing  heat  flux,  10%  oil  concentration,  test  5)  was  used  for  the 
sample  calculations  in  order  to  validate  the  program  used  for  data 
acquisition  DRP4RH.  The  working  fluid  was  R-114. 

1.  Test  tube  dimensions 

Dtc  -  11.60  mm 
Dc  -  14.15  mm 
D,-  —  12.70  mm 
L  -  203.2  mm 
Ly  -  25 .4  mm 

2 .  Measured  Parameters 


T1  - 

10. 

.62 

°C 

T2  - 

10. 

.94 

°C 

T3  - 

9. 

.96 

°c 

T4  - 

10. 

.90 

°c 

T5  - 

10. 

,04 

°c 

T6  - 

9. 

.07 

°c 

Tldl 

/ 

■  4 

>.27 

'  °c 

Tld2 

* 

■  4 

>.21 

.  °c 

Aas  -  3 . 

.513 

1  V 
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Vas  -  3.189  V 


3.  Calculations 

The  heaters  power  is  first  calculated  for 
q  -  Vas(V)xAas(V)x60(V/V)xl(A/V) 

Note:  The  multiplication  factors  of  volts  and  amp  sensors  are  60  and 
1,  respectfully. 

Therefore : 

q  -  (3 . 189) (3.513) (60V/V) (1A/V) 

q  -  672.19  Watts 

The  tube  inside  wall  temperature  is  obtained  from  the  average  of  all 
six  thermocouple  readings. 

6 

E  Tn 

6 

-  1/6(10.62  +  10.94  +  9.96  +  10.90  +  10.04  +  9.07) 

-  10.25  °C 

The  tube  outside  temperature  is  calculated  by  knowing  the  inside  wall 
temperature  using  Fourier's  Conduction  Law.  Uniform  radial  conduction  is 
assumed. 

rp  •  T  _ UtC 

wi~  2it  (kcu)  (L) 

where  the  second  term  on  the  right  hand  side  is  the  Fourier  conduction 
term.  If  we  define  this  term  as 


131 


2  it  {k^)  (L) 


and 

O*  ”  Two  -  Tsatc 

where  is  the  thermal  conductivity  of  copper  and  is  calculated 

follows 


kcu  -  434.0  -  [0 . 112  ( Tvi)  ] 


kcu  -  434.0  -  [0.112(283.25)] 


kcu  -  402.28  W/mK 


now 


672.15  [In  ( 44-4I- )  ] 


jt  B  t  _  11.60 

"  vi  2  It  (402 . 28)  (.2032 


Tuo  -  (10.25  -  .2601)  0  C 


Tvo  -  9.98*  C 


The  liquid  saturation  temperature  at  the  top  of  the  tube  bundle  is 


as 
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Tsat  - 


tldl  +  Tld2 
2 


Tsat  - 


2.27  +  2.21 
2 


Tsat  -  2.24 °C 

In  order  to  calculate  the  local  saturation  temperature  for  each  tube, 
correction  factors  are  needed  to  account  for  hydrostatic  pressure 
differences  between  the  tube  locations  and  the  liquid  free  surface.  This 
difference  is  calculated  by: 

Ap  -  p  ( g )  (ht) 

For  the  top  tube  in  the  bundle  which  is  0.124  m  below  the  thermocouple 
measuring  pool  temperature. 

AP  -  1523.12  (9.81)  (0.124) 

AP  -  1852.78  Pa 

For  1852.78  Pa  pressure  difference,  corrected  saturation  temperature  is 
obtained  by  adding  0.04  °C  (from  standard  tables  for  R-114)  to  Tsat. 
Corrected  Tsat  is: 


Tsa te  -  (2.24  +  0.04)  °  C 
Tsa tc  -  2.28  0  C 
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Therefore,  the  wall  superheat  can  be  obtained  by  the  following: 

Bb  -  Two  -  Tsatc 

9b  -  (9.98  -  2.28 )  °  C 

0i-7.7O°C 

Now  that  the  wall  superheat  is  known,  we  need  to  calculate  the  heat 
flux  and  the  heat-transfer  coefficient.  To  do  this,  we  know  that  the  tube 
is  12  inches  long  and  is  heated  in  a  eight  inch  center  portion  of  the 
tube.  The  unheated  lengths  of  the  tube  are  a  one  inch  and  a  three  inch 
section  on  opposite  ends  of  the  tube.  These  unheated  lengths  have  a  fin 
effect  during  the  heat  transfer  process  to  the  evaporating  refrigerant. 
In  order  to  account  for  this,  the  following  procedure  was  adopted  for  both 
one  and  three  inch  sections.  Calculations  are  shown  below  for  the  one 
inch  section.  Heat  transfer  from  the  unheated  end  is  calculated  as  heat 
from  the  base  of  the  fin: 

qt  -  [<hb)  (p)  (k^)  (Ac)]°-i(eb)  (tanh  [  (n)  (Le)J  ) 
where 

P  -  *  (D0) 

-  n  (.01415)  m 

-  .04445  m 

now 

Ac  -  w/4(D20  -  DZj ) 

-  jr/4(  . 014152  -  . 012702) 

-  3.0578  x  10‘5  m2 
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The  corrected  length  of  unenhanced  surface  at  the  end  was  calculated  as 
follows 

4  -  Lu  +  (t/2) 

-  0.0254  +  [(0. 01415-0. 0127)/2] 

-  0.0258  m 

hb  is  the  natural  convection  heat  transfer  coefficient  of  the  fin  like 
ends  and  was  calculated  by  using  Churchill-Chu  [Ref.  22]  correlation  for 
natural  convection  for  a  smooth  cylinder,  as  modified  by  Pulido  [Ref. 27], 


h*  - 


t10 


[ 


6+ . 387 


[gr< P>  {Dl)  (0b)  (tanh(nLc) )  ] 
v  (a)  UJ  (n) 


1/6 


Cl*  l- 


.559 

Pi 


]9/l6] 


S/27 


where 


n 


,  lhb)(P)  °'5 
<*«)  <AC) 


Therefore  an  iterative  technique  was  necessary  to  calculate  hb.  The 
iterative  technique  used  was  to  assume  hb  was  190  W/m2K  and  continue  the 
iteration  until  successive  values  are  within  0.001  of  each  other.  The 
fluid  physical  properties  are  calculated  at  the  vapor  mean  film 
temperature,  given  by  the  following  equation. 


Tlilm 


Tsatc  +  Two 


2 
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^  _  2.28+9.98 

1  film  "  *5 


Tfii„  -  6 . 13  °  C  -  279 .13  °  K 


For  R-114,  the  physical  properties  are  given  in  the  program  by 
Dynamic  viscosity.  Tfnm  in  °K 

/i  -  exp [-4. 4636  +  (1011 ,47/Tfilm)  JxlO"3 
fi  —  430.927xl0~6  kg/m  s 
Specific  heat.  in  °K 

Cp  -  [0.40188  +  1.65007xl0'3(Tmm)  +  1 . 51494xlO~6(Tf  Uni2)  - 
6 . 67853xlO-10(Tf nm3)  ]xl03 
Cp  -  966.31  J/kgK 
Density,  Tfilm  in  °K 

p  -  16.0184533  (36.32  +  61 . 146414<|i1/3  + 

17 .47  6  838i|r1/2  «*•  1 . 119828i|r2) 


where 


*  -  1 


[1. 8 (Tfilm) ] 
753.95 


and 


♦  -  .333 

p  -  1512.09  kg 
m3 

Thermal  conductivity  of  R-114,  Tf j in  °C 
k  -  0.071  -(0.000261) (TfUffl) 
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k  -  6.936xlO-2  W/mK 


Prandtl  Number 

Pr  -  [  (Cp)/i]/k 
Pr  -  6.003 

Thermal  Expansion  Coefficient 
P  -  -  (1/p)  (Ap/A T) 

p279.03  -  15  1  2.3  95 

m3 

P279.21  “  1511.824  -^2 
m3 

P  -  -(1/1512.395) [(.571)/ (0.2)] 

P  -  1.89 xlO'3  (1/JO 

Kinematic  viscosity 


v  -  -t 
P 

v  -  430-927 xlO'* 
1512.09 

v  -  2.849X10'7  m2/s 

Thermal  Diffusivitv 
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a  - 


k 

(P  )cp 

6 . 9  36x1 0~2 
(1512.09)966.31 

a  -  4  .747x10'®  m2/s 

Knowing  the  above  properties,  the  heat-transfer  coefficient  hb,  can  be 
obtained  by  iteration 

hb  -  362.57  W/m2K 

Knowing  this  we  can  calculate  n 

n  ,  (P) 

1  <*„)  Uc) 


n  _  ,  (362.57)  (44.45xl0~3)  j0'5 

1  (402.28)  (3.0578X10'6) 

n  -  36.19 

then  we  can  obtain  qf 

qf  -  [(hb)  (P)  (*„,)  (Ac)]°-5(6i>)  (tanh  [  (n)  (Lc)  ]  > 

qf  -  [(362.57)  (.04445)  (402.28)  (3 . 0578X10-5) )  0  5 
(7.70) (tanh [(36. 196) (0.0258)] ) 
qr,  -  2.51  W 

The  corresponding  results  for  the  three  inch  section  are 
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hb  -  289.47  W/m2K 
qf  -  1.24  W 


Therefore,  the  heat  transfer  through  the  heated  length  of  the  tube  is 

9s  “  3  “  qf  (1  inch  section)  -  qf  (3  inch  section) 
qs  -  (672.19  -  1.24  -  2.51)  W 
qs  -  668.35  W 

and  the  heat  flux  and  the  heat  transfer  coefficient  are  as  follows 

q"  -  qs/As 

-  qs/((*)<Do)(L)) 

-  (668. 35)/((ir)(0.01415)(.  2032)) 

-  7 . 398x10*  W/m2 

and  finally  the  heat  transfer  coefficient 


A,(TW  -  Tsat) 


_ 668,35 

9 . 033x1  O'3  (7 .70) 


h  -  9.6OSX103  W/m7K 
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APPENDIX  C:  UNCERTAINTY  ANALYSIS 


The  same  data  run  (TBD1005)  was  chosen  for  the  uncertainty  analysis. 
Therefore,  the  measured  and  calculated  parameters  found  in  the  sample 
calculation  were  used  in  this  section.  The  unceratainty  analysis 
performed  was  for  a  high  heat  flux,  but  the  procedure  could  be  performed 
at  any  heat  flux  to  determine  the  uncertainty  bands.  All  uncertainties 
are  presented  as  a  percentage  of  the  calculated  parameter.  The 
uncertainty  associated  with  the  experimental  parameters  is  calculated  from 
the  equation  suggested  by  Kline  and  McClintock  [Ref.  35].  For  example: 

R  -  R(xltx2 . xn) 

then 

bR  ’  (<-|r8xi)2+(-irfixJ)2+ . +  <-|J-**n)ai0-5 

where 

5R  -  uncertainty  of  the  desired  dependant  variable 
xn  -  measured  variables 
6xn  -  uncertainty  in  measured  variables 
The  boiling  heat-transfer  coefficient  is  given  by 

h - ** _ 

At(Tvo-Tsat) 

where 


<7(ln  ( )  ] 
"  Twi~  2x(kcll)t(L) 
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In  the  above  equation,  the  second  term  on  the  right  hand  side  is 
usually  called  the  Fourier  heat-transfer  conduction  term.  we  define 
this  as 


Qrtln(-^p-)  ] 

_ UtC 

2n  (k^)  (L) 


and 

6„  -  Tw-Tsatc 

With  this  notation,  the  uncertainty  in  the  heat-transfer  coefficient 
is  obtained  using  the  following  equation. 

asat:)2]0-5 

h  g  As  Qb  0b 

where 
q-  Vxl 

q-  V(V)xI(V)x60(V)xl(A/V) 
and  the  uncertainty  is 


bg 

Q 


[<^V- 


The  accuracy  in  the  voltage  and  current  sensors  are  as  follows 
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tVas  -  ±0.05  V 


t>Aas  -  ±0.025  A 


From  the  sample  calculation  section 

Vas  -  3.189  V  Aas  -  3.513  V 

Therefore , 


6<?  .  [  ( A¥aS)2+  ( AA££)2]°5 

g  Vas  Aas  1 


A3  .  U  °-05  (  0-025, 210. 5 

q  3.189  3.513 


Aa  -  1.72  percent 
Q 


Calculation  of  the  uncertainty  of  the  surface  area  is  as  follows 
^  -  *(D0)(L) 


0.5 


Knowing  the  dimensions  of  the  tube  from  the  manufacturer  and  estimated 
inaccuraccies  from  work  shop  tools  and  human  error,  the  unceratinty  was 
calculated. 

Dimensions 

D0  -  14.15  mm  L  -  203.2  mm 

Inaccuracies  in  measurements 
6D0-  0.1mm  6L  -  0.2  mm 

Uncertainty  analysis  performed 
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Ml  _  [  (Ml)2*  (A£:)  *]  o  s 


Ml  -  [  (  . -JLi—  )2+  (  °-2  )  2]  0-5 

A,  14.15  203.2 


0.7135  percent 


The  uncertainty  calculation  for  the  Fourier  conduction  term  given  below 


_  f  (  &<7\  j2f  i  &  4  \  21  0.5 

4 >  <1  *cu  L 

kcu  was  calcluated  using 


-  434.0-  [0.112  (Tvi)] 
k, m  -  434.0- [0.112(283.25)] 
km  -  402.28  W/mK 

and  its  uncertainty 

bkcu  -  [  (0.112  llTvi))2]0i 

6TMj  and  fiTsat  are  obtained  using  uncertainties  in  the  thermocouple 
readings.  Average  wall  inside  temperature  TMj  was  obtained  taking  the 
average  of  six  thermocouple  readings  inside  the  tube  wall.  The 
uncertainty  associated  with  this  variable  is 

&Twi  -  [6  (  )  ]«•» 

O 
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where  6Twj  for  each  thermocouple  was  obtained  by  taking  the  difference 
between  the  measured  wall  temperature  and  the  average  wall  temperature. 
Using  this  method,  it  has  been  attempted  to  try  and  take  into  acount  some 
the  uncertainty  introduced  by  the  fabrication  procedure  for  the  tube  (ie. 
air  gap)  .  For  this  particular  heat  flux  the  following  6Twj  were  found  for 
tube  thermocouples  1,  2,  3,  4,  5,  and  6. 


bfwi  - 

o  o 


)  *( 


0.25 

6 


0.75 

6 


)2+( 


0.21 

6 


>2+( 


1.18 

6 


0.5 


brvi  -  .274°  c 


The  uncertainty  level  for  all  remaining  thermocouple  readings  (ie. 
excluding  those  in  the  tube  wall  which  was  considered  to  have  a  higher 
uncertainty)  was  estimated  to  be  +/-  0.1  °C  corresponding  to  an  emf  of 
approximately  4  /*V.  Saturation  temperature  was  obtained  by  taking  the 
average  of  two  thermocouple  readings  and  the  uncertainty  in  this 
temperature  was  calculated  from  the  following  equation. 


&T  2 

bTsat  -  [2  ( — ~ )  ]°  5 

bTsat  - 

a  Tsa t  -  0 . 07  °  C 

Knowing  the  uncertainty  in  the  temperatures,  we  can  now  calculate  the 
uncertaines  in  the  following: 
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bkcu  -  [  (0.112  (&Tvl)  )2]  0  5 
&kcu  -  [<0. 112  (283. 25)>2]0-5 
5Jccu  -  31.724  W/mK 

Now  we  can  calculate  the  uncertainty  in  the  Fourier  conduction  term 


-  [(l2)2+(^£H)2+(-®i;)2]' 
<t>  <1  km  L 


M.  [(0.0172)J.(-|— 


-  8.072  percent 
<P 


from  the  sample  calculations  we  know  that 


<J>  -  0.2679  °  C  6$  -  0.2679(0.08072)  -  0.0216°  C 


therefore : 


ftrsatyjo.s 

h  q  As  0  b  0.,  0* 


-  [(o.0172)2+(0.07135)M44£>2+<-^|^>2+(^£>2]0'5 

h  7.70  7.70  7.70 


-  7  .42  percent 
h 


Finally  the  calculation  of  wall  superheat  temperature 


Qb  -  Tvo  -  Tsat 


i8„  ..87.  1 


b  _  it  VIvo\  j.  /  6rsat j 2] o.5 


e»  -  I(-ef>  *  <-e 
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-  [(  0  ■•!?*.  )*  +  (^07  pio.s 

eb  11  7.7  1  K  7.7 ‘  J 

se. 

-  3.67  percent 

Table  4  shows  the  results  of  the  uncertainty  analysis  performed.  The  high 
and  low  heat  flux  correspond  to  the  approximate  values  of  7.5xlOA  W/m2  and 
lxlO3  W/m2  respectively.  Note  that  the  highest  uncertainty  (over  50%)  is 
in  the  wall  superheat  at  very  low  heat  flux  (1  kW/m2)  .  This  is  due  to  the 
very  low  measured  value  of  wall  superheat  (0.54  °C)  which  can  not  be 
accurately  measured.  Thus  higher  uncertainty  occurs  at  very  low  heat 
fluxes.  However,  once  the  wall  superheat  gets  higher  (at  higher  heat 
fluxes)  the  uncertainty  in  wall  superheat  decreases  significantly  (to 
about  4%)  indicative  of  the  fact  that  the  measure  wall  temperature  is 
relatively  more  accurate. 
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Table  4.  UNCERTAINTY  ANALYSIS  RESULTS 


VARIABLE 

HIGH  HEAT  FLUX 

LOW  HEAT  FLUX 

0* 

7.7 

0.54 

Tvl 

10.25 

2.80 

oh 

2.24 

2.20 

6V as/Vas 

1.57% 

14.6% 

5Aas/Aas 

0.712% 

7.9% 

6q/q 

1.72% 

16.6% 

5D0/D0 

0.707% 

0.707% 

5L/L 

0.098% 

0.098% 

^Ag/Ag 

0.714% 

0.714% 

^^cu/^cu 

7.89% 

7.89% 

60,,/e* 

3.67% 

52.37% 

6h/h 

7.42% 

28.93% 
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APPENDIX  D:  OPERATING  PROCEDURE 


A.  SYSTEM  STARTUP 

1.  Power  to  the  28  kW  (8  ton)  refrigeration  unit  is  provided  by  the 
breakers  located  in  the  main  distribution  panel  located  in  the  laboratory. 
These  breakers  were  never  secured.  However,  if  power  to  this  panel  was 
lost,  then  these  breakers  must  be  reset. 

2.  Turn  the  switch  on  the  refrigeration  unit  control  panel,  located 
in  front  of  the  refrigeration  unit  to  the  "auto"  position  after  passing 
through  "on"  position.  This  switch  is  always  left  on,  unless  unit  was 
taken  down  for  long  repairs . 

3.  Push  the  start  button  in  the  control  box  for  the  recirculation 
pump.  This  control  box  is  located  on  the  bulkhead  above  the  recirculation 
pump  in  the  outside  area  adjacent  to  the  refrigeration  unit. 

4.  Set  the  desired  temperature  on  the  roughly  graduated  Fahrenheit 
scale  on  the  control  panel  thermostat.  It  requires  approximately  one  hour 
to  chill  the  sump  to  -  15  °C.  The  thermometer  in  the  ethylene  glycol/water 
mixture  (sump)  must  be  monitored  to  ensure  the  desired  sump  temperature  is 
attained  and  maintained.  Slight  adjustments  in  the  refrigeration  unit 
thermostat  can  be  expected  due  to  the  coarseness  of  its  scale. 

5.  Energize  the  desired  pumps  by  switching  on  the  breakers  in  the  main 
distribution  power  panel.  Once  the  power  is  energized  in  the  main  power 
panel  and  after  ensuring  the  pump  suction  valves  are  open,  turn  on  the 
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pump  motors  by  pushing  down  on  the  arm  of  the  appropriate  breaker  box  for 
the  pumps  located  on  the  bulkhead  next  to  the  ethylene  glycol/water  sump. 
The  pumps  are  marked  "auxiliary  condenser"  (Pump  #2)  and  "instrumented 
tube  condenser"  (Pump  #1),  respectively  in  the  breaker  box.  Flow  in  the 
auxiliary  condensate  coils  can  be  controlled  with  the  individual  globe 
valves  located  the  coil  penetrations  on  the  apparatus.  The  auxiliary 
condenser  coils  will  produce  the  fastest  adjustments  to  the  system 
pressure . 

6.  Energize  the  heater  variac(s)  desired  by  switching  on  the  breakers 
(Bundle,  and  Simulation  (for  test  7  only))  in  the  main  distribution  power 
panel  and  individual  breakers  for  each  of  these  in  the  power  distribution 
box  (near  apparatus). 

7.  After  ensuring  that  the  breakers  for  the  heated  tubes  desired  are 
in  the  "on"  position,  follow  the  experimental  procedures  for  normal 
operation  outlined  in  Chapter  IV. 

B.  SYSTEM  SHUTDOWN 

1.  Turn  all  variacs  to  the  zero  position  and  switch  off  all  breakers 
in  the  power  panels. 

2.  If  apparatus  will  not  be  operated  for  an  extended  period,  turn  the 
switch  on  the  refrigeration  control  panel  to  the  "off"  position  after 
passing  through  "on". 

3.  Allow  the  recirculation  pump  to  operate  for  at  least  five  minutes 
after  switching  off  the  refrigeration  pump  unit  to  dissipate  any  back 
pressure  in  the  system. 
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4.  Turn  the  breakers  for  the  pumps  to  the  off  position  at  the  switch 


boxes,  and  then  secure  the  power  at  the  main  distribution  power  panel. 


C .  EMERGENCY  SHUTDOWN 

1.  Secure  all  power  at  the  main  distribution  power  panel 

2 .  Evacuate  building 

3.  Call  Fire  Department 
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APPENDIX  E:  PROGRAM  DRF4RH 


10001  file  name  ORP4RH 

1304'  OATE  November  22,  1988 

1008'  REVISED  FES  1992  (R.  HAAS) 

1012' 

1016  BEEP 

1023  PRINTER  IS  1 

1024  IdP‘0 
1023' 

1032  PRINT  USING  *4X,  '"Select  option  default  is  3 

1026  PRINT  'JSIN6  "6x,"0  Taking  data  or  re-drocesaing  previous  data"' 

1040  PRINT  USING  '6X,"1  Plotting  data  on  Leg-Log 

1044  PRINT  USING  '6X,"2  Plotting  data  on  Linear"' 

1048  PRINT  USING  '6X,"3  Purge'" 

1052  PRINT  USIN6  '6X."4  Fixup*" 

1056  PRINT  USING  '6X."S  Move"* 

1063  PRINT  USING  '6X,"6  Comb'" 

1064  PRINT  USING  '6X,"7  Read  Plot"* 

1068' 

1072'  IOP  IS  A  PROGRAM  UARIA8LE  TO  SELECT  A  SUBROUTINE 
1076  INPUT  I  dp 

1080  IF  Idp-0  THEN  CALL  Main 
1084  IF  Idp-1  THEN  CALL  Plot 
1088  IF  Idp-2  THEN  CALL  PI  in 
1092  IF  Idp-3  THEN  CALL  Purg 
1096  IF  Idp-4  THEN  CALL  Fixup 
1100  IF  Idp-S  Then  call  Move 
1104  IF  Idp-6  THEN  CALL  Comb 
1108  IF  Idp-7  THEN  CALL  Readplot 
1112  END 

1116' 

1120  SU8  Main 

1124'  ICAL-THERMOCOUPLE  CALIBRATION 
1128  COM  /Cc/  C <  7  ) 

1132  0 IM  Enf(35),T<35),0la<6),O2a<6),Oia<S ) , Doa< 6  > ,La( 6 ) , Lua< 6  )  ,Kcua( 6  > ,St  < 19), 

Ldte(20),Uolt<2  ) , Amp(  1 1  ),Tua(5),Tul5).Theta(5),Thetab(5),O(5>,0l<5),Qdp(5) 

1136  OIM  Htube(5),Tn(5),Tp(6) 

1140' 

1144'  THERMOCOUPLE  ARRAY  <C<  ))  INITIALIZATION 
1148  OATA  0.10086091,25727.94369,-767345.8295,78025595.81 
1152  OATA  -9247486589, 6. 97688E+11. -2. 56192E+13. 3. 94078E+14 
11S6  REAO  C( •  ) 

1160' 

1164'  PRINT  HEAOER  ANO  INITIALIZE  TIME  CLOCK 
1168  PRINTER  IS  701 
1172  BEEP 

1176  INPUT  'ENTER  MONTH,  OATE  ANO  TIME  ( MM : 00 = HH : MM : SS  )* , Da t e* 

1180'  OUTPUT  OIRECTEO  TO  OATA  AOUISITION  SYSTEM  (HP  3497A) 

1184  OUTPUT  709, 'TO* ;0ateS 
1188  OUTPUT  709. 'TO* 

1192  ENTER  709;Oate* 

1196  PRINT 

1230  PRINT  '  Montn,  date  and  time  -'.Date* 

1204  PRINT 

1208  PRINT  USING  *  10X .  *  *  NOTE  Program  name  0RP4RH”* 

1212  SEEP 
1216' 

1220'  ON  IS  THE  VARIABLE  FOR  OISC  NUMBER  FOR  RECORD  KEEPING  ONLY 

1224  INPUT  'ENTER  DISK  NUMBER', On 

1228  PRINT  USING  '16X,"0isk  number  -  ",2Z*.0n 
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1232  BEEF 
1236  lm-0 

1240  INPUT  "ENTER  INPUT  MODE  < 0-34S7A . 1-FIlE  0-OEFAULT '  ,  Ir 
1244' 

1246'  INPUT  MODE  ZERO  IS  FROM  THE  DATA  AQUISITION  SYSTEM 
12S2  IF  ln-e  THEN 
12S6  BEEP 

1260  INPUT  "6IUE  A  NAME  FOR  THE  RAW  DATA  FILE"  ,D2f lie* 

1264  PRINT  USING  "16K,""File  name  * " . 14A" ; D2f l le* 

1268' 

1272'  CREATE  BOAT  FILE  ON  THE  MASS  STORAGE  MEDIA 
1276  CREATE  BOAT  O2file*,60 

1280'  CREATE  an  INPUT/OUTPUT  LINK  TO  OPEN  FILES 
1284  ASSIGN  #File2  TO  02file* 

1288' 

12321  CREATE  DUMMY  FILE  UNTIL  Nrun  KNOWN 

1296  D1  f  i le*” "DUMMY" 

1300  CREATE  BOAT  Olfile*,60 

1304  ASSIGN  CFilel  TO  Olfile* 

1308  OUTPUT  •Filel'Oete* 

1312' 

1316'  CREATE  A  PLOT  FILE 

1320  BEEP 

1324  INPUT  "6IUE  A  NAME  FOR  THE  PLOT  FILE", Pf lie* 

1328  CREATE  BOAT  Pfile«,30 

1332  ASSI6N  0Plot  TO  Pfile* 

1336  BEEP 

1340' 

1344!  IOTC  -  NUMBER  (TOTAL)  OF  OEFECTIUE  THERMOCOUPLES 
1348  INPUT  "ENTER  NUMBER  OF  DEFECTIVE  TCS  ( 0-OEFAULT > " , I dtc 

1352'  LDTC  -  LOCATION  OF  DEFECTIVE  THERMOCOUPLE 

1356' 

1360  IF  Idtc-0  THEN 

1364  PRINT  USIN6  " 16X , " "No  defective  TCs  exist""" 

1368  ELSE 

1372  PRINT  USING  " 16X ," "Defect ive  Thermocouple:  Indicated  by  -99. 99""" 

1376  END  IF 

1380' 

1384  BEEP 

1388'  DEFECTIVE  THERMOCOUPLES  MAY  BE  IN  CHANNELS  40-E9 

1392'  THEMOCOUPLES  ARE  ENTERED  AS  DEFECTIVE  BY  COMPUTER  CHANNEL  NO. 

1396'  JDTC-COUNTER  IN  LOOP  FOR  DEFECTIVE  THERMOCOUPLES 

1400' 

1404  IF  ldtc>0  THEN 

1408  FOR  Jdtc-0  TO  Idtc-1 

1412  INPUT  "ENTEF  DEFECTIVE  TC  LOCATION  (BY  COMPUTER  CHANNEL  NUMBER) 

" ,Ldtc<  Jdte  ) 

1416  BEEP 

1420  NEXT  Jdtc 

1424  END  IF 

1428  PRINTER  IS  701 

1432  OUTPUT  9F 1 1  el ; Ldte( •  ) 

1436  ' 

1440  I  In* 1  option  (THIS  OPTION  ALLOWS  DATA  ENTRY  WITH  DATA  C'lE’ 

1444  else 

1448  BEEP 

1452  INPUT  "GIVE  THE  NAME  OF  THE  EX1STIN6  DATA  FILE' , D2f ; le* 

1456  PRINT  USIN6  "16x,""File  name  . 14A" ;D2f i le* 

1460  ASSIGN  @File2  TO  DZfiie* 

1464  ENTEP  9F i 1 e2 • Nrun 
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1468  ENTER  $F  i  lei -Odd*  ,Ldtc'.  •  1 , 1 1 t ,8oo  ,Nht  ,Nato  ,Nrt ,  Carr 

14'Z  SEEP 

1 4'S  INPUT  ”  G I UE  A  NAME  FOR  PLOT  FILE'.PfileS 

1480  CREATE  BOAT  P?ile*,30 

1484  ASSIGN  SPlot  TO  PTa«* 

1488  PRINT  USING  -l6X,--This  data  set  taken  on  "  , 14A* .Gold* 

149C  BEEP 

1426  PRINTER  IS  1 

1500  PRINT  USING  "4X , " "SELECT  TUBE  TYFE "  *  “ 

1504  PRINT  USIN6  "S<  , "  "0  SMOOTH"* 

1508  PRINT  USING  -SX,-*1  FINNEOt  1 9/ IN  )  --- 

1512  PRINT  USING  -6X,--:  HIGH  FLUX 

1516  PRINT  USING  "  6  X  ,  -  *  3  TUR80-B  --- 
1520  INPUT  Itt 

1524  END  IF 

1528  IF  In-1  THEN  SOTO  1768 

1532  PRINTER  IS  1 

1536' 

1540  IF  In-0  THEN 


1544 

PRINT 

USIN6  -4X, “Select  tube  type-*- 

1548 

PRINT 

USING  -6X, **  0 

Snooth 

15S2 

PRINT 

USING  -SX  ,  *  *  1 

FINNED  19/IN  < DEFAULT >--* 

1556 

PRINT 

USING  -6X , *  *  2 

HI6H  FLUX--- 

1560 

PRINT 

USING  -6X  ,  *  *  3 

TURBO-B--- 

1564 

PRINT 

USING  -6X,--  4 

GROWTH--- 

1583 

PRINT 

USING  "SX  , -  -  5 

SROWTH--* 

1572 

PRINT 

USIN6  -6X , "  “  6 

GROWTH- - " 

1576' 

ITT-TUBE  TYPE 

1580 

INPUT 

Itt 

1584  OUTPUT  tFilel.Itt 

1588  ENO  IF 
1592  PRINTER  IS  701 
1596'  Itt-2 

1600  PRINT  USIN6  "  16X , ' "Tube  Type:  --,0D-;Itt 

16041 

1608  BEEP 
1612  Bop-0 

1616  INPUT  "ENTER  BULK  OIL  X  1OEFAULT-0X)  -,Bop 
1620  OUTPUT  CFi lei ; Bop 

1624  PRINT  USIN6  "  16X, ’"Bulk  o i IX- - - ,DO- . Bop 
1628' 

1632  BEEP 

1633  Ipo-l 

1634  INPUT  -ENTER  POOL  HEIGHT  ABOUE  TOP  TUBE  10-LES5  THEN  Sen,  l*Scn  OR  6REATER 
( DEFAULT >-,Ipo 

1635  OUTPUT  0Fi lei ; loo 

1636  PRINT  USING  -16X,--Pool  height--*  .DO";  Ipo 

1637  BEEP 

16391  NHT-NUMBER  OF  HEATED  TUBES 
1640  Nht-S 

1644  INPUT  "Enter  nunber  of  heated  inatrunented  tubes! default-5  )" ,Nht 
1648  OUTPUT  SFiIel.Nht 

1652  PRINT  USING  "  16X ,* "Nunber  of  heated  mstrunented  tubes-" * ,DD‘ .Nht 

1653  BEEP 

1655  Ipos-1 

1656  IF  Nht-1  THEN  INPUT  "WHICH  POSITION  IS  THIS  SINGLE  TUBE  (1  TO  5  0EC- 1  ) * , I 
pos 

1657  Ipos*Ipos-l 

1658  BEEP 
I860' 
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166-  Natb-humber  of  active  dummy  pairs 
loot  Natp-0 

1671  !N°UT  "Enter  number  of  active  dummy  pairs  i  De  f  aul  t  -0  I  ■ ,  Na  t  c 
I67E  OUTPUT  GFilel.Nato 

1680  PRINT  USIN6  '  16X "Number  of  active  dummy  oa ir s» ■ ■ ,0C" . Na tc 
168-*  BEEP 
1668  > 

169:'  NRT-NUHBER  op  added  heated  TUBES  to  ENHANCE  bundle  EcFECT 
1696  N-t-0 

1700  INPUT  "Enter  number  of  added  heated  tube:  from  simulation  hea t ers ( De f aul t » 
0  ) " ,  Nr  t 

1704  OUTPUT  #Fi lei ;Nrt 

1708  PRINT  US1N6  ’ 16X ," "Number  of  added  heated  tubeslfrom  simulation  heaters!*" 

",00'iNrt 

1712  BEEP 

1716' 

1720'  CORR  IS  CORRECTION  FOR  INSTRUMENTED  TUBE  HEI6HT 
1724  Corr*0 

1728  INPUT  "UANT  TO  CORRECT  TSAT  FOR  TUBE  HEI6HT  (  0-YES(  DEFAULT ). 1-NO >" ,Corr 
1732  IF  Corr-0  THEN  PRINT  USIN6  "16X,“"T5AT  is  corrected  instrumented  heat 

ed  tube  heioht""" 

1736  IF  Corr-1  THEN  PRINT  USIN6  "16X,""TSAT  is  NOT  corrected  for  instrumen 

ted  heated  tube  heiQht""" 

1740  OUTPUT  CFi lei ; Corr 
1744  BEEP 

1748'  ILQU-1NPUT  MODE-  LIQUID,  UAPOR.OR  LIQUID  VAPOR  AVERAGE 
1752  Ilqv-0 

1756  INPUT  "SELECT  < 0-LIQ1 default  ) , 1 -VAF , 2-( LIQ+VAP  >/2  ) " , I Iqv 
17601 

1764'  OlA-Oiameter  at  thermocouple  positions  (meters) 

1768  DATA  . 0122 , 0. 0098 ,0. 01 0£ , 0. 01 16 , 0 . 0 . 0 
1772  READ  Diet*) 

1776  Dl-Dlatltt) 

1 780 ' 

1784'  D2*Diameter  to  base  of  fins  (outside  die  for  smooth Xneters > 

1788  DATA  .  0158 ,0 . 0125 , 0 . 0 156 .0 . 01 4 15 , 0 .0 ,0 
1792  READ  D2a( •  ) 

179E  0:*D2e: Itt  > 

1800' 

1804'  Di-Inside  diameter  of  unenhanced  ends  (meters) 

1806  DATA  .0132,0. 0109, e. 0116,0. 0127, 0,e,0 
1812  READ  Diet •> 

1816  Di*0ia( Itt  ) 

1820' 

1824'  Do-Outside  diameter  of  unenhanced  ends  (meters) 

ie:e  data  .015675,0.0125,0.015875,0.01415,0,0,0 

1 632  READ  Doa( •  ) 

1636  Oo*Doe! Itt) 

1840' 

1844'  L“Lenoth  of  enhanced  surface  (meters) 

1848  DATA  .2032, .2032, .2032, .2032, .2032. .2032. .2032 
ieS2  READLa'*) 
ieSE  L-La(Itt) 

I860' 

1864'  Lu-CORRECTEC  Length  of  unenhanced  surface  at  the  end:  tMETEFS 

1863'LU»LCIN  *  THICKNESS/2 

1672  DATA  .0261, .0254, .0264. C. 0256, C, 0,0 

1 E7£  READ  Luat •  ) 

I860  Lu-Lue- Itt  ) 

18  =  4  ' 
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1 88  S  ■  L'.— corrected  length  o'  2  inch  finneC  lile  enC 
1692  01"  wve : 6  i 

I S9E  OrTA  . 0769 , . 0762 , . C?"Z ,0 . 076S .0 .0 .0 
1900  READ  Lva'.  •  ) 

1904  Lv-LvaUtt) 

1906  Kcua-Thermai  Conductivity  of  tube 

191C  data  40i ,c.c,e.e,c.c 

1916 '  READ  Kcua  •  ' 

1920'  Kcu-Rcua' 1 1 1  > 

1924  A-P I • ( 0c ‘2-Di ‘21/4 

1928  F-PI*0o 

1932  J-l 

1936  Sx-0 

1940  Sy-0 

1944  SxS-0 

1946  Sxy-0 

1952  Repeat :  ' 

1956' 

1960  IF  Im-0  THEN 

1964'  Otld-desired  temperature  of  liquid 

1968'  Dtld-47.5  'R-113 

1969  Dtld-2.2  'R-114 

1972  Ido*2 

1976  ON  KEV  0,15  RECOVER  1952 

1960  PRINTER  IS  1 

1994  PRINT  USING  *4X  ,  *  "SELECT  OPTION 

1988  PRINT  USING  "6X.*"0-TAKE  OATA""" 

1992  PRINT  USING  "6X,""1-SET  HEAT  FLUX" 

1996  PRINT  USING  "6X,""2-SET  Taat  (DEFAULT  SET  FOR  R-114)""" 

2000  PRINT  USING  "4X , " "NOTE :  KEY  0  *  ESCAPE" "" 

2004'  Ido-desired  option 

2008  BEEP 

2-12  INPUT  Ido 

2016' 

2020  BEEP 

2024'  Set  default  value  for  input 

2028  Ic  Ido>2  Then  Ido-2 

2632'  Tafe  data  option 

2235  IF  Ido-0  THEN  2440 

2740' 

2,. 44'  LOOP  TO  SET  HEAT  flux  (FOR  TOP  INSTRUHENTED  TUBE) 

2348  IF  Ido-1  THEN 

2e52  Dq dp -10000C 

2l56  PRINT  USING  "4X,""0dp  QDPsim  Nrt 

Qtot""" 

2P50  PRINT  USING  "4X , " " ( V/n"2 >  (W/m"2) 

<U  )""• 


Qdpaux 

(U/m'2  > 


2.64  Err- 1 

2066'  Reset.read  channel  25-30 , automat ic  scaling 

2(  32 i  Channel  25*au>  amps,26*sir  anpi,27*mst  volts, 28-sim  volts, 29-au» 

vol ts , 30-34* ins t  amps 

2-76  OUTPUT  709. -AR  AF25  Al34  VR5" 

2P30  FOR  1-10  TO  11 

2094  0UT°UT  70S  "AS  SA" 

2036  ENTER  709. Amp' I  , 

2092  NEXT  I 

2056  FOR  1-0  TC  2 

;iee  output  705 -as  sa- 

2104  Enter  709. Volt ( I  ’ 

210E  next  1 
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2ies 

2112 

rue 


2120  1 

2124 

:i:e 

Ziz: 

:ize 

:i40 
2144 
2146 
2152 
2155 
2160 
2154 
2168 
2172 
2176 
2180 
2184 
2188 
2192 
(0  J.Qdp 
2196 
2200 
2204 
2208 
2212' 
2216' 
2220 
2224 
2228 
2232! 
2233 
2236 
2240 
2244 
2246 
2252 
2256 
2260 
2264 

Tva 

2268 
2272 ' 
2276 
2280' 
rnocoup 
2284 

228e 

2292 
2296 
2300 
2304 
2306 
2212 
2316 
232  e 
2324' 


606  i-C  TO  4 
0L'T°UT  709.  ‘AS  SA" 

ENTEF  70S. Anp< I  ) 

NEXT  1 

Calculate  actual  heat  flu* 

0(0  >«SP*Uol  t  <  0  >*Anp<  lpos  > 

QoplC  )-O(0  )/<PI*02*L  ) 

Os  im-S0*20*Uo  1 1  <  1  '*Amp<  11  ) 

Qdps in-Qs in/  (FI*C2*.2032*3  1 
Oau»-60* 20*Uo 1 t ( 2  >*Anp< 10) 

Qdpau* -Qaux/ < PI*. 0160*. 1778*4  ) 

Qtot-010 )*Nht*0sin+0au« 

Nrt*QdpSin/Qdp< 0  ) 

IF  ABSi Aqdp-Dqdo >>Err  THEN 
IF  Aqdp>Dqdp  THEN 
BEEP  4000.. 2 

ELSE 

BEEP  250.. 2 
END  IF 

IF  Nr t < . 1  THEN  Nrt-0 
IF  Qdpaux<100  THEN  Qdpaux-0 
IF  Odps in< 1 00  THEN  Qdpsin-0 

PRINT  USING  -AX.ZtMZ.ZDE.DXl.DX.lMDD.DDl^X.ZtHZ.SDE.DXl^Odp 
sin.Nrt .Odpaux.Otot 
WAIT  2 
60T0  2076 
END  IF 
END  IF 

LOOP  TO  SET  Tsat 
IF  ldo-2  THEN 

IF  Ikdt-1  THEN  2240 
BEEP 

INPUT  'ENTER  DESIRED  Tsat  ( 0EFAULT-47.S  C  -  R-lljC.Dtld 
INPUT  -ENTER  OESIRED  Tsat  (DEFAULT-2.2  C  -  R-114>-,Dtld 
Hdt»l 
01 dl-0 
01 d2-0 
Nn-1 

Nrs*Nn  MOD  15 
Nn»Nn* 1 
IF  Nrs-1  THEN 

PRINT  USING  " 4X , " *  DT sat  Tldl  Tld2  TlbC 

Tvab  Tlav 

END  IF 

Read  thernocouple  voltages  for  vapor,  liquid 
OUTPUT  70S.-RF  AF0  AL5  UR5 " 

Sanple  each  thermocouple  20  tines  and  report  temp  for  each  the 
le,  vapor-0, 1,2;  liquid*3&4 
FOR  1-0  TO  E 
Sun-0 

OUTPUT  70S. "AS  SA¬ 
POR  Ji-1  TO  20 

ENTER  709 . E 1  iq 
Sun-Sun+El iq 
NEXT  Ji 
£m  f  <  I  '-Sun  '  20 
T  < I  >*CNT  vs  v ( Enf  < I  ’  ■ 

NEXT  I 

Compute  average  temperature  cf  liquid 
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;-;£ 

T  1  a v- '  T  .  2  '*T'  4  >•  .5 

-  ~~2  * 

Conoute  everagt  tenpe-etu'"e  of  vepo' 

zzzz 

Tva*  T<0  ><T<  :  >  >/2 

2334 

Tve«2*T(2 ) 

2336 

Tvav-i  T<0  )*T(  1  )  +  T<  2  >  i-  2 

zz*e 

IF  ABS( Tlav-Dt Id  >  ’  .2  THEN 

C  34  4 

IF  T 1 a v : D 1 1 C  THEN 

C348 

8EEP  4000, .2 

;35C 

ELSE 

2356 

8EEP  25e..2 

2360 

END  IF 

2364 

ELSE 

2368 

IF  ABS! Tlav-Dt Id )>. 1  THEN 

kj  r*. 

IF  At  ld>Dt  1  d  THEN 

2376 

BEEP  3000,. 2 

2380 

ELSE 

2384 

BEEP  800, .2 

2388 

END  IF 

2392 

END  IF 

2396 

END  IF 

2400 

Errl*T lav-Oldl 

2404 

Oldl-Tlav 

2408 

£rr2-T  vav-01d2 

2412 

01d2-Tvav 

2416 

PRINT  US1N6  "4X ,7( MOOD . DO ,3X  ) ' . Dt 1 d , T  < 

V 

2420 

WAIT  2 

2424 

GOTO  2252 

2428  END 

IF 

2432' 

2436!  TAKE  DATA  IF  In-0  LOOP 

2440  IF 

Ikol-1  THEN  2452 

2444 

BEEF 

2448 

Ikol-1 

2452 

OUTPUT  709. *AR  AF0  AL5  UR5‘ 

2456 

FOR  1-0  TO  5 

2460 

OUTPUT  709. "AS  SA" 

2464 

Sun-0 

2468 

FOR  Ji-1  TO  20 

2472 

ENTER  709 • E 

2476 

Sun-Sun+E 

2480 

IF  I >2  THEN  Et! Ji-1  )-E 

2484 

NEXT  Ji 

2488 

Kdl-0 

2492 

IF  I >2  THEN 

2496 

Eave-Sun/20 

2500 

Sun-0 . 

2504 

FOP  Jk-0  TO  19 

250e 

IF  ABSfEtf  JI  >-Eave  K5.0E-I 

2512 

5un-Sun+Et<  Jk ) 

2516 

ELSE 

2520 

Kdl-Kdl+1 

2524 

END  IF 

2328 

NEXT  J ! 

2532 

1^  I  2  THEN  FRINt  USING  ■4>.," 

2535 

IF  !  dl  10  then 

2S40 

EEEP 

2544 

EEEP 

2548 

PRINT  USING  1 4 X  ,  " " T o c 

oeat  oete  set 

.... 

■K  d  1 


.DC 


.  kdl 
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:ss: 

:5S6 

:sBe 

2564 

2568 

2s?: 

2S7B 

2580 

2584 

2586 

2592 

2595 

2600 

2B04 

2608 

2612' 

2616' 

2620 

2624 

2628 

2632 

2636 

2640' 

2644' 

2648 

2652 

2656 

2660 

2664 

2668' 

2672 

2676 

2680 

2684 

2688 

2692 

2696 

2700 

2704' 

2708' 

2712 

2716 

2720 

2724 

2728 

2732 

2736 

2740 

2744' 

2748 

2749 

2750 

2751 

2752 

2754 

2755 

2756 

2757 
27S5 
2760 


GOTO  1580 
6  NO  16 
END  IF 

Emf :  1  '-Sum.-  v  20-Kdl  1 
NEXT  1 

OUTPUT  709. "AR  AF4C  AL69  UPS* 

FOR  1-6  TO  35 

OUTPUT  709.  *A£  SA* 

Sum-0 

FOR  Ji-1  TO  5 
ENTER  709.  E 
Sum-Sum+E 
NEXT  Ji 
£mf !  1  )-Sum/5 
NEXT  I 

READ  U0LTA6EE  (  27- Inst  ,  28-Sin,  29-Au>.  ) 

OUTPUT  709. *AR  AF27  AL2S  UR5* 

FOR  1-0  TO  2 

OUTPUT  709; *AS  SA* 

ENTER  709; Volt! 1 > 

NEXT  I 

READ  CURRENTS  (30-34-Inst  tubes ; 35-39-ACT IUE  Dummy) 
OUTPUT  709. *AR  AF30  AL39  URS* 

FOR  1-0  TO  9 

OUTPUT  709. *AS  SA* 

ENTER  709. Amp ( I ) 

NEXT  I 

Read  Currents! 25-Au«  amps,26-Sim  amps) 

OUTPUT  709; *AR  AF25  AL2E  VR5* 

FOR  1-10  TO  11 

OUTPUT  709; *AS  SA* 

ENTER  709. Amp! I  ) 

NEXT  I 
ELSE 

ENTER  0C 1 1 e2 ; Em  f  < •  ) , Vo  1 1 !  •  ), Amp ( •  ) 

ENO  IF 

CONUERT  EMr'S  TO  TEMP , VOLT , CURRENT 
FOR  1-0  TO  35 

T ( I  i-FNT vsv! Emf < I )) 

IF  I >5  AND  ldtc>0  THEN 
FOR  I i-0  TO  Idte-1 

IF  Ldtc!  Ii  )-I-5+39  THEN  T(I>— 99.99 
NEXT  I. 

END  IF 
NEXT  I 

Ntc-nr  of  thermocouples 
Ntc-6 

IF  lpos>0  THEN 

Q(  Ipos  >-50»Vol  t  (  0  )»Amp(  Ipos  > 

Twe' Ipos  ;-0 


Ndtc-8 

FOR  1-1  TO  Ntc 

Np»Ipos‘6*o*J  j 
J  j  ■  J  ,i  ♦  1 

IF  A8S( T( Nn  )  ;  99  THEN 
T(Nn>— 99.99 
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Ndtc-Ndt  e  + 1 
ELSE 

Twa(  Ipos  >»Twai  Ipos  >  +  T<  Nn 
END  IF 
NEXT  I 

Tui a(  Ipos  >«Tua<  Ipos  )/(6-Ndtc  > 

goto  :s:e 

END  IF 

FOR  Ii-e  TC  A 

O'  Ii  >-60+Volt(0  >+Anp<  li  ) 

Tua-Averafie  temperature  of  the  wall 

Tuia(  1 1  )-0 

Ndtc-0 

FOR  I-I  TO  Ntc 

Nn  is  counter  in  temp  array,  start  at  E  (this  is  the  first  th 


ernocouple  in  the  tube  bent  ) 


Nn-Ii«G+l+S 
IF  ABSl  T ( Nn ) )>99  THEN 
T<Nn>— 99.99 
Ndte-Ndtc+1 

ELSE 

Tua(  1 1  )-Tua<  Ii  )+T<Nn  ) 
END  IF 
NEXT  I 

Twa(  I  i  >“T*ia<  1 1  >/  ( 6-Ndtc  > 

NEXT  Ii 

Tlav-(T(3)+T(4))/2 
Twav«T ( 2  ) 

Tvav-(  T<0  >+T<  I  >+T<2  )  )/3 


Tlav*T ( S ) 
Teu*T«e(0 1 
Kcu-FNKcu<  Tcu > 


'THERMAL  CONDUCTIVITY  OF  COPPER 

'IF  CURVE  FIT  NOT  AVAIL  USE  ARRAY  KCU( • > 


2844'  FOURIER  CONDUCTION  EOUATION  WITH  CONTACT  RESISTANCE  NEGLECTED 
2848  for  I-e  TO  4 

2852  Tg<  I  >»Tua<  I  !-0<  I  >*L0G<02/Dl  1/ (  2+PI +Kcu+L  > 

285G  IF  Ilqv0  THEN  Texs-Tlav 

28E0  IF  Ilqv-1  THEN  Texs-Tvav 

28G4  IF  I  lqv»2  THEN  Texs-<  T1  ev+T<  2  >  >•  .  5 

2868  IF  Corr-1  THEN  Thetab( I  )»Tw< I >-Texs 

2872  IF  Corr-0  THEN  Thetabf I  )»Tg( I  >-(Texs+. 056+1*. 123)  'R-114 

2876'  IF  Corr-0  AND  Ipo-1  THEN  The taDU  >-Tw<  1  >-<  Te*s+ . 054+1  *.  144  >  'R-ll 

3 

2877'  IF  Corr-0  AND  Ipo-0  THEN  The t e&( I >-Tw( I  )-< Texs- 1 . 078+ . 147. 1)  'P- 


NEXT  I 


COMPUTE  VARIOUS  PROPERTIES 

Tf  iln-(Tw<0  >+Te»s  )♦  .5  'FILM  TEMPERATURE 

Rho-FNRho' Tf 1 1 n  )  'DENSITY 


Mu-FNMu(Tfiln) 
y  -FNK  (  Tf  i  Ip.  ) 
Cp-FNCp ( Tf i In ) 

Be  te*FN6e ta( Tf i In ) 
Ni-Mu/Rhc 
A lpha-K / ( Rho+Cp > 
Pr-Hi /Alpha 


'DENSITY 

'VISCOSITY 

'Thermal  conductivity 
'SPECIFIC  HEAT 

'Thermal  expansion 

'►  INEMATIC  VISCOSITY 

'Thermal  ciffusivity 
1  prandtl 


COMPUTE  natural-convective  heat-transfer  COEFFICIENT 
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2936' 

:9*e 

2944 

:94s 

:ss; 

2956 

2960 

2964 

2966 

29?2 

2976 

2980 

2984 

2988' 


FOP  UNENHANCEC  ENO<  5  : 

Lu*tus >  1 1 1  ' 

Hber-190 

Fe*  (  Hbar *P  '  ( h cu* A  > ) *  .5*Lu 
Tanh-FNTenh( Fe 1 

Thetai Idos  l-Thetab(Ioos >«Tanh/Fe 

X**(S.6L*6ete«Thetab'  loos  '•Gc*3«Tenh'(Fe-Ni-Alpha  ;  -  _  i£6SE"’ 
Vy-<  l  +  <  .S59/Pr  f(  9.'16  >  )*(  8/27  ) 

Hbarc-K/Go+i  .6-* .  387+Xx/Vy  >'2 
IF  AB5i  (Hba<--Hbare  )/Hbarc  )'  .00!  THEM 
Hbar-<  Hbar+Hberc  1+ . 5 
60TC  2948 
END  IF 


2992 1  COMPUTE  HEAT  LOSS  RATE  THR0U6H  UNENHANCED  END! S  ) 

2996  01(0 '» (The  tab'.  Ipo5)*Tanh)»((Hba'-«P»Kcu*A;‘.S) 

3000  Qq-Ql(0>+Oq 

3004  2-2+1 

3006  IF  2-1  THEN 

3012  Lu-Lv 

3016  60T0  2944 

3020  END  IF 

3024  2-0 

3028  01pct-0q/0( loos  ) 

3032  Oq-0 

3036  As-PI +02+L 

3040  FOR  li-0  TO  4 

30*<  01 < I i >-Qlpct +Q< I i  > 

30*8  Odp< li >«(Q< lx )-01( li  )  )/As 

3052  Htube( h  )-Qdp( It )/Thetab( It  ) 

3056  NEXT  1j 

3060  PRINTER  IS  701 


3064' 

3068' 

RECORD  TIME  OF  DATA  TAKING 

3072 

IF  Im-0  THEN 

3076 

OUTPUT  709 ; -TO* 

3080 

ENTER  ?0S;Told* 

3084 

END  IF 

30es> 

CHURCHIIL/CHU  CORRELATION 

3087 

Raa-9. 81 +Bet a«Thetab(  Ipos ) 

3088 

Denom-( 1 .♦(  . 559/Pr  )“( 9/16  ! 

3089 

Nueh»< . 6+ . 387+ ( Raa/Denon ! ' 

3090 

Och*K  +  Nuch«Thetab(  Ipos )/<D 

30911 

3e92' 

OUTPUT  DATA  TO  PRINTER 

3096 

PRINTER  IS  702 

3100 

PRINT 

3104 

PRINT  USING  " 10X , " ’Da t a  Se 

3108 

PRINT 

3112 

PRINT  USING  -10X,"  Tvl 

Tvav 

TIdav 

3116 

PRINT  USING  ■lOX.etMDD.OO,; 

3120 

PRINT 

3124 

PrINT  USING  "6>: ,  "Tube  1 

H 

Thetat  ‘  ’  ‘ 

3128 

PRINT  USING  ”6X , *  *t  1 

2)  ( U/n 

"  2  .  K  )  ( K  > - 

3129 

IF  lpos'O  then 

3130 

jj-e 

3132 

FOP  Ji-0  TO  5 

Number  -  "  .DDD.2X , 14A' . J  ,Toldt 
Tv2  Tv3  Tidl  Tld2 
.'*.T<0>,T(  j  i,T(2  ),T(3),T(4),T(S) 
11  Temperatures  (Deg  C<  Tr.a 

:  4  S  E  :Des 


Tld3 

■ Tvav , T 1 ev 

ve  Odp 

C  J  (  ui r  “ 
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3133  *C' Ji  '-Ti ltos*E»J j +6  1 

3134  .'j-Jj  +  l 

3136  NEXT  Ji 

3137  Tnn-lpos+1 

3138  PRINT  USING  "EX  ,D,  IX  ,7(  HDD. DO  1 , 1X,2<  MZ .30E  ) .  IX  .  1  <  HOD. DO  ' .  Tnr. .  Tp<  C 
,  Tp!  1  ),Tp(3  ),Tp(  3  i,Tp<  4  )  ,Tp(  5  1 ,  Tu«(  Ipos  l,Odp(  Ipos  :■  , Htubel  Ipos  ,  The  tab-,  loos  > 

3139  60TC  3177 

3140  ENC  IF 

3141  Jj-0 

3143  FOP  1 l *0  TO  Nht-1 

3143  FOP  Ji-0  TO  S 

3144  Tp< J i  )*T( 1 1 *S+J j+6  ) 

3148  Jj-Jj+1 

3152  NEXT  Ji 

3156  Jj*n  +  i 

3160  FOR  Ji-e  TO  4 

3164  TnlJi >*1+Ji 

3168  NEXT  Ji 

3172  PRINT  USING  *6X ,0. IX ,?( nOD.OO  >, IX ,2<  HZ .30E  ) , IX, 1 1 MOD. 00  )' . Tn(  1 1 > , T 

p(0),Tp(l),Tp(2),Tp(3),Tp<4),Tp(5),Tua(Ii  )  ,Qdp<  1 1  )  ,Htube( 1 1  > , Thetabt 1 l  ) 

3176  NEXT  Ii 

3177  PRINT 

3179  PRINT  USIN6  *6X,"  Heat  Flux  and  Tdel  fron  Churchi 1 1 /Chu  Correlation  i 
s  ‘ ,  1  <  HZ  .  3DE  1 ,2X  ,  1<  fIDD.OO  )  ‘ ; Qch , The  tab ( Ipos  > 

3180  PRINT 

3182  OX-l 

3184  IF  In-0  THEN 

3188  BEEP 

3192  INPUT  ‘OK  TO  STORE  THIS  ORTA  SET  ( 1-YC default  ), 0-N >7 ‘ , Ok 

3196  END  IF 

32001  J-the  counter  Tor  data  sets 

3204  IF  OX-1  OR  In-1  THEN  J-J+l 

3208  IF  ox-1  AND  In-0  THEN  OUTPUT  #F i lei ;Enf < • > , Vol t ( • >, Anp< • 1 

3212  IF  In- 1  OR  Ox-1  THEN  OUTPUT  #Plot , Qdp< • ) ,Htube< • > , Thetao< • > 


3216  Go  on-1 

3220  IF  In-0  THEN 

3224  BEEP 

3228  INPUT  "WILL  THERE  BE  ANOTHER  RUN  ( 1*Y< default  /  , 0-N  l9 " ,Gc_on 

3232  Nrun-J 


3236  IF  6o_on-G  THEN  3272 

3240  IF  Go_on<>0  THEN  Repeat 

3244  ELSE 

3248  IF  J<Nrun+l  THEN  Repeat 

3252  END  IF 

3256  St-1 

3260  BEEP 

3264  IN°UT  "ARE  YOU  SURE  YOUR  READY  TO  TERMINATE  t 1  —  Y  < DEFAULT  i , 0-NO  )7 ‘ , S t 

3268  Go_on-l 

3272  IF  St>e  THEN  3280 

3276  IF  St-e  THEN  GOTO  3240 

3280  IF  In-0  THEN 

3284  BEEP 

3288  PRINT 

3292  PRINT  USING  ‘  10X , '  'NOTE  data  runs  were  stored  ir  T;le 

, 10A‘ . J-l  ,D2f ; let 

329E  ASSIGN  @Fi iel  TO  • 

3300  OJTPUT  0File2.N-un-l 

3304  ASSIGN  ®ri lei  TO  Dlfile* 

3306  ENTER  ®c  i  lei  .  Datet  ,Ldtc<  •  ; ,  1 1 1  ,6od  ,Nnt  ,Natc  .  Nr*  ,Co>‘' 

3312  OUTPUT  0Fi le2 . Datet .Ldtci  • ! , 1 1 1 .Bop ,Nht ,Natp ,Nrt , 6c" 
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3216  FOF  1-1  TC  Nrur-1 

3320  £NT£F  ffFilel.E-tT'.'.Vclt'./.Mmp  ♦ 

3324  OUTP'JT  S-  1 1  e2  ;  Enf  (  •  / ,  Vcl  t  '•  •  < . «np'  •  • 

3226  NEXT  I 

3332  ASS16N  0Fjlel  TG  • 

3336  PURGE  'OUMMv 

3340  END  IP 

3344  BEER 

334e  PRINT 

3352  PRINT  USING  '10X,"'NOTE  ",Z2.‘"  X-V  pairs  ue+e  sto-ed  in  plot  oata  f 

lie  " " , 10A' ; J-l ,Pf i let 

3356  ASSIGN  0File2  TO  • 

3360  ASSIGN  ©Plot  TO  ♦ 

3364  BEEP 

3366  SU8EN0 
3372' 

3376 'XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX XXX xxxxxxxxxxxxxxxxxxxxxxx 
3380' 

3384  OEF  FNKcu(Tcu) 

338B '  OFHC  COPPER 

3392  Tk-Tcu+273 . IS  !C  TO  K 

3396  Kcu-434-.ll2*Tk  '2S0-300K  USE  FOR  R-114  02.2  C 

3400!  Kcu-433 . 0- . 1 *Tk  '200-400K  USE  FOR  R-113  047.5  C 

3404  RETURN  Kcu 
3408  FNEND 
3412! 

3416  DEF  FNhu(T) 

3420!  CURUE  FIT  OF  VISCOSITY 
3424  Tk-T+273. 15  'C  TO  K 

3428  Mu“EXP<  -4 . 4636+<  1011 . 47/Tk  )  )•  l  . 0£-3  'R-114  170-360  K 

3432'  Hu- .0000134»(  1 0 - < 503 /<  Tk-2. IS ) ) )  'R113 

3436  RETURN  Mu 

3440  FNEND 

34441 

3448  DEF  FNCp(T) 

3462'  CURUE  FIT  OF  Cp 

3456  Tk-T+272. 15  'C  TC  k 

3460  Cp-. 40188+1. 6S007E-3*Tk+l.S1494E-6*Tk*:-6.678S3E-ie*Tk'3  'R-114  180-400  K 

3464'  Cp-<  929*1 .C3*T )«. 001  'R-113 

3468  Cp-Cp*1000 

3472  RETURN  Cp 

3476  FNEND 

3480' 

3484  OEF  FNRho(T) 

3488  Tv-T+273. 15  'C  TO  K 

3492  X-I-( l.8*Tk/753.95)  'K  TC  R 

3496  Ro-36. 32+61 . 146414* X*( 1/3 >♦ 16 . 41 80 15* X+ 17. 476836.X' .5+1 . 119828*X'2 

3500  Ro-Ro/ .062428  'R-114 

3504'  Rc*1.6207479E+3-T.(2.2186346+T.2.357B291E-3>  'R-113 

3506  RETURN  Ro 
3512  FNEND 
3516' 

2520  DEF  FNPr(T)  'GOOD  FOF  R-114/R-113 
3524  Pr-FNCp(T  i.FNMu'T  >/FNM  T  ! 

3528  RETURN  cr 
3532  FNEND 
3536' 

3540  DEF  FNK(T) 

3544 i  T <360  y  UITh  T  IN  C 
2546  K- . 071 - . 00026 1 *T 
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355:  RETURN  t 
.556  FNEND 
3560  1 

3564  DEF  FNTanhlFei 
3566  P-E>P(Fe) 

3573  O-EXP' -F«: 

3576  Tann-<  P-0  ).' :  P*Q  ) 

3580  RETURN  Tanr 
3584  FNEND 
3588' 

353:  OEF  FN'Tvsv(U) 

3596  COM  /Cc/  Cl  ?! 

3600  T-C(0> 

3604  FOP  1-1  TO  7 
3608  T-T+Cl I  )*V‘ 2 
3613  NEXT  I 
3616  RETURN  T 
3620  FNENO 
3624  1 

3628  DEF  FNBeta(T) 

3632  Rop-FNRhol T+ .  1  > 

3636  P.oo-FNRhol  T-.  1  ) 

3640  6eta--2/( Rop+Ron  >• (  Rop-Rom  )  /  .  2 

3644  RETURN  Beta 

3648  FNEND 

36S2  DEF  FNPolyl X  ) 

3656  COM  /Cply/  A< 10 , 10 ) , C< 10  )  ,B< 4  )  ,Nop , Iprnt ,0po , I Ioq 
3660  Xl-X 
3664  Poly-B<0) 

3668  FOR  1-1  TO  Nop 

3672  IF  IIoq-1  THEN  Xl-LOS(X) 

3676  Poly-Poly+R'.  I  )*X1 '  I 

3680  NEXT  I 

3684  IF  IIoq-1  THEN  Poly-EXP(Poly  ) 

3688  RETURN  Poly 
3692  FNEND 
3696' 

3700'XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

3704' 

3708  SUE  Roly* Df l leSI •  ) , Np , I tn  ) 

3712  DIM  R<  10), S<  lei.Syl  12  J,S*<  12  ),X«(  100  ),Yy(  1001, Xy(  17  > 

3716  COM  /Cply/  A( 10 . 10 ) ,C( 10  )  ,B< 4 ) ,N , Iprnt , Opo, I  loo 
3720  COM  /Xxyy/  Xp(5),Vp(5) 

3724  FOR  1-0  TO  4 
3728  B<I>*0 

3732  NEXT  I 
3736  Ir.»l 
3740  BEEP 

3744  INPUT  -ENTER  DATiR  FILE  NAME  ’  ,0' 1 1  et <  0  > 

3746  BEEP 

3752  INPUT  "ENTER  NUMBER  OF  X-v  PAIRS", Np 
3756  BEER 

3760  INPUT  "LIKE  TO  EXCLUDE  DATA  FAIRS  < I -X , 0-N : DEFAULT ! . 1 eC 
3764  IF  Ied-1  THEN 
3768  BEER 

2772  INRUT  "ENTEF  NUMBER  OF  PAIRS  TO  EE  EXCLUDED' . Ipe' 

3776  END  IF 

3780  ASSIGN  8r : 1 e  TC  Cfile*-0> 

"784  N»2 

278S  BEER 
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379:  INPUT  -ENTEF  THE  ORDER  Of  PC-YNOM1A-  'DE^AUlT-D'  ‘.N 

--’9E  FOP.  1*0  TO  N 

3800  Syl  I  >*0 

3804  S»!l**0 

3806  NEXT  I 

361?  IF  Ied-1  AND  l»'l  THEN 
3816  FOR  1*1  TO  Ipe> 

3630  ENTER  •File.Xyt-’ 

3834  NEXT  I 

3636  END  IF 

3633  FOR  1-1  TO  Np-Ipex 
3836  ENTER  0F  i le ; Xy( •  ) 

3840  IF  Opo-0  THEN 

3944  Y»Xyt Itn-1 > 

3845  X-Xy<ll  +  Itn> 

3e53  END  IF 

3856  IF  Opo-1  THEN 

3860  Y-Xy(5+Itn> 

3864  X*Xy( 11+Itn) 

3868  END  IF 

3872  IF  Opo*2  THEN 

3876  Y*Xy(S+Itn) 

3880  X-Xy< Itn-1) 

3884  END  IF 

3888  IF  1 1 oa*l  THEN 

3892  X-L06!  X  ) 

3896  Y-L06IY) 

3900  END  IF 

3904  Xx(I)»X 

3908  Yy(  I  )*Y 

3912  R<0)-Y 

3916  Sy(0)»Sy<0)+Y 

3920  S(l)-X 

3924  Sx< 1  >-Sx< 1  )+X 

3928  FOR  J-l  TO  N 

3932  R<  J  )-R(  J-l  )*X 

3936  Sy<  J  >-Sy<  J  )+P' J  > 

3940  NEXT  J 

3944  FOR  J*2  TO  N*3 

3946  S'.  J  >-S<  J-l  )*X 

3952  Sx(  J  >-Sx<  J  >+S< J  ) 

3956  NEXT  J 

3960  NEXT  I 
3964  Sx!0)-Np 

3966  FOR  1-0  TO  N 
3972  Cli-Syd) 

3576  FOR  J-0  TO  N 

3980  A!  I ,  J  >-S*(  I  +  J  1 

3984  NEXT  J 

Z98E  NEXT  I 

2992  FOR  1-0  TO  N-l 

3996  CALL  Divide! I  > 

4000  CALL  Subtract! !  +  !  ! 

4004  NEXT  I 

4006  6! N  >»C<  N  )/A<  N,N  ) 

4012  FOP  1-0  TO  N-l 

4016  B! N'l-I  »— Cf  N- 1 - J  / 

4C20  FOP  J-0  to  I 

4524  6' N'l-I  i-B'N-:-l  '-A' N-; 'I  .N-J  l-6'N-J  - 

4026  NEXT  J 
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403:  b< n- i - :  '-BtN-i-i i.  a,n-i-;,n-i-i  > 

4C3S  NEXT  I 

4040  IF  Iprnt-C  THEN 

4044  PRINT  USING  “13>  , --EVFCNENT  COEFFICIENT"" 

404E  FOP  I-C  TO  N 

4053  PRINT  USINC-  "  1SX' ,  DC ,  SX  ,  MD .  70E" .  I  ,B<  I  ) 

405E  NEXT  I 

4080  PRINT  "  ' 

4064  PRINT  USING  -12X."DATA  POINT  >  X  v:  CALCULATED 01 

SCREPANCY "  " 

4066  FOR  1-1  TO  No 

4073  Yc-B(0) 

4076  FOR  J-l  TO  N 

4080  Yc*Yc+8<  J  )*X>:(  I  )'J 

4084  NEXT  J 

40e6  D*Yy(I)-Yc 

4093  PRINT  USIN6  ‘ 15X . 3D ,4X ,4( MO .50E , IX  ) " ; I , Xx< I  ) .  Yy< 1  ! ,  Yc .0 

4096  NEXT  I 

4100  END  IF 

4104  ASSIGN  CFile  TO  • 

4108  SUBEND 
4112! 

4116  SUE  Divide(M) 

4120  COM  /Cply/  A! 10,10 >,C<10),B<4 i.N.Iprnt .Opo.Ilofl 
4124  FOR  1*M  TO  N 
4128  Ao-A(I,M) 

4132  FOR  J-M  TO  N 

4136  A(I,J)-A(I,J)/Ao 

4140  NEXT  J 

4144  C(I)-C(I)/Ap 

4148  NEXT  I 

4152  SUBEND 

4156' 

4160  SUE  Subtract ( K  ) 

4164  COM  /Cply/  A< 10,10 ),C(  10 ),B(4 i.N.Iprnt .Opo.Iloo 

4168  FOR  I-K  TO  N 

4172  FOR  J-K-l  TO  N 

4176  A(  I .  J  >-A(K-l,J  >-A(I,J  ) 

4180  NEXT  J 

4184  C(  1  )-Ct K-l  )-C(  1  ) 

4188  NEXT  I 
4192  SUBEND 
4196' 

4200  SUB  PI  in 

4204  COM  /Cply/  A(  10, 10  ),C(  10  ),B(4  ),N,  Iprnt  ,0po.  Hob 
4208  COM  /Xryy/  Xx(5),Yy<5> 

4212  PRINTER  IS  70S 
4216  BEEP 

4220  INPUT  "SELECT  (0-h/h0X  sane  tube.  l-h<  HF  )/h(  sr.  >- ,  Irt 
4224  BEEP 

4328  INPUT  "WHICH  Tsat  (  1-6 . 7,0— 2 .2  .  Iset 

4232  Xr,  i  n-0 

4236  X  _/-ie 

4240  X»tep»2 

4244  ip  Irt-0  THEN 

4248  Yr*;n“0 

4252  Yrnax- 1 . 4 

4256  Ystep*.2 

426e  ELSE 

4  264  Yr.  in*0 
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4  2  6  E  Y  na . *  1 5 

42~2  Ystep-5 

4276  END  IF 
4280  6EEP 

4264  print  *  in. spi . if  23ce.2200.s30e ,680e- " 

4268  PRINT  'SC  0 , 1 00 , e , 100 . Tl  2,0." 

4252  S f* ■ 1 00/ t X«a> -Xmi n  / 

4296  S ?y“100  Yma* - Ymir  > 

4300  PRINT  "PU  0,0  PD" 

4304  FOR  Xa“Xmn  TO  Xna>  STEP  xstep 

4308  X«(  Xe-Xm n  )«Sf  x 

4312  PRINT  -PA".X,",0,  XT." 

4316  NEXT  Xa 

4320  PRINT  "PA  100,0;PU; " 

4324  PRINT  "PU  PA  0,0  PO" 

4326  FOR  Ya*Ym n  TO  Yna-  STEP  Ystep 
4332  Y«( Ya-Ynin  )»Sfy 

4336  PRINT  -PA  0 ,  “ ;  Y , " YT" 

4340  NEXT  Ya 

4344  PRINT  *PA  0,100  Tl  0  2" 

4348  FOR  Xa-Xnin  TO  Xmax  STEP  Xstep 
4352  X»t  Xa-Xnin  >»Sf x 

4356  PRINT  "PA" ; X , * , 100 ;  XT- 

4360  NEXT  Xa 

4364  PRINT  -PA  100,100  PU  PA  100,0  PD‘ 

4368  FOR  Ya*Y«in  TO  Yma*  STEP  Ystep 

4372  Y»< Ya-Ynin )*Sfy 

4376  PRINT  "PO  PA  100,*,Y,-YT" 

4380  NEXT  Ya 

4384  PRINT  -PA  100,100  PU" 

4388  PRINT  "PA  0,-2  SR  1.5,2" 

4392  FOR  Xa-Xnin  TO  Xnax  STEP  Xstep 
4396  X»(  Xe-Xmm  >*Sfx 

4400  print  -PA- ; X , * ,0. • 

4404  PRINT  -CP  -2,-1 .LB-iXa;  ” 

4408  NEXT  Xa 

4412  PRINT  "PU  PA  0,0- 

4416  FOR  Ye»Yr,in  TO  Yna*  STEP  Ysteo 

4420  IF  AES t  Ya  )<1. E-5  THEN  Ya-0 

4424  Y»! Ye-Ynin  >*Sf y 

4428  PRINT  "PA  0,- .Y,-- 

4432  PRINT  ‘CP  -4,-.2S;LB"; Ya. ** 

4436  NEXT  Ya 

4440  Xlabel**"0i  1  Percent* 

4444  IF  Irt-0  THEN 
444e  Ylabel»*-h/h0S- 

4452  ELSE 

4456  Ylabelf * "h/nsnoot h " 

4460  END  IF 

4464  PRINT  ‘SF  1.5,2;PU  PA  50,-10  CP '. -l£N( X j abe li  )/2 . ‘0 . LP" . X  label  * 
4468  PRINT  -pa  -11,50  CF  0 , " • -LEN( Y labels  )/2«5/6.*DI  0 , 1 ; lE " ; Y1 abe  1$ 
4472  PRINT  "CP  e,0" 

447£  lpn»0 
448C  eEEc 

44f4  INpu:  -UANT  T0  p.L07  Dat f.  pR0~  4  cji_E  i  ,0-N  )T-  ,Oi  p 

A  Aft  l  -  r 

44=2  IF  Ol P*1  THEN 
449=  BEEP 

4500  INPUT  "ENTEP  THE  NAME  OF  THE  DATA  =ItE',t:  f|]e» 

4504  SEEP 
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4S0t  1N*>VT  -SELECT  '0-.INEAF.  1-_0G'X,\  r.IIo; 
i £::  ASSIST  ®F;le  TC  C_file* 

4S1E  6EEF 

4=:e  INPUT  "ENTEP  THE  BEGINNING  RUN  NUMBER' , MB 
45C4  BEEF 

45*5  INPUT  "ENTER  THE  NUMBER  OF  X-Y  PAIRS  STORED" .Noair 
4532  BEEP 

453c  INPU*  "ENT£R  DESIRED  HEAT  FLUX".G 

4540  BEEF 

4544  PRINTER  IS  1 

4546  PRINT  USING  "4X, ""Select  a  symbol 

4553  PRINT  USING  "4X."1  Star  3  Plui  sign*" 

4556  PRINT  USING  "4X,"3  Circle  4  Souare"" 

456e  PRINT  USING  "4X,""5  Rombus'" 

4564  PRINT  USING  "4X.""6  Right-side-up  triangle"* 

4568  PRINT  USING  *4X,*'7  Up-s 1 de-doun  triangle"" 

4572  INPUT  Syn 

4576  PRINTER  IS  70S 
4580  PRINT  'PU  DI' 

4584  IF  Sym-1  THEN  PRINT  *SM‘* 

4588  IF  Sym-2  THEN  PRINT  'SM+ “ 

4592  IF  Sym-3  THEN  PRINT  'SMo' 

4598  Nn-4 

4600  IF  Ilog-1  THEN  Nn-1 

4804  IF  Hd>!  THEN 

4608  FOR  1-1  TO  (Md-1  ) 

4812  ENTER  9File.Xa,Ya 

4616  NEXT  I 

4620  END  IF 
4624  Q1 “0 

4628  IF  Ilog-1  THEN  Q«L06<Q> 

4632  FOR  1-1  TO  Npairs 
4636  ENTER  (File.Xa.Bl*  > 

4640  Ya-E<0> 

4644  FOR  K- 1  TO  Nn 

4648  Ya*Ya+B<  K  )»2"K 

4652  NEXT  K 

4656  IF  Iiog-I  THEN  Ye-EXP(Ye) 

4660  IF  1 log-B  THEN  Ya-Bl/Ya 

4664  IF  Irt-0  THEN 

4668  IF  Xa-0  THEN 

4672  Yc-Ya 

4676  Ya-1 

4680  ELSE 

4684  Ya-Ya/Yo 

468c  END  IF 

46S2  ELSE 

4696  Hsn-FNHsnoct M C ■ Xa , isat  ) 

4700  Ya-Ya/Hsm 

4704  END  Ic 

4708  Xx< I-I  )-Xa 
4712  Yy !  1-1  )-Ya 
4716  X*<  Xa-Xr.m  )»Sf  > 

4”20  Y«'  Ye-Yr.ir  Y  •  S  F  > 

4724  IF  Syr,. 3  THEN  PRINT  "SM" 

4~3c  IF  Sy"  4  THEN  PRINT  "  Sc  1.4, 2. 4" 

4722  PRINT  -PA".Y,Y," 

4'36  I-  Syr.  3  THEN  PRINT  "SF  1.2. 1.6' 

4740  IF  Syn-4  THEN  CP!NT  "UC2. 4.95,0. -6, *4,0,0. 6.4,0. 
4744  ;f  5 yr-S  THEN  print  "U23 .0 , 95 , -3 .-£ .-2, £ ,2,E , 2 . *6 
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<-4e 

<?  s: 

4756 

4768 

4764 

4766 

4772 

4776 

4780 

4784 

4788 

4792 

4796 

4808 

4804 

4ee8 
4812 
4816 
4620 
4824 
4828 
4832 
4836 
4840 
4844 
4848 
4852 
48S6 
4860 
4864 
4868 
4872 
4876 
4880 
4884 
4888 
4eS2 
4895 
4900 
4  904 
4908 
4912 
4916 
4920 
4924 
4928 
4932 
4936 
4940 
4944 
4948 
4952 
4  956 
4  960 

4  964 
496e 
4  972 
4976 
4980 
4984 


16  Syl—E  THEN  PPINT  - uCC , 5 . 2 . 95 , 3 , -6  . -6 . 0 . 2 , 6  .  ~ 

IP  Sy"-7  THEN  PRINT  -UCC.-S.Z.99.-3.6.t,e.-2.-6. ' 

NEXT  I 

BEEP 

ASSI6N  SPile  TO  • 

END  IF 

PRINT  -PU  SM" 

SEEP 

INPUT  -WANT  TO  PLOT  ft  POLYNOMIAL  <  1-Y  ,0-N;'’  - .01  c 
IF  Ol.p*  1  then 
BEEP 

INPUT  "SELECT  ( 0-L INEAfi ,  1 -L06<  X , Y  ) ) "  ,  1 1 oq 
Ipr-nt*  1 

CALL  Folyl I tn  ) 

FOR  Xe-Xmir.  TO  Xma>  ST£C  Xstep/2S 
lcn-lcn+1 
Ya-FNPplyl Xa ) 

Y-(  Ye-Ymin  )«Sf  y 
X-(  Xa-Xnin  )«Sf* 

IF  Y <0  THEN  Y-e 

IF  Y> 100  THEN  GOTO  4868 

Pu-0 

IF  Ipn«l  THEN  Idf-Icn  MOD  2 
IF  Ipn-2  THEN  Idf-Icn  MOD  4 
IF  Ipn-3  THEN  Idf-Icn  MOD  8 

IF  Ipn-4  THEN  Idf-Icn  MOD  16 

IF  Ipn-5  THEN  Idf-Icn  MOO  32 

IF  Idf-1  THEN  Fu-1 

IF  Pu-0  THEN  PRINT  "PA-.X.Y.'PD" 

IF  Pu-1  THEN  PRINT  "PA" ,X.Y,"PU" 

NEXT  Xa 
PRINT  *PU" 
lpn-Ipn+1 
60T0  4480 
END  IF 
BEEP 

INPUT  "UANT  TO  QUIT  <  l-Y.O-NJP" .Iquit 
IF  Iqult*i  THEN  4904 
GOTO  4480 
PRINT  "PU  Se0" 

SUBEND 

SUB  Stats 

PRINTER  IS  701 

J-0 

K-0 

BEEP 

input  "plot  file  to  analyze’’ .File* 

ASSIGN  *rile  TG  Pile* 

BEEP 

INPUT  "LAST  RUN  No’' 0-OUIT  )" ,Nn 
IF  Nn-e  THEN  5092 
Nn-Nn- J 
S»-0 

sy-e 

Sr-0 

s«i-e 

Sv*»C 
S;  j-0 

POP  1-1  TO  Nr. 

J  - 1  *  1 
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493=  ENTgf  Srile.C.T 
49°:  M-Q/1 

499c  S.--S-  +  G 

SC00  S*s*S>s*C‘2 
5004  Sy*5y*T 
5003  Sys-Sys*T"2 
50i:  Sr*S:+H 
5016  Sis-Srs-i-H'S 
5020  NEXT  I 
5024  Qave-S.-./Nn 
5028  Tave*Sy/Nn 
5032  Have*S:/Nn 

5036  Sdevq-SORi ABS<  <Nn»S^.s-Sx"2)/(Nn*<  Nn-1  ) )  1  :■ 

5040  Sdevt-SOR>  ABS(<Nn«Sys-Sy*2  >/<  Nn«<  Nn-1 ) ) ) ) 

5044  Sdevh*SQR<  ABS<  <Nn«S:s-Sz'2>/<Nn*<  Nri-1  ) ) ) ) 

5048  Sh* 1 00*Sdevn/Have 
5052  Sq-100«Sdevq/Qave 
5056  5 t * i 00« Saev t /Tave 

5060  IF  K-l  THEN  5084 
5064  PRINT 

5068  PRINT  USING  ’IIX.-'DATA  FILE •"■,14ft" ;Fi  let 
5072  PRINT 

5076  PRINT  USING  "1 IX, ""RUN  Htube  SdevH  Qdp  SdevQ  Thetab  SdevT"" 

5080  K-l 

5084  PRINT  USING  ' 1 IX ,00 , 2< 2X ,0 .30E, IX ,30. 20 > , 2X ,0D. 3D , IX , 30 . 20" ; J .Have , Sh , Oave 

,  Sq  ,  Tave ,  St 

5088  GOTO  4940 

5092  ASSIGN  BFilel  TO  • 

5096  PRINTER  IS  1 
5100  SUBEND 
5104  SUB  Coet 

5108  COM  /Cply/  AU0.101.Ci  10) ,B< 4 ) ,N, Iprnt .Opo, Iloe 
5112  BEEP 

5116  INPUT  -GIUE  A  NAHE  FOR  CROSS-PLOT  FILE’.Cpfl 
5120  BEEP 

5124  INPUT  -OUTPUT  TYPE  ( 0-q  vs  Ot .  1-h  vs  Ot .  2-h  vs  qr.Opc 
5128  CREATE  BOAT  Cpft.E 
5132  ASSIGN  OFi le  TO  Cpft 
5136  BEEP 

5140  input  -SELECT  ( 0-L INEAR , l -LOGi  X , Y ) ) * , I  log 
5144  BEEP 

5148  INPUT  -ENTER  OIL  PERCENT  ( -1-STOP  1" , Bop 
5152  BEEP 

5156  INPUT  -ENTER  TUBE  NUMBER  (l,  2.  3,  4,  OR  5  >’  .  1 tn 
5160  IF  Bop  :0  THEN  5176 
5164  CALL  Polyiltn) 

5168  OUTPUT  0Fiie.6oc.E'.  •  ■ 

5172  GOTO  5144 

SI 76  ASSIGN  0=  :  le  TO  • 

5180  SUBENC 
5184  i 

5188'  xxxxxxxxxxxxxxxxxxxxYxxxxx/xxxxxxxxyxxxxxxxxxxxxxxyxyy.xxxxy 

5192' 

5196  SUE  Plot 

5-00  COM  /Cc  1  y  t  A  1 C .  1 C  ,0  1 0  . ,  E  ’  4  ; , Noc  ,  Ip--,'. ,  Op z  ,  1 1  cp 

5204  DIM  > :  17  i 

5203  INTEGER  Ii 

5212  PRINTER  IS  i 

5216  SEER 
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5220 

lavl 

5224 

INPUT 

"LIKE 

OEFAULT 

UALUES  fop  PLO"  1 i-Y< DEFAULT  -,0-N 

5226 

Opo-0 

5232 

BEEP 

5238 

PRINT 

USING 

■4X,""Select  Option  -*■ 

5240 

PRINT 

USIN6 

•EX,  --0 

Q  versus  delta-T< DEFAULT " 

5244 

PRINT 

USING 

*6x,--l 

h  versus  delta-T""" 

5246 

PRINT 

USING 

"6X , ‘  "2 

h  ve-sus  q " " " 

S252 

INPUT 

Opo 

5256 

BEEP 

5260 

INPUT 

"SELECT  UNITS 

<  0-SI > OEFAULT  ), I-ENGLISH  )• , Iun 

5264  PRINTER  IS  70S 

s:ee  ip  iavoi  then 

5:7:  BEEP 

5276  INPUT  -ENTER  NUMBER  OP  CYCLES  POR  X-AXIS". C» 

5280  BEEP 

5284  INPUT  -ENTER  NUMBER  OP  CYCLES  FOR  Y-AXIS-, Cy 
5288  BEEP 

5292  INPUT  -ENTER  MIN  X-UALUE  (MULTIPLE  OF  10)-.Xnin 

5296  BEEP 

5300  INPUT  -ENTER  MIN  Y-UALUE  (MULTIPLE  OP  10>-,Ynin 

5304  Else 

5308  IF  Opo-0  THEN 

5312  Cy-2 

5316  Cx-2 

5320  Xnin-1 

5324  Ymin-1000 

5328  ENO  IF 
5332  IF  Opo-1  THEN 
5336  Cy-2 

5340  Cx-2 

5344  Xxun-1 

5348  Y is  in- 100 

5252  ENC  IF 
5356  IF  Opo-2  THEN 
5360  Cy-2 

5364  C«-2 

5368  Xxun-1000 

5372  Yrun-100 

5376  ENC  IF 
5330  ENC  IF 
5384  EEEF 

5388  PRINT  *IN;SP1;IP  2300,2200.6300,6800; ’ 

5392  print  -sc  e.iee.e, 100. tl  2,0.* 

5396  Sfx-100/Cx 
5400  S»y-10B/Cy 
5404  BEEP 

5408  in°ut  -uant  to  by-pass  cage  u-y,  e-N0(0£FAULT>",jbp 
5412  IF  I bp* I  then  5908 
5416  PRINT  -Pu  0,0  PD- 
5420  Nn-9 

5424  FOP  1-1  TO  C* ♦ 1 
5428  Xet-Xnin«l0" :  I-l  > 

5432  IF  3  — C x ♦ 1  THEN  Nn-l 

543E  =0F  J-l  TC  Nr. 

5440  IF  .1-;  THEN  PF I NT  "TL  2  0" 

5444  IF  J-;  THEN  -tl  1  e- 

544  £  xe-Xat-J 

5432  > -L5T '  xa 'xr.  1  r.  >•£  f , 

54SE  RF IN"  -PA" . >  ,C .  >*.  " 
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546C  NEXT  ; 

5464  N£V*  ; 

5468  PE  1  NT  "PE  10C.C.PU; " 

54?:  print  -pu  pe  c.e  p d” 

5476  Nn-9 

548C  C0P  1-1  TO  Cy+ 1 
5464  Vat -Ymin*  10"<  I- 1  ) 

5486  Ir  I  -Cy* 1  THEN  Nn-1 

54$:  FOF  J-l  TO  Nn 

5498  IF  J-l  THEN  PRINT  "TL  2  C 

ssae  if  j-:  then  print  -tl  i  e- 

5504  Ye-Yat  •  J 

5508  Y-L6T(Ya/Ypun)«Sfy 

5512  PRINT  "PE  0." .Y.'YT" 

5516  NEXT  J 

5520  NEXT  I 

5524  PRINT  "PE  0,100  TL  C  2" 

5528  Nn-9 

5532  FOR  I-l  TO  C*+l 
5536  Xat-X«in»10'( I-l  ) 

5S40  IF  I-Ck+1  THEN  Nn-1 

5544  FOR  J-l  TO  Nn 

5548  IF  J-l  THEN  PRINT  " TL  0  2" 

5552  IF  J>1  THEN  PRINT  "TL  0  1" 

5556  Xa-Xat«J 

5560  X-L6T( Xa/Xnin  )*Sfx 

5564  PRINT  "PE" i X , * , 100 ;  XT- 

5568  NEXT  J 

5572  NEXT  I 

5576  PRINT  "PE  100,100  PU  PE  100,0  PO" 

5580  Nn-9 

5584  FOR  I-l  TO  Cy+1 

5588  Yat-Yrun- 10" i I-l  ) 

5592  IF  I-Cy+1  THEN  Nn-1 

5596  FOP  J-l  TO  Nn 

5600  IF  J-l  THEN  PRINT  "Tl  0  2" 

5604  IF  J,1  THEN  PRINT  "TL  0  1" 

5608  Ya-YaW 

5612  Y-LGT(Ya/Ynm»«Sfy 

5616  PRINT  "PD  PE  100, '.Y.'YT' 

5620  NEXT  J 

5524  NEXT  1 

5S2e  PRINT  "PE  100,100  PU" 

5632  PRINT  "PE  e,-2  SR  1.5,2" 

5636  I  i«i_GT<  Xru  n  ) 

5640  FOR  I-l  TC  C**l 
5644  Xa-Xn;nM0'<  I-l  > 

564e  X-LGT(  Xa/Xrun  )* S f  * 

5552  PRINT  "PE" ; X, ' ,0. ' 

5656  IF  1 1  >-0  THEN  PRINT  "CP  -2  , -2  .  L810  .  PF  -2 . 2  .  !_E  ‘ .  1 1  . 

5660  IF  1 1 <e  THEN  PRINT  "C°  -2 . -2 . LB  10 . PR  0 ,2 . LB"  • 1 1 . "  " 

5664  1 1 - 1 1  +  1 

5566  NEXT  I 

5672  PRINT  "Pu  FE  0,0" 

5676  I  i*_GT (  Ymr.  - 
5680  v 1 0* 1 G 
5564  F0-  1-1  TC  Cv4  1 

5686  Ye-Ynjn*10- < I-l  , 

5692  v-LGTt  Ya 'Yni n  »-5  T  y 

5696  PRINT  -PE  0,-.Y,-- 
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s~oc  print  -CF  -4,-.:s.lei0.ff  :le* . ; i . • - 

£704  I  1 ■ 1 1  +  1 

S~OS  NEXT  I 
5-::  BEEP 
5  ""16  Idl-1 

S-;e  INPUT  "WANT  USE  DEFAULT  LABELS  (1-Y<  default  ),0-N  )•>-,  Idl 
S?:4  IP  Idl  : M  then 


st:e  beep 

5732  IN°UT  ‘ENTER  X-LABEL" .Xiabel* 

S*3E  BEEP 

S~40  INPUT  -ENTER  Y -lABEL " , Vlobelt 

5714  END  IP 

5746  IF  Ooo<2  THEN 

5752  PRINT  ‘SR  1.2. PL1  PA  40,-14;- 

5756  PRINT  -LB(T,PR  -1.6,3  PD  PR  1.2,0  PU.PF  . 5 , -4 , LBuo ; PR  .5,1;- 
5760  PRINT  "LE-T ; PR  .5,-1 ; LBsat ; PR  .5.1;* 

5764  IF  Iun-0  THEN 

5768  PRINT  ‘LB  )  /  (K  >* 

5772  ELSE 

5776  PRINT  "LB  )  /  <  F  ) " 

5780  END  IF 

5784  END  IF 
5788  IF  Opo-2  THEN 
5732  IF  Iun-0  THEN 

5796  PRINT  "SR  1.5,2;PU  PA  40,-14;LBq  /  <U/n;SR  1,1. 5. PR  0 . 5 , 1 ;  LB2 ;  SP.  1 

. 5 , 2 ; PR  0. 5,-1. LB)“ 

5800  ELSE 

5804  PRINT  ‘SR  1.5,2-PU  PA  34,-14;L8q  /  <8tu/hr;PR  ,S,.5.LB.;PP  .5, -.5; 

• 

5808  PRINT  -L8ft;PR  ,5,1:SR  1,1.S;LB2;SR  1.5, 2. PR  ,5,-l;LB>.- 

sei2  END  IF 

5816  END  IF 

5820  IF  Opo-0  THEN 

5824  IF  Iun-0  THEN 

5828  PRINT  -SR  1 .5 ,2 ; PU  PA  -12,40;OI  0,l;L8q  /  <U/n;PR  -1,0.S;SR  1,1. 5;L 

82. SR  1 .5, 2 ; PR  1,.5;LB>- 
5632  ELSE 

5636  PRINT  -SR  l.S,2;PU  PA  -12,32.01  0,l;LBq  /  (Btu/hr.PR  -  .5 , . 5 . LB . .PR 

.5. .5. ■ 

5640  PRINT  'LBftiSR  1,1. 5. PR  -1,.5;LB2.PP  1,.5;SR  1 .5,2 ;LB  )* 

5844  END  IF 

5848  END  IF 
5852  IF  Opo>0  THEN 
5856  IF  Iun-0  THEN 

5860  PRINT  -SR  1 . 5 , 2 ; PU  PA  - 12 , 38 : DI  0,l;L8h  /  (U'n;PR  -1,.S;SR  1.1. 5. LB 

2 . SR  1 .5,2 .PP  .5. .5. • 

5864  PRINT  ’LB.. PR  . 5.0.L6K  )" 

5868  ELSE 

5872  PRINT  -SR  1 .5,2 ;pU  PA  -12,28.01  e.l.LBh  /  <etu/hr;PR  - .5, .5; LE . . PF 

.5,  .5;  - 

5676  PRINT  ”LBf t ; PR  -1..5.SR  1,1.5;LB2.SR  1.5.2;PR  ,5,.5;L6..PR  .5, .5. 

LB-  >" 

5e80  END  IF 

5884  END  IF 
5888  IF  Idl-e  THEN 

Se92  PRINT  -SF  1.S.2.PU  PA  5C.-16  CP"  .  -LEN  >.  label*  i/2  .  -0,  lE"  •>  lebel*  .  ‘  “ 

5696  PRINT  -PA  -14.SC  CP  C.-.-LEN, Ylabelf  ••2oE/E. “DI  0, 1 ;lB‘ • Ylabelt . "  ‘ 

5S00  PRINT  -cr  C,0  ci- 

5904  END  IF 
5908  Ion-0 
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SSi:  Repeat  > 

S91E  Xll-l.E+E 
5920  xui--:.E»6 
=924  Icn-e 
5928  BEEP 
59;:  04-1 

593E  INPUT  "WANT  TC  PLOT  DATA  FROM  A  FILE  (  1  -Y '  DEFAULT  ; , 0-N  >’ 1  , 0( 
S94e  IF  OH*!  THEN 
5944  BEEF 

5948  INPUT  "ENTER  THE  NAME  OF  THE  OftTft  FILE ” .Df 1 leS! 0 > 

5950  ASSIGN  @Fi la  TO  OfileStE) 

5956  8EEP 
5968  Npairs-20 

5964  INPUT  "ENTER  THE  NUMBER  OF  X-Y  PAIRS  STORED! DEFAULT-20  Npa ir 
5968  BEEP 
5972  Itn«Itn+l 

5976  INPUT  "ENTER  TUBE  NUMBER  <1,  2.  3,  4,  OP  5>".Itn 

5988  BEEP 

5984  PRINTER  IS  1 

5988  INPUT  "WANT  DEFAULT  SYMBOLS’  (YES-0  ( DEFAULT  )  ,N0-1 1  * ,  Synb 

5992  Sym-Itn+2 

5996  IF  Synb-0  THEN 

6000  60T0  6036 

6004  END  IF 

6008  PRINT  USIN6  "4X, " "Select  a  symbol  """ 

6012  PRINT  USIN6  “6X,““1  Star  2  Plus  sign"" 

6016  PRINT  USING  "6X,“"3  Circle  4  Square""" 

6028  PRINT  USING  "6X,""5  Rombua""" 

6024  PRINT  USING  "6X,’"6  Right-side-up  triangle""" 

6028  PRINT  USING  *6X , " "7  Up-s 1 de-doun  triangle""" 

6032  INPUT  Syn 

6036  PRINTER  IS  705 
6040  PRINT  "PU  DI" 

6044  IF  Sy«*l  THEN  PRINT  "SM-" 

6048  IF  Syn-2  THEN  PRINT  " SM+  * 

6052  IF  Syis-3  THEN  PRINT  "SMo" 

6056  FOR  1-1  TC  Npair- 
6060  ENTEF  ®F;ie.Xy(*  ) 

6064  IF  Opo-0  THEN 
6058  7 a - X y !  Itn-1  ) 

6072  Xa-Xy( 11+Itn  1 

6077  END  IF 
6080  IF  Opo-1  THEN 
6084  Ya-Xy( 5*!tn  ) 

6088  Xa-Xy( 11  +  Itn  ) 

6092  END  IF 

6096  IF  Opo-2  THEN 

6108  Ya-Xy!5+Itn> 

6104  Xa-X  y( Itn-1  ) 

6108  END  IF 

6112  IF  Xa<X 1 1  THEN  Xll-Xe 

6116  IF  Xa 'Xu  1  THEN  Xul-Xa 

612e  IF  Iun-1  THEN 

6124  IF  Opo  .2  THEN  Xe*Xe*l.E 

6126  IF  Opc-0  THEN  Ye*ie».17El 

6122  IF  Gpc*C  THEN  >e*<e-.217 

E13E  Ic  Opo-2  THEN  Xa->e*.  31- 

El  48  END  IF 

6144  > -LOT  ;  Xa  'Xm  1  n  >•  S  ’  ■ 

£  1 4 E  y.^5T' Ya  'Yr.ir.  J-Sf  y 
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6152 

j*0 

6156 

CALL  Syrtbl  >  .'i  ,Sve,,  Icl 

,L  j  ) 

6160 

60T0  6212 

6164 

IF  Sym> 2  THEN  PRINT  - 

SM" 

6168 

IF  Syw<4  THEN  FR1NT  ‘ 

SR  1.4 

6172 

IF  Icl-0  THEN 

6176 

PRINT  "PA' ,x,Y, "■ 

6180 

else 

6184 

PRINT  ‘FA’.X.Y.-FD 

Else  END  IF 

Bis:  IF  Sym;3  Then  PRINT  "SR  1. 2,1.6" 

6 1 96  IF  Sym-4  THEN  PRINT  ‘UC2 , 4 , 99,0  ,-6 , -4 . 0,0  .  e ,  4 , 0  ;  - 
6200  IF  Sym-S  THEN  PRINT  ‘UC3 . 8 . 99 , -3 , -6 , -2, 6 ,  3  ,6 , 3 . -6 . 
6204  IF  Sym-E  THEN  PRINT  ‘UC0 , S . 3 , 99 , 3 . -6 . -6 . 0 , 3 , 6 , " 
6208  IF  Sym-7  THEN  PRINT  ‘UC0 , -6 .3 , 99 , -3 , 8 .6 ,0  ,-3 . -8 . " 
6212  NEXT  I 
6216  PRINT  *PU’ 

6220  BEEP 


6224 

6228 

6232 

6236 

6240 

6244 

6248 

6252 

6256 

6260 

6264 

6268 

6272 

6276 

6280 

6284 

6288 

6292 

6296 

6300 


I lab-1 

INPUT  "WANT  TO  LABEL’  < 1-Y ( DEFAULT ), 0-N  ) ■ , I lab 
IF  I  lab-1  THEN 

print  ‘spo.sp:- 

BEEP 

IF  K lab-8  THEN 
Xlab-ES 
Ylab*8S 

INPUT  ‘ENTER  INITIAL  X.Y  LOCATIONS' ,Xlab , vlab 

Xtt-Xlab-5 

Ytt-Ylab+8 

PRINT  ‘SR  1,1.5* 

PRINT  ‘SM;PA‘,Xtt,Ytt,"LB  Tube  1  File' 
Ytt-Ytt-3 

PRINT  -PA-,Xtt,Ytt,*LB  No  Oil  Nane" 

IF  Syn«l  THEN  PRINT  ‘SM*‘ 

IF  Sym-2  THEN  PRINT  *SN+‘ 

IF  Syn-3  THEN  PRINT  ‘Sn0‘ 

Klab-I 
ENO  IF 


6304  K'j-1 

6308  CALL  Sy«b< X lat , Ylab , Sym , Ic 1 ,K j > 

6312  PRINT  ‘SR  1,1.5.SM‘ 

6316  IF  Syn<4  THEN  PRINT  ‘PR  2,0* 

6320  PRINT  ‘PR  2,0;LB* . Itn.  “ 

6324  BEEP 

6328  INPUT  ‘ENTER  BOP< 0-PCFAULT  )‘ , Bop 
6332  IF  Bop  10  THEN  PRINT  ‘PR  3 ,0 ;Lfi‘ ; Bop . ‘ ‘ 
6336  IF  Bop >9  THEN  PRINT  ‘PR  1 . 5 , 0 i LB ’ : Bop .  ‘ 
6340  PRINT  ‘PR  2 ,0 ; LB' i Of  1 le*<  0  ) ;  “ 

6344  PRINT  ‘S°0 . SP1 ; SR  1 .5,2" 

6346  Ylab-Ylab-S 

E352  END  IF 
6356  BEEP 

6360  ASSI6N  (File  TC  • 

E364  Xll-Xll/1.2 
E368  Xul -Vul *  1 . 2 
6Z~2  '  GOTO  8040 
£376  END  IF 
E380  PRINT  ‘PU  SN“ 

6364  BEEP 
6386  6c_on-l 
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E:  SI  INPUT  “WAN’  TO  PLOT  A  POLYNOMIAL  (  1-Y  '  DEPAULT  ,  ,0-N  :•»“  ,Gc_on 

6296  IF  6c_on-i  THEN 

6400  BEEF 

B404  PRINTEP  IS  1 

6406  INPUT  “WANT  DEFAULT  LINE  TYPE'?  , YES-0  <  DEFAULT  ),  NO- 1  Ln 

6411  lpn*Itn 

6416  IF  Ln-0  THEN 

6410  60TC-  6448 

6414  ENC  IF 

6426  PRINT  USING  “4X, "Select  line  type-'"" 

6431  PRINT  USINS  "BX , " "0  Solid  line"'" 

6436  PRINT  USING  *BX,""1  Dashed-" 

6440  PRINT  USING  *6X  , "  “I , ,  ,S  Longer  line  -  dash”" 

6444  INPUT  Ipn 

6448  PRINTER  IS  70S 

6452  BEEP 

6456  Ilog*l 

6460  INPUT  "SELECT  < 0-LIN. 1-L06< DEFAULT >>", 1  log 

6464  Iprnt-1 

6468  CALL  Poly( Of ileS< •  )  ,Npa irs , I tn  ) 

6472  FOR  Xx-0  TO  Cx  STEP  Cx/200 

6476  Xa*Xmin«10"Xx 

6480  IF  Xa<X  1 1  OR  Xa>Xul  THEN  6572 

6484  lcn»lcn+l 

6488  Pu-0 

6492  IF  Ipn-1  THEN  Idf-lcn  MOD  2 

6496  IF  Ipn-2  THEN  Idf-lcn  MOO  4 

6500  IF  Ipn-3  THEN  Idf-lcn  MOD  8 

6504  IF  Ipn-4  THEN  Idf-lcn  MOO  16 

6508  IF  Ipn-5  THEN  Idf-lcn  MOO  28 

6512  IF  Idf-1  THEN  Pu-1 

6516  Ya-FNPol y<  Xa  ) 

6520  IF  Ya<Ynm  THEN  6572 

6524  IF  Iun-1  THEN 

6528  IF  Opo<2  THEN  Xa-Xa*1.8 

6532  IF  Opo>0  THEN  Ya-Ya».1761 

6536  IF  Opo-0  THEN  Ya-Ya*.3l7 

6540  IF  Opo-2  THEN  Xa-Xa*.317 

6544  ENO  IF 

6548  Y-L6T< Ya/Ynin  ) • S f y 

6552  X-L6T  (Xa/Xmin)»Sf* 

6556  IF  Y <0  THEN  Y-0 

6560  IF  Y>10e  THEN  60TC  6572 

6564  IF  Pu-0  THEN  PRINT  "PA" , X , Y , "PO" 

6568  IF  Pu*l  THEN  PRINT  "PA" ,X, Y, "PU" 

E572  NEXT  X>. 

6576  PRINT  "PU" 

6580  END  IF 
6584  BEEP 

6588  INPUT  “WANT  TO  QUIT  ( 1 -Y ,0-N( DEFAULT >)", lot 

6592  IF  I q t  —  1  THEN  6600 

6596  GOTO  5916 

6600  PRINT  "PU  PA  0,0  SP0“ 

6604  SUBEND 
6606  1 

6612 1  xxxxxxxxyxxxxxxxxxxxxxxxxxxxxxxxxxx/xxxxxxxxxxxy xxxxxxxxxxxxxxx 
6616' 

6620  SUE  Synbf  X  ,  Y  ,  Svn ,  1  c  I  ,f  j  •’ 

6624  IF  Syr.--;  THEN  PRINT  "SM” 

6628  IF  Syn<4  THEN  PRINT  “SR  1.4, 2. 4“ 
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6632  Yad-0 

6636  I F  ► J-1  THEN  YaP- . 6 

6640  IF  Icl-0  THEN 

6644  PRINT  "PA" ,  X  ,  V  +  Y  a  d , 

664E  ELSE 

6652  PRINT  "PA" ,X ,Y+Yad, "PD" 

6E5E  END  IP 

6660  IF  Syr, >3  THEN  PRINT  *SF  1.2.1.E" 

6664  IF  Syr«4  THEN  PRINT  "UC2 . 4 . 9S . 0, -E . -4 . 0 . G , E . 4 ,0 : ' 

6666  IF  Syc-S  THEN  PRINT  "UC3 . 0 , 9S , -3 , -6 , -3 . E  . 3 .6 . 3 , -6 ; " 

6672  IF  Sym-E  THEN  PRINT  "UCB .5 .3 . 95, 3, -6 . -6 . 0. 3 . 8 . " 

6676  IF  Syn-7  THEN  PRINT  "UC0 ,-S. 3, 99, -3, 8. 6. 0.-3, -8.  " 

6680  IF  Kj-1  THEN  PRINT  "SM;PR  0,-.8" 

6684  SUBEND 
6688' 

6692'XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

6696' 

6708  SUE  Fixup 
6704'  FILE ■  FIXUP 
6708' 

6712  DIN  Enf <34  >,Amp< 11 >,Uolt<4 ),Ldtc(  4 ! 

6716  BEEP 

6720  INPUT  'OLD  FILE  TO  FIXUP" ,02f l leS 
6724  ASSI6N  8F i 1 e2  TO  02filet 
6728  D1 f l leS“" TEST " 

6732  CREATE  BDAT  Olfile*,60 
6736  ASSIGN  8FUel  TO  Olfile* 

6740  ENTER  8Fi 1*2 ; Nrun.DateS ,Ldtc< « ) , I t t , Bop ,Nht ,Natp ,Nrt ,Corr 
6744  OUTPUT  0Fi lei ; Nrun , OateS ,Ldtc(  •), It t ,6op ,Nht ,Natp ,Nrt , Corr 
6748  FOP  1-1  TO  Nrun 

6752  ENTER  8Fi Ie2 ; Told*, Emf( •) .Volt <♦) ,Anp( •  ) 

6756  IF  1-1  THEN  6764 

6760  OUTPUT  8F i  lei ; Bop , Tol dt , Enf  < • ),Uolt( • ), Amp< •  ) 

6764  NEXT  I 

6768  ASSI6N  8File2  TO  • 

6772  ASSIGN  8Fj lei  TO  • 

6776'  RENANE  "TEST"  TO  02_file* 

678e  BEEP  2000,-2 
6784  BEEP  4000,-2 
678B  BEEF  4000,-2 
6792  SUBEND 

6795' 

6800 'XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX XXX XXX 
6804  1 

6808  SUE  Move 

6612'  FILE  NAME  MOUE 

6816' 

6620  DIM  A<  66  ) , B<  66  !,C< 66  -  ,D< 66  !  ,E<  6£  > ,F'.  66  > ,G(  66  > ,H<  66  i ,  J'  66  :  ,M  66  1  ,L< 66  ;  ,M!  66 

) 

6824  DIM  N(  66  ),Enf  (34  :-,Uolt(2  i.Anpi  1 1  >,Ldtc' 4  ) 

6828  BEEF 

6832  IN°UT  "OLD  CILE  TC  MOUE " , C2_» : 1  el 
6636  ASSIGN  8Fjle2  TC  02_f;let 

6640  ENTER  8F 1 1  e2  .  Nrur.  .Toldi.Ldlc'  •  . , !  1 1 .  Bop  ,  NM'. .  Na  t  p  ,  N't ,  Cc-r 

6644  FOP  1-1  TC  Nrur 

6S4E  ENTER  @f , ] e: , Tol dt 

6652  ENTER  0F ■  ie2 .  A-  1  / ,Bi  I  > , C’  I  ■  ,D(  1  1  ,E(  I  :  ,F':  I  1 ,6(  I  ’  .hi  I  .  J '  I  ■  :  I  :  ,l '  I  ,M( 

i  ;  ,N'  i ) 

6956  ENTEF  0f i - Ewf ' •  i , Up  1 1 ' •  ; . Amp( •  ) 

6ese  nex t  ; 
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ES64  ASSIGN  ?:.ie:  TC  ♦ 

6666  EEEC 

66':  INPU*  "SHIFT  DISK  AND  HIT  CONTINUE-. 0* 

66'6  E£Ee 

668C  INPUT  -INPUT  80AT  SIZE", Size 
6864  CREATE  BOAT  02_f 1 1 e* . S 1  re 
6886  ASSIGN  0=Uel  TO  Dl.fileS 

689C  OUTPUT  @F ; lei ' Nrun, Da te* . Ldt c ! • ) . 1 1 1 , Bop . Nht . Na tc , Nr t . Co-r 

6896  FOR  1*1  TC  Nrun 

6900  OUTPUT  @p j lei ; Tel d* 

6904  OUTPUT  9Filel ,  At  I  >,Bi  I  >,C<  I  >.0<  I  ).E<  I  >.F(  I  ),6<  1  J.H!  I  ),  J<  I  ;.K(  I  ),L<  I  ),M 

(  I  >,N<  I) 

6908  OUTPUT  tf i lei ; Enf ( ♦  ) , Uol t < •  ) . Amp ( •  ) 

6912  NEXT  I 

6916  ASSIGN  @F 1 1  el  TO  • 

6923'  RENAME  ‘TEST-  TO  02_file* 

6924  BEEP  2000,-2 
6928  BEEP  4000,-2 
6932  BEEP  400e,.2 
6936  PRINT  -DATA  FILE  MOUEO' 

6940  SUBEND 
6944! 

6948'XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

6952' 

6956  SUB  Purg 
6960  BEEP 

6964  INPUT  -ENTER  FILE  NAME  TO  BE  DELETED* , F i le* 

6968  PURGE  File! 

6972  GOTO  6960 
6976  SUBENO 
6980' 

6984  IXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
6988' 

6992  sue  Conb 

6996'  FILE  NAME:  COMB 

7000' 

7004  0 1 M  Enf ( 34 ) , Uol t ( 2 1 ,  Anp< 1 1  ) ,Ldt c <  4  ) 

700e  BEEP 

7012  INPUT  -OLD  FILE  TC  FIXUP' , 02_f i le* 

7016  ASSI6N  *File2  TO  02_fiie* 

7020  Dl_file*--TEST- 

7024  CREATE  BOAT  Ol_flle*,30 

7028  ASSIGN  *Filel  TO  Dl.file* 

7032  ENTER  0File2.  Nrun,  Date*  ,Ldtc<  «  l.Iti  ,  Bop  ,  Nnt ,  Na  t  p ,  Nr  t ,  Com 

7036  IF  K *0  THEN  OUTPUT  *F  j  lei ;  Nrur. .Date* ,Ldtc<  «  ! ,  1 1 1 ,6oo  .Nht  ,Natp  ,Nr t  ,Corr 

7040  FOP  1-1  TC  Nrun 

7044  ENTER  9F;le2- Bop, Told*, Enf<« ;,Volt< •  l.Ampi • > 

7048  OUTPUT  ?F i lei ; Bop , Told* , Enf < •  ) , Vol t <  •  1 , Anp( «  ) 

70S2  NEXT  I 

705E  ASSIGN  @':le:  TO  • 

7060'  RENAME  "TEST"  TO  D2_fileS 
7064  BEEP  400e. .2 
7068  BEEP 
7072  Ql  a*  1 

7e78  Input  "UANt  TO  ADC  another  FILE  (  l-Y,0-N>  defaui  •-  ■■  ^  "  , 0»  6 
7080  IF  Ole-I  THEN 
7064  K-: 

708e  BEEP 

"’092  IN°UT  “GIUE  NEW  FILE  NAME  *  , Nf  i  1  e* 

7098  ASSIGN  0cile:  TO  Nfile* 
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7 IOC  GOTO  7832 

710-  END  IF 

7108  ASSIGN  ®File2  TO  • 

■’ll-  SUBEND 
71  IE  1 

7120IXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXVXXXXXXXXX 

7134’ 

7126  SUE  Readplot 

7132  DIM  Qdp(  5  >,Htube(S),ThetadiS> 

713E  PRINTER  IS  701 

7140  INPUT  ’ENTER  FILE  NAME’, File* 

7144  INPUT  ’ENTER  THE  NUMBER  OF  DATA  PAIRS’, Nrun 
7146  ASSI6N  JF i lei  TO  File* 

7152  FOR  1-1  TO  Nrun 

7156  ENTER  ®F i le 1 ; Odp ( •  ) , Htubel «  ) , The t ab< • > 

7180  PRINT  Qdp<«  ) 

7164  PRINT 

7166  PRINT  Htubel •  ) 

7172  PRINT 

7176  PRINT  Thetabl ♦  ) 

7180  PRINT 

7184  PRINT 

7188  NEXT  I 
7192  SUBEND 
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